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FOREWORD 


This  section  presents  a detailed  account  of  the  project  activities  m the  area  of 
cylinder  dynamic  testing.  A mathematical  analysis  of  locked-rod  or  impulse  testing 
is  developed  and  the  effects  of  changing  test  setup  parameters  are  also  demonstrated. 

A computer  program  for  simulating  locked-red  testing  is  presented  along  with  a program 
user's  guide.  Draft  procedures  for  both  stroking  (endurance)  and  locked-rod  (impulse) 
tests  are  included  in  an  appendix.  Results  of  experimentation  with  the  locked-rod 
test  setup  are  also  presented. 
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CHAPTER  I 
INTRODUCTION 

The  reliability  of  hydraulic  equipment  depends  to  a considerable  extent  on  the  proper 
selection  of  individual  components.  For  various  reasons,  most  of  which  were  indicated 
by  Hopler  at  the  29th  National  Conference  on  Fluid  Power,  Ref.  [1],  the  end-item  user 
should  exercise  significant  control  over  individual  component  specifications.  Such 
specifications  should  address  themselves  to  the  performance  of  the  component  rather 
than  design  detai I s . The  objective  of  the  U . S . Army  Mobility  Equipment  Research  and 
Development  Command  (MERADCOM)  - Fluid  Power  Research  Center  program  over  the 
past  five  years  has  been  directed  towards  implementing  this  philosophy  by  developing  and 
appraising  test  methods  and  procurement  specifications  for  hydraulic  components. 

Hydraulic  cylinders  are  important  components  on  a vast  majority  of  mobile  equipment 
and  function  not  only  as  fluid  power  components  but  also  as  structural  elements.  Recog- 
nizing the  importance  of  assessing  the  structural  integrity  of  hydraulic  cylinders  under 
dynamic  loading  the  U.S.  Army  MERADCOM  initiated  a study  of  two  currently  used  test 
procedures  in  1975.  Results  of  the  first  year's  effort  were  presented  in  the  Annual 
Report  on  MFRD/OSU  Hydraulic  System  Reliability  Program,  Ref.  [2]. 

The  current  project  was  undertaken  as  the  continuation  of  the  previous  year's 
effort  to  establish  a rational  basis  for  appraising  the  structural  integrity  of  hydraulic 
cylinders  under  static  and  dynamic  loading.  The  effort  has  been  done  as  a Technology 
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Development  Project  and  has  gained  the  active  participation  of  six  industrial  companies 
having  strong  interests  in  mobile  hydraulic  cylinders.  The  following  is  a list  of  the 
industrial  sponsors  and  their  monitors  for  the  Technology  Development  Project  on 
Cylinder  Structural  Integrity  Assessment. 


Allis  Chalmers 
Bruning  Division  of  Gould 
Deere  and  Company 
Eaton  Corporation 
J.  I.  Case  Company 
Koehring 


L.  Stikeleather 
K.  Koch 

D.  Malm 

R.  Lindgren 

E.  Falendysz 
J.  Parrett 


Since  some  of  the  monitors  were  assigned  primarily  for  administrative  control,  the 
sponsoring  companies  have  maintained  technical  liaison  through  persons  actively  involved 
in  cylinder  design  and  selection.  These  persons,  who  have  made  valuable  technical  in- 
puts to  the  program,  both  at  review  meetings  as  well  as  in  correspondence  include 
Dr.  G.  Ekstrom,  C.  Brundidge  and  E.  Esser  (Eaton) , M.  Beck  (Drott),  D.  Hancock 
(Allis-Chalmers) , W.  Snyder  (Deere  & Co.),  andJ.  Fairbairn  (Koehring).  The  U.S. 
Army  MERADCOM  was  represented  at  the  technical  meetings  by  H.  Y.  Smith. 

The  effort  has  been  divided  into  two  phases,  namely,  the  Static  Analysis  Phase  and 
the  Dynamic  Analysis  Phase.  The  results  of  the  Static  Analysis  Phase,  which  consists 
of  verification  fo  the  joint  flexure  theory  of  hydraulic  cylinder  deflection,  will  be  pre- 
sented in  a self-contained  report.  This  section  is  devoted  to  the  Dynamic  Analysis 
Phase  whose  objective  was  to  examine  the  effect  of  test  parameters  on  locked-rod  testing. 
The  next  chapter  presents  the  results  of  dynamic  stress  analysis  as  well  as  a parameter 
sensitivity  analysis.  Chapter  I I I includes  the  user's  guide  for  the  digital  computer 
program  for  simulating  the  locked-rod  test  in  addition  to  a listing  of  the  FORTRAN 
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source  code  for  the  program. 

Appendix  A contains  the  draft  procedures  for  stroking  and  locked-rod  testinq. 
These  procedures  incorporate  modification  to  reflect  t tie  results  of  the  parameter  sensi- 
tivity analysis  and  simulation  results. 


I 1-3 


CHAPTER  I I 


DYNAMIC  STRESS  ANALYSIS 


INTRODUCTION 


The  main  objective  of  this  effort  was  to  examine  the  effects  of  pin  clearances , frame 
stiffness,  and  flow  rate  on  stress  cycles  in  stroking  and  impulse  testing  of  cylinde-  s. 
Experimental  effort  was  to  be  devoted  to  measuring  dynamic  strain  and  p-  ssure  during 
cyclic  loading  on  a locked-rod  test  setup.  Based  on  the  analysis  and  review  of  experi- 
mental data,  revised  test  procedures  for  conducting  stroking  and  locked-rod  tests  were 
to  be  drafted.  The  following  sections  will  discuss  the  effect  of  test  frame  parameters  and 
present  experimental  data  on  locked-rod  testing.  General  conclusions  on  both  facets 
will  be  presented  at  the  end.  Appendices  A and  B contain  drafts  of  the  revised  proce- 
dures for  both  impulse  and  stroking  tests. 


EFFECT  OF  TEST  FRAME  PARAMETERS 


The  effect  of  test  frame  parameters  (e.g.  , frame  stiffness,  pin  clearances,  pur 
flow  rate)  all  can  be  best  explained  by  developing  mathematical  models  of  the  setups 
under  discussion.  The  locked-rod  setup  will  be  analyzed  below.  A parallel  analysis 
for  stroking  tests  was  furnished  in  the  report  on  the  previous  year's  effort,  Ref.  [2]. 


nmGEuma  page  not  fllmlu 
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The  following  symbols  will  be  used  in  the  various  equations  comprising  the  mathe- 
matical model.  Time  derivatives  of  quantities  will  be  indicated  by  a dot  over  the  symbol. 


P = System  pressure 

A = Cross  section  area  of  the  cylinder 

Q.  = Flow  rate  of  the  pump 

in 

V = Volume  of  the  test  cylinder  which  is  being  pressurized  during  one  half  of 
the  test  cycle 

x = Displacement  of  the  piston  with  respect  to  a given  reference  mark  on  the 
cyl  inder 

v = Velocity  of  the  piston  with  respect  to  the  cylinder 

I = Inertia  of  the  piston  and  attached  moving  parts 

K = Stiffness  of  the  test  frame 

B = Drag  coefficient  for  the  piston  rod 

CL  = Leakage  coefficient  for  the  cylinder 

P = Relief  valve  set  pressure 

r 

8 = Effective  bulk  modulus  of  system  fluid 


It  will  be  assumed  that  the  directional  control  valve  used  for  cycling  the  pressure 
switches  instantaneously  and  that  the  system  relief  valve  opens  when  the  set  pressure 
is  attained  and  insures  that  the  pressure  is  not  exceeded.  It  will  also  be  assured  that 
the  directional  control  valve  and  piping  offer  no  resistance  to  flow.  The  net  result  of 
these  assumptions  is  to  slightly  overestimate  the  pressure  rise  and  decay  rates  for  the 
cycle.  Though  the  assumptions  can  be  lifted  by  developing  more  elaborate  mathematical 
models,  such  a course  of  action  is  not  considered  necessary  for  the  majority  of  testing 
situat  ions . 
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It  is  sufficient  to  consider  only  half  of  a complete  cycle  since  the  analysis  is 
precisely  the  same  for  the  other  half.  The  half  cycle  under  consideration  consists  of 
a pressure  rise  period,  a dwell  period,  and  a pressure  decay  period  (i.e.  , the  wave- 
form is  approximately  trapezoidal). 

The  pressure  rise  is  not  instantaneous  for  three  reasons:  firstly,  the  fluid  and  the 
conductors  (i  .e.  , hoses,  tubing)  are  elastic  and  a finite  amount  of  fluid  is  needed  to 
make  up  for  the  compressibility  of  the  system  volume  as  well  as  the  dilation  of  the  hoses, 
tubing,  etc.;  secondly,  the  cylinder,  especially  if  it  has  worn-out  seals,  will  leak  at  a 
rate  generally  increasing  with  pressure  and  this  leakage  fluid  reduces  the  flow  available 
to  raise  the  system  pressure;  lastly,  because  of  the  flexure  and  slack  of  the  frame  the 
piston  rod  will  move  just  enough  to  take  up  flexure  and  slack  in  the  frame. 


The  following  equation  establishes  the  relationship  between  the  pump  flow  and  the 
various  quantities  referred  to  above: 


Q.  = Av  + V P + PC  (2-1) 

' 1 

The  first  term  on  the  right  hand  side  of  the  equation  is  the  flow  rate  needed  to  fill 
up  the  volume  rendered  void  by  movement  of  the  piston  rod.  The  second  term  expresses 
the  flow  rate  needed  to  raise  system  pressure  (i.e.,  it  is  the  "compressibility  flow  rate"), 
while  the  last  term  accounts  for  cylinder  leakage.  It  should  be  noted  that  in  this  equa- 
tion  Q.^  is  the  inflow  to  the  system  volume  and  is  obtained  by  subtracting  from  the  pump 
outflow  whatever  passes  over  the  relief  valve.  It  may  also  be  noted  that  once  the  re- 
lief valve  is  opened  and  the  frame  has  flexed  sufficiently  to  develop  the  necessary  reaction, 

. 

the  pressure  rise  rate,  P,  goes  to  zero  and  the  only  inflow  is  that  needed  to  handle 
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cylinder  leakage. 


The  second  equation  needed  to  describe  the  locked-rod  setup  is  a force  balance  equa- 
tion for  the  piston  rod  and  attached  moving  parts,  as  follows: 


PA  = Iv  + Bv  + Kx 


The  first  term  on  the  right  hand  side  is  the  force  needed  to  overcome  the  inertia 
of  the  moving  parts  (piston  + attachments,  if  any),  while  the  second  is  the  force  needed 
to  overcome  the  friction  drag  at  the  piston  and  gland.  The  last  term  is  the  reaction  of 
the  frame  due  to  flexure. 


PARAMETER  SENSITIVITY  ANALYSIS 


The  following  parameters  can  be  generally  changed  with  more  or  less  ease  for  any 

test  setup:  Q , P , I , B,  C,  , and  K.  Q.  can  be  changed  by  using  different  pumps , 

in  r L in 

operating  a pump  at  various  speeds  or  using  a variable  volume  pump.  P^  can  be  changed 
easily  if  an  adjustable  relief  valve  is  used  in  the  circuit.  The  last  four  parameters  cited 
above  are  characteristic  of  the  test  cylinder  and  the  test  frame.  Of  these,  the  leakage 
coefficient,  C^,  will  generally  increase  in  the  course  of  a test,  while  changes  in  the 
drag  coefficient , B,  depend  on  the  design  of  the  seals.  Frame  stiffness  is,  perhaps, 
the  most  difficult  to  change. 


It  is  instructive  to  examine  the  effect  of  changing  the  above  parameters  on  the  test 
cycle.  Even  though  the  effects  can  be  quantitized  by  simulating  the  system  using  different 
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corrbinations  of  parameter  values,  a better  understanding  of  the  influence  of  a specific 
parameter  is  obtained  by  examining  the  pertinent  equations.  Thus,  an  inspection  of 
Eq.  (2-1)  indicates  that: 


T 


(i)  A larger  pump  flow  rate  will  result  in  higher  pressure  rise  rates,  as  shown  in 
Fig.  2-1  . Too  low  a pump  flow  may  result  in  the  directional  control  valve 
switching  flow  before  the  relief  valve  setting  is  reached. 


Fig.  2-1  . Effect  of  Pump  Flow  Rate 
on  Pressure  Rise  Rate  and  Final 
Pressure 


(ii)  As  may  be  anticipated,  the  relief  valve  setting,  P^,  controls  the  pressure 
level,  and  consequently,  the  stress  level  to  which  the  test  cylinder  will  be 
subjected. 

(iii)  The  effect  of  piston  rod  inertia,  I , is  hard  to  ascertain  directly  from  Eqs.  (2-1) 
and  (2-2).  However,  by  differentiating  Eq.  (2-2)  with  respect  to  time  and  sub- 
stituting for  P and  P in  terms  of  x , v,  v,  and  v,  the  following  equation  is  ob- 
tained: 


Q. 


in 


(2-3) 

K 


If  v is  zero,  it  is  seen  that,  for  a constant,  Q.^,  any  increase  in  I will  result  in 
a reduction  of  v,  other  quantities  being  constant.  From  Eq.  (2-2)  it  can  be  seen  that 
a reduction  of  v implies  a lesser  value  of  P.  However,  Eq.  (2-1)  requires  that  a re- 
duction  in  P will  need  to  be  compensated  by  increases  in  either  v or  P or  both.  The 
overall  effect,  thus,  cannot  be  summarized  as  an  increase  or  decrease  in  the  rise  rate. 
Simulation  results  have  indicated  that  for  the  common  range  size  of  cylinders  and  flow 
rates,  piston  rod  inertia  has  a rather  small  impact  on  both  the  pressure  rise  rate  and 
the  final  pressure. 

(iv)  The  effect  of  the  cylinder  drag  is  equally  complex  and  hard  to  generalize  without 
numerical  values  for  other  parameters.  Simulation  using  selected  sizes  of 
cylinders  has  indicated  that  this  parameter  has  little  impact  on  the  waveform  for 
inpulse  testing. 

(v)  The  cylinder  leakage  coefficient,  C^,  is  an  important  parameter  since  it  affects 
both  the  pressure  rise  rate  and  the  final  pressure  attained  in  the  test  cycle. 
Inspection  of  Eq.  (2-1)  shows  that  an  increase  of  reduces  both  the  pressure 
rise  rate  as  well  as  the  final  pressure.  If  the  cylinder  leakage  is  excessive, 
all  the  pump  flow  may  leak  past  the  seal  and  none  over  the  system  relief  valve, 

(vi)  The  frame  stiffness,  K,  mainly  affects  the  pressure  rise  rate — a stiffer  frame 
leading  to  higher  rise  rates.  This  effect  is,  however,  not  very  significant  if 
the  test  frame  is  made  of  material  having  a high  Youngs  modulus  and  designed 
so  that  lateral  deflection  of  load-carrying  elements  is  minimized. 

The  effect  of  pin  clearances  can  be  ascertained  for  any  specified  test  setup  by  assign 
ing  a value  of  zero  to  the  frame  stiffness  until  the  slack  is  taken  up.  From  Eq.  (2-2), 
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it  can  be  seen  that  if  K is  zero,  a given  pressure  will  result  in  a higher  acceleration 
and  velocity  for  the  piston  rod.  Qualitatively,  the  pressure  rise  rate  will  be  small 
until  the  piston  rod  has  moved  far  enough  to  overcome  the  slack  due  to  pin  clearances, 
and  after  that  the  pressure  will  rise  exactly  as  before.  The  computer  program  written 
for  digital  simulation  and  discussed  in  Chapter  III,  includes  provisions  for  simulating  the 
behavior  of  a test  frame  with  variable  stiffness  and  a variety  of  pin  clearances.  Figure 
2-2  presents  the  effect  of  chagning  pin  clearances  on  the  pressure  rise  rate.  It  is  seen 
that  large  clearances  lead  to  lower  pressure  rise  rates,  while  very  small  clearances  can 
result  in  pressure  overshoots. 


TIME  — ► 

(secs ) 


Fig.  2-2.  Effect  of  Pin  Clearance  on  Pressure  Pise  During 
Locked-Rod  Test  (Simulated) 

* 


The  mathematical  model  presented  in  Eqs.  (2-1)  and  (2-2)  can  be  used  to  predict 
the  value  of  pressure  at  any  time  in  the  course  of  an  impulse  test  on  a cylinder  with  a 
rod  fixed  at  any  specified  position.  Using  this  value  of  pressure,  the  stresses  in  the 


cylinder,  rod  or  any  other  load  carrying  member  can  be  established  using  standard  for- 
mulas, Ref.  [3],  Thus,  for  example,  the  hoop  stress  at  the  inner  and  outer  surfaces  of 
the  cylinder  wall , in  terms  of  the  internal  pressure  and  the  cylinder  dimensions  are  given 
by: 


Su  A p 

Hoop  int . — 


(R2  + r2) 
(R2  - r2) 


and: 


Su  . A 

Hoop  ext . — 


2r 


(R‘ 


2. 

r ) 


(2-4) 


(2-5) 


where  R and  r are  the  external  and  internal  radii,  respectively,  of  the  cylinder. 

Since  the  cylinder  wall  thickness  is  almost  always  within  the  range  0.1"  to  1"  (2.5  mm 
to  25  mm)  the  time  required  for  a stress  wave  to  travel  from  the  internal  surface  to  the 
outer,  or  vice  versa  is  of  the  order  of  10  ^ seconds.  Since  pressure  rise  rates  in 
inpulse  tests  rarely  exceed  100,000  psi/sec  (690  MPa/s),  the  time  needed  for  pressure 
to  rise  from  0 to  3000  psi  (0-21  MPa)  is  0.03  seconds  which  is  orders  of  magnitude  larger 
than  the  time  needed  for  the  stress  wave  to  travel  from  the  inside  to  the  outside  of  the 
cylinder  or  vice  versa.  Consequently,  for  all  practical  purposes,  the  stress  rise  accom- 
panying the  pressure  rise  occurs  simultaneously.  Experimental  data  presented  subse- 
quently confirm  this  deduction. 


EXPERIMENTAL  EFFORT 


The  final  report  on  the  previous  year's  activity  in  the  area  of  cylinder  structural 


integrity  assessment  presented  the  results  of  dynamic  pressure  and  strain  measurements 
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on  a locked-rod  test  setup.  Ref.  [2].  It  showed  that  the  pressure  rise  rate  correlated 
well  with  strain  rate  in  the  cylinder  in  both  the  axial  and  hoop  directions.  Consequently, 
in  conducting  fatigue  tests  (either  impulse  or  stroking)  it  is  sufficient  to  measure  cylinder 
pressure  dynamically  and  to  calculate  the  hoop  and  longitudinal  stresses  using  standard 
formulas,  Ref.  [3]. 

Figure  2-3  presents  traces  of  cylinder  hoop  strain  and  pressure  for  one  half  of  a 
complete  impulse  cycle  in  a locked-rod  test.  The  strain  rise  is  seen  to  occur  synchro- 
nously with  the  pressure  rise,  which  affirms  the  validity  of  Eqs.  (2-4)  and  (2-5). 

Figure  2-4  presents  similar  traces  of  pressure  and  axial  strain.  Once  again  the  corre- 
lation between  pressure  rise  and  axial  strain  rise  is  seen  to  be  very  good,  thus,  obviating 
the  need  for  strain  measurements. 

It  should  be  noted  that  the  actual  rise  rate  and  the  final  pressure  for  a cycle  depend 
on  the  pump  flow  rate,  characteristics  of  the  system  relief  valve,  switching  speed  of  the 
directional  control  valve,  and  capacitance  (i.e. , elasticity)  of  the  tubing  hoses,  etc.  , 
connecting  the  test  cylinder  to  the  directional  control  valve  and  pump.  In  addition, 
it  also  depends  on  the  stiffness  of  the  test  frame — a stiffer  frame  resulting  in  higher 
pressure  rise  rates.  Table  2-1  summarizes  the  effect  of  changing  some  of  the  para- 
meters mentioned  above.  The  following  parameters  were  kept  constant: 

(i)  pump  flow  rate 

(ii)  test  frame  stiffness 

(iii)  cylinder 

(iv)  fluid  bulk  modulus 

(v)  relief  valve  gradient  (slope  of  pressure  flow  characteristics) 

(vi)  upstream  volume 
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Fig.  2-4.  Dynamic  Pressure  and  Axial  Strain  for  Locked-Rod  Test 


Since  the  locked-rod  test  setup  is  essentially  a first-order  dynamic  system,  its 
transient  behavior  can  be  adequately  described  by  its  time  constant — a larger  time  con- 
stant implying  a "spongy"  system  with  (ow  pressure  rise  rate  and  a small  time  constant 
a "stiff"  system.  Table  2-1  indicates  that  a 250%  change  in  the  slack  in  the  frame 
(caused  mainly  by  pin  clearances)  changes  the  time  constant  only  7%.  Change  in  the 


relief  valve  cracking  pressure,  and  consequently,  the  pressure  Eit  which  the  cycling  is 
performed,  on  the  other  hand,  significantly  affects  the  pressure  rise  rate.  Consequently, 
high  pressure  cylinders  tested  at  the  same  flow  rate  will  be  subjected  to  higher  pressure 
rise  rates  than  low  pressure  cylinders.  An  inspection  of  Pigs.  2-3  and  2-4  confirm 
this  conclusion.  It  is  seen  that  in  both  cases  the  time  taken  to  reach  final  pressure 
is  approx i mate1  y one  second.  Since  the  final  pressure  is  3000  psi  (20.691  MFa)  in  the 
first  case  and  1000  psi  (6.897  MPa)  in  the  second,  the  rise  rate  is  approximately  tripled 
due  to  the  increase  in  relief  valve  setting. 
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CHAPTER 


USER'S  GUIDE  FOR  LCKROD 


LOKROD  is  a digital  computer  program,  written  in  FORTRAN  IV,  for  simulating 
a locked-rod  hydraulic  cylinder  test.  It  takes  into  account  the  compressibility  of  the 
system  fluid,  inertia  and  drag  of  moving  parts,  leakage  in  the  test  cylinder  and  pin 
clearances  in  establishing  the  pressure  waveform.  The  program  can  be  used  to  examine 
the  influence  of  changing  not  only  the  above  parameters  but  also  the  hydraulic  test  circuit 
pump  and  system  relief  valve.  Mathematical  details  of  the  analysis  are  reported  in 
Chapter  I I . The  remainder  of  this  chapter  is  devoted  to  explaining  the  structure  and 
usage  of  the  computer  program , which  is  included  in  Appendix  C. 


PROGRAM  STRUCTURE 


The  locked-rod  hydraulic  cylinder  test  setup  is  treated  in  the  program  as  a lumped 
parameter  dynamic  system.  It  is  modeled  by  a third— order  vector  nonlinear  differential 
equation  with  a step  change  in  flow  (caused  by  switching  of  the  directional  valve)  being 
the  forcing  function.  Any  other  type  of  forcing  function  can  be  imposed  on  the  mathema- 
tical model  by  appropriate  changes  in  one  of  the  subroutines  to  described  subsequentK  . 
Since  the  system  is  nonlinear,  a numerical  integration  method  is  used  to  develop  the 
time-histories  of  the  dynamic  quantities  (i.e.,  system  pressure,  piston  rod  displacement, 
and  its  velocity).  The  numerical  integration  method  calculates  the  values  of  the  dynamic 
quantities  at  a specific  point  in  time  by  extrapolating  from  the  values  of  the  points  in 
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time  prior  to  the  time  under  consideration.  The  process  is  iterated  to  develop  the  com- 


plete time-history. 

The  time  step  between  two  consecutive  points  on  the  trajectories  of  the  dynamic 
variables  depends  not  only  on  the  dynamics  of  the  system  being  simulated  but  also  on  the 
numerical  integration  method.  Since  the  time  step,  which  may  be  variable,  influences 
the  speed  of  simulation,  it  is  important  to  match  the  integration  method  to  the  system  being 
simulated.  LOKROD  uses  Gear's  algorithm  for  numerical  integration,  a fairly  recent 
development  which  is  suited  for  the  system  under  consideration,  Ref.  [4],  In  order  to 
use  Gear's  algorithm,  which  is  included  in  a subroutine  called  DIFSUB,  it  is  necessary 
to  furnish  a main  program  as  well  as  certain  auxiliary  programs.  These  will  now  be 
discussed. 

MA  I N 

The  main  program  is  used  primarily  for  controlling  input  and  output  (I/O)  for  the 
computer.  All  the  system  parameters  as  well  as  the  integration  parameters  are  read 
off  data  cards  and  printed  out  by  the  main  program.  The  initialization  of  the  dynamic 
variables  is  also  performed  in  the  main  program  prior  to  calling  the  integration  subrou- 
tine DIFSUB.  Each  call  to  DIFSUB  advances  the  trajectories  of  the  dynamic  variables 
by  one  time  step.  Between  repeated  calls  to  time  DIFSUB,  the  main  program  check 
if  the  integration  is  proceeding  satisfactorily  and  stores  the  results  of  the  numerical 
integration  after  a specified  number  of  time  steps  have  been  taken.  Finally,  the  main 
program  prints  out  the  time-histories  of  the  dynamic  variables  and  any  messages  pertain- 
ing to  abnormal  termination  of  the  integration,  if  any  such  terminations  occur. 
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DIFFUN 


This  subroutine  contains  the  mathematical  model  for  the  system  in  the  form  pre- 
scribed by  Gear's  integration  subroutine  DIFSUB.  The  dynamic  quantities  are  referred 
to  in  this  subroutine  as  "State  Variables,"  in  conformity  with  the  terminology  of  modern 
control  theory.  The  listing  for  this  subroutine  shows  how  the  mathematical  model  for 
the  system  is  written  in  the  form  of  FORTRAN  statements. 

It  is  instructive  to  examine  how  information  about  the  behavior  of  different  components 
in  the  test  circuit  is  incorporated  in  the  mathematical  model.  This  will  be  done  briefly 
below.  Reference  will  be  made  to  the  line  numbers  of  the  FORTRAN  program,  which 
are  printed  at  the  right  edge  of  the  listing. 

Lines  2890  through  2920  of  the  listing  in  Appendix  C contain  the  model  for  the  system 
relief  valve.  If  the  system  relief  valve  cracking  pressure,  there  is  no  flow  through  the 
relief  valve — i .e.  , all  the  pump  flow  is  directed  into  the  cylinder.  However,  if  the 
relief  valve  cracking  pressure  (RELFCR)  is  exceeded  by  the  system  pressure  (Y  (1,1)), 
part  of  the  pump  flow  goes  over  the  relief  valve  (RVFLO).  The  quantity  of  this  flow  is 
determined  by  the  relief  valve  gradient  (XM).  An  ideal  valve  would  have  a infinite  grad- 
ient (i.e.,  the  pressure  is  always  equal  to  tne  cracking  pressure,  no  matter  what  the  flow 
through  the  valve).  In  actual  practice  relief  valves  have  a finite  gradient  which  can  be 
ascertained  during  tests  of  regulation  characteristics,  Refs.  [5,6]. 

FORTRAN  statements  in  lines  2960  through  3050  describe  the  dynamic  behavior 
of  the  test  cylinder  as  well  as  the  frame.  This  is  done  by  furnishing  expressions  for 
the  rate  of  change  of  the  state  variables  (system  pressure,  rod  displacement  and  rod 
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velocity,  respectively).  Thus,  in  line  2960,  the  rate  of  change  of  pressure  is  expressed 
in  terms  of  the  inflow,  rod  velocity,  and  cylinder  leakage.  Note  that  the  fluid  bulk  modu- 
lus (BETA)  as  well  as  the  upstream  volume  of  the  test  setup  (VA)  enter  in  the  expression 
It  is  seen  that  increasing  the  bulk  modulus  and  pump  inflow  tends  to  increase  the  pressure 
rise  rate  while  increasing  the  cylinder  volume  tends  to  reduce  it.  It  may  also  be  noted 
that  when  the  inflow  (QA)  is  exactly  equal  to  the  leakage  flow,  and  the  rod  velocity  is 
zero,  the  pressure  rise  rate  falls  to  zero. 

Statements  in  lines  2970  through  3000  describe  the  dynamic  flexure  of  the  test 
frame.  The  pressure  in  the  test  cylinder  causes  the  test  frame,  pins,  etc.,  to  deflect — 
the  exact  amount  of  deflection  depending  upon  the  effective  stiffness  of  the  test  setup. 

For  a given  pressure  a stiff  system  would  deflect  less  than  one  which  is  less  stiff.  In 
order  to  account  for  the  fact  that  there  will  generally  be  clearances  in  the  cylinder  rod- 
end  and  eye,  the  mathematical  model  incorporates  expressions  for  taking  such  clearances 
into  account.  Basically,  the  mathematical  model  assumes  that  the  frame  stiffness  is 
zero  until  the  slack  due  to  pin  clearance  is  taken  up  (i  .e.  , the  cylinder  rod  moves  freely 
until  the  clearance  is  taken  up).  Beyond  that  point  the  frame  deflects,  building  up  a 
resisting  force  which  ultimately  balances  the  cylinder  pressure. 

DIFFUN  is  called  by  DIFSUB  a number  of  times  during  the  numerical  integration 
process.  After  each  call  DIFFUN  returns  the  values  of  the  derivative  vector  DY  to 
DIFSUB.  Parameters  needed  for  calculating  these  derivatives  are  communicated  via 
the  labeled  COMMON  block  BLK1  from  the  MAIN  program. 
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D I F SUB 


This  subroutine  performs  the  numerical  integration  to  develop  the  time-histories  of 
the  state  variables  ( i . e . , systerr  pressure,  rod  di splacerrent , and  rod  velocity)  using 
the  mathematical  model  furnished  in  the  subroutine  DIFPUN.  It  is  called  b>  the  main 
program  repeatedly  until  the  desired  time-history  is  developed  or  integration  fails  due  to 
abnormal  conditions.  The  attached  listing  explains  the  significance  of  all  quantities 
of  interest.  The  subroutine  has  three  options  for  performing  numerical  integration, 
of  which  only  the  third  (M  = 2)  is  actually  used  in  the  locked-rod  simulation  program. 
Complete  details  of  the  mathematical  basis  for  DIFSUB  is  given  by  Gear,  Ref.  [4], 
and  w i 1 1 not  be  discussed  here . In  addition  to  DIFFUN,  DIFSUB  requires  two  sub- 
routines called  MAT  I NV  and  PEDERV.  The  first  is  needed  for  inverting  a matrix 
generated  by  DIFSUB.  The  second  is  really  needed  only  for  the  option  in  numerical 
integration  method  corresponding  to  MF  - 1 and  is  consequently  a dummy  subroutine 

in  the  present  package. 

MAT  I NV 

This  subroutine  is  needed  by  DIFSUB  to  invert  a matrix.  The  actual  inversion 
is  carried  out  in  a subroutine  called  MINV,  contained  in  the  IBM  Scientific  Subroutine 
Package  (SSP)  and  written  in  FORTRAN.  The  SSP  subroutine  called  ARRAY  is  used 
to  arrange  both  the  matrix  to  be  inverted  and  the  matrix  after  inversion  in  forms  suitable 
to1-  further  processing.  Both  MINV  and  ARRAY  can  be  substituted  by  functionally 
equivalent  subroutines. 
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PROGRAM  USAGE 


The  program  is  setup  to  read  pertinent  data  off  the  card  reader,  designated  as 
until  5 in  the  READ  statements.  Output  is  directed  to  unit  6 in  the  WRITE  state- 
ments. No  scratch  files,  tapes  or  auxiliary  storage  facilities  are  required.  It  is 
assumed  that  the  SSP  subroutines  MINV  and  ARRAY  are  available  in  the  system 
library.  If  this  is  not  so,  source  or  object  decks  of  these  subprograms  can  be  concate- 
nated to  the  program  or  substituted  by  functionally  equivalent  subprograms.  Source 
listings  of  these  subprograms  can  be  found  in  Ref.  [7].  The  remainder  of  this  section 
is  devoted  to  explaining  the  data  input  for  the  program  and  interpreting  the  output. 

Table  I presents  the  description,  units,  and  computer  names  anf  format  specifications 
for  all  variables  and  parameters  needed  for  simulating  the  locked-rod  cylinder  test 
setup.  In  accordance  with  the  standard  practice  in  using  FORTRAN,  all  real  quanti- 
ties (i.e.,  those  whose  computer  names  start  with  letters  A through  H,  and  O through 
Z)  should  conform  to  the  F,  E,  or  D formats,  and  all  integer  quantities  (i.e.,  those 
whose  computer  names  start  with  letters  I through  N)  should  conform  to  the  I format 
and  be  right  justified  in  the  field.  The  computer  program,  written  for  use  on  the 
IBM  360/370  or  equivalent  systems  uses  double  precision  for  all  real  quantities  except 
those  starting  with  "P". 

Figure  3-1  includes  excerpts  from  the  computer  print-out  for  an  example  simulation. 
The  first  part  gives  the  title  of  the  test  system  and  its  parameters.  The  next  part  indi- 
cates the  numerical  integration  method  used  and  the  integration  parameters.  The  third 
part  of  the  output  includes  the  initial  values  for  the  state  variables  (i.e.,  system  pressure, 
piston  rod  displacement  and  piston  rod  velocity).  All  of  the  above  mentioned  parts  of  the 


TABLE  3-1.  INPUTS  TO  LOKROD 


value:  same  as  TPRESS) 
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LOCKED  ROD  CYLINDER  TEST  SIMULATION  770113 


SYSTEM  PARAMETERS 


INFLOW 

FRAME  STIFFNESS 
DRAG  COEFF  IC  IENT 
CY  L I NOE  R AREA 
LEAKAGE  COEFFICIENT, 

UPSTREAM  VOLUME 

FLUID  BULK  MODULUS 

TOTAL  SLACK  I N FRAME 

RELIEF  VALVE  CRACKING  PRESSURE 

RELIEF  VALVE  CRACKING  GRADIENT 


*.0*250*01 
1.570*0*06 
l.  33**D*  OZ 
1.96300*01 
l.*l 70D-03 
3.  63000*  02 
2 • 0000D*0  5 
8.00000-03 
1.90000*03 
*.0*250-01 


CU.  IN.  /SEC 
LBS -F/ IN, 
LBS.F*S£C/!N 
SO.  IN.  _ 

IN  ./  S EC  *P  S I 
CU.  IN. 

PS! 

IN..F*SEC**2/ IN,/ 
PSI 

CU.  IN. /SEC *P SI 


INTEGRATION  METHOD:  GEAR'S  IWlTH  NUMERICAL  DIFFERENTIATION! 

STIFF  SYSTEM  INTEGRATION 


INTEGRATION  PARAMETERS 


ORDER  OF  SYSTEM 

MAXIMUM  URCER  OF  INTEGRATION 

PRINT-OUT  INTERVALS 

I N IT  IAL  STEP  SIZE 

MINIMUM  STEP  SIZE 

MAXIMUM  STEP  SIZE 

ALLOWABLE  ERROR 

MAX.  NUMBER  OF  STEPS 

INITIAL  T |ME 

TIME  AT  PRESSURE 

FINAL  TIME 


3 

3 

5 

l .00000-06 
1 .00000-07 
2 .50000-01 
l. 00000-01 
10250000 
0.0 

1.00000*01 
1 .00000*01 


INITIAL  state  vector 


PRESSURE 

ROD  DISPLACEMENT 
ROD  VElOC  I IV 


l.  -,  7000*01  PSI 
C.r  IN. 

0.0  IN/SEC 


Fig.  3-1.  Computer  Print-Out  for  Locked-Rod 
Test  Simulation 
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Fig.  3-1.  Continued 


RAJECT0R1ES 

NUMBER  OF  POINTS: 

193 

T IME 

PRESSURE 

ROD  DISPLACEMENT 

ROD  VELOC 1 1 

l SEC) 

( PSI  1 

( 1NS1 

I INS/SEC i 

1. 00000-06 

1 .4 1220*01 

2. 88940- 1 0 

2. 88950-04 

5.  23040-04 

2 .5498U+01 

6.2  7890-05 

2.03200-01 

1.23360-03 

3. 7o4l0+01 

3. 69970-04 

5.95590-01 

2.16680-03 

4. 78270*01 

1.29140-03 

1.26780+00 

3.44?  30- 03 

5. 12420*01 

3.51 110-03 

2 . 1694D+00 

4.32280-03 

4.90770*01 

5. 46160-03 

2.05130*00 

5.029/0-03 

5.40250*01 

6.40960-03 

4. 32630-01 

5.55560-03 

6.  72 640  * 0 1 

6.22390-03 

-1 .06130*0 0 

6.13350-03 

8. 84460+01 

5.  3 91  aO-  0 3 

-1.60650*00 

6. 83190-03 

1.12890+02 

4.46960-03 

-8. 74280-01 

7. 59850-03 

1. 29990+02 

4. 33420-03 

5. 59190-01 

8.62410-03 

1 .35600*02 

5. 73290-03 

1. 87000*00 

9. 37530-03 

1.38810*02 

6.83770-0  3 

6.73410-01 

9.84700-03 

1.48710*02 

6. 79840-03 

-7.27050-01 

1.03110-02 

1 .63710*02 

6. 26420-03 

-1.42230+00 

1.06590-02 

1 . 75850  + 02 

5.7  743  D— 0 3 

-1.31800*00 

1.10860-02 

1.8d61D*02 

5.  362  00-03 

-7.18770-01 

1.14890-02 

1 .97420*02 

5.26770-03 

1.  15950-01 

1. 17380-02 

2. 00990*02 

5.3795 D-03 

6.54350-01 

1. 18980-02 

2.02600+02 

5.  51390-03 

9.62950-01 

1. 23820-02 

2 .04910+02 

6.15710-03 

1. 53920*00 

1.29070-02 

2. 06690*02 

6.92080-03 

1.18930*00 

1.31960-02 

2 .09660*02 

7.  15580-03 

4.35230-01 

1.34230-02 

2.13960*02 

7. 15740-03 

-2.56040-01 

1.37100-02 

2.21530+  02 

6. 96050-03 

-9.0260D-01 

1.41370-02 

2.34700*02 

6.48310-03 

-1.20240*00 

1.48210-02 

2.53260+02 

5.93890-03 

-2.02670-01 

1.54720-02 

2.61820*02 

6.25100-03 

1.08270*00 

1.60200-02 

2.64370+02 

6.93800-03 

1. 20300*00 

1. 6580D— 02 

2.  7053D+  02 

7.31060-03 

3.51820-03 

1.69220-02 

2. 7833D+02 

7. 15970-03 

-7.  18970-01 

l. 72560-02 

2 . d7580+02 

6.85870-03 

-9. 84180-01 

1. 76060-02 

2. 97190*02 

6. 54500-03 

-7.53680-01 

1.81870-02 

3.09070+02 

6. 39300-03 

2.97100-01 

1. 85190-02 

3.12770*02 

6.58890-03 

8.52900-01 

1.90350-02 

3.  16040*  02 

7.  11920-03 

1.05240+00 

3.55170-02 

5.64080+02 

7. 71 140-03 

3. 2710D-01 

3. 69340-02 

5. 79930*02 

8.00900-03 

-2.7640D-02 

3.31 190-02 

5.97410*02 

7. 84040-03 

-9.59810-02 

4.  10910-02 

6.33690*02 

7.98570-03 

5.60590-02 

1.  19540-01 

1.21730+03 

6.25880-03 

1.32190-01 

1.19700-01 

1.21o30*03 

8.40290-03 

1 . 44900*00 

1.20020-01 

l .2097D+03 

9.12620-0  3 

2.65160*00 

1.  20370-01 

1. 20450*03 

9.  72  700-0  3 

6.32630-01 

1.20530-01 

1.20560*03 

9.68100-03 

-9. 66680-01 

1.20710-01 

1 .209  30  + 03 

9.39710-03 

-2.04770*00 

1.20920-01 

1.21543*03 

8.91100-03 

-2.241 70*00 

1.21090-01 

1.21930*03 

8.60110-03 

-1.53090*00 

1.21420-01 

1. 22180*03 

8.48670-03 

9.16140-01 

1 .21850-01 

1.2 1460*03 

9.  28880-03 

2.12540*00 

1.22320-01 

1.21290*03 

9.61080-03 

- 1. 09630*00 

1.22440-01 

1.22020*03 

9.0442D— 03 

-2.17750*00 

1.23030-01 

1.22770+03 

8.47950-03 

-2.  33150-01 

1.23440-01 

1 .22420*03 

8.92660-03 

2.10510*00 

1.23  700-01 

l. 21950*03 

9.  45360-03 

1.70090+00 

1.23880-01 

1.21790*03 

9.  65530-03 

4.02040-01 

1 . 24  C50-01 

l .21920*03 

9. 58950-03 

-8.62230-01 

1.24230-01 

1.22280*03 

9.  3 2 3 1 D-  0 3 

-1.74750*00 

1.24470-01 

1.22830*03 

8. 8858D-03 

-1. 74750*00 

1.  24730-01 

1.23240*03 

8.58910-03 

-5.02850-01 

1.24840-01 

1.23290*03 

8.  57820-03 

2.16600-01 

1.25010-01 

1 .23210*03 

8.71230-03 

1. 13750+00 

1.25250-01 

1.  22890*03 

9.07930-03 

1.77390*00 

1.25580-01 

1.22510+  03 

9. 54090-03 

6.888  70-01 

1.25850-01 

l .22660*03 

9.4864D— 0 3 

-1.05880*00 

1.26170-01 

1. 23300*03 

8.99190-03 

-1.63420*00 

1.26470-01 

1 .23760*03 

8.66240-03 

-3. 86840-01 

1.26830-01 

1 .23650*03 

8.87280-03 

1.40720 *00 

1.27160-01 

1.23210*03 

9.  3806D-0  3 

1.29770+00 

1.27590-01 

1.23260*03 

9.  46  890-  03 

-9.  5*540-01 
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print-out  present  the  input  information  so  that  the  program  user  may  easily  verify  that 
the  correct  information  is  being  furnished  to  the  computer. 

The  last  part  of  the  print-out  includes  the  time-histories  for  the  state  variables 

printed  out  in  increments  of  the  print-out  interval  MULTSP.  The  total  number  of 

points  in  the  trajectory  is  also  indicated.  Digital  simulations  is  essentially  a discrete 

process  and  establishes  the  values  of  the  state  variables  at  discrete  points  in  time. 

Values  at  other  points  in  time  can  be  easily  obtained  by  interpolation.  The  common 

practice  is  to  plot  the  data  points  graphically,  either  manually  or  via  a digital  computer. 

From  the  print-out  it  can  be  seen  that  the  time  increments  between  the  data  points 

-3  -3 

increases  from  .523  x 10  seconds  at  the  beginning  to  2.9  x 10  seconds  at  0.3811 
seconds.  This  change  in  step  size  is  automatically  brought  about  by  the  integration 
algorithm  and  ensures  that  simulation  is  done  efficiently  and  quickly. 
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CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 


Hydraulic  cylinders  are  important  components  of  mobile  hydraulic  systems  and 
need  to  be  carefully  evaluated  to  ensure  reliable  operation  in  the  field.  Apart  from 
being  hydraulic  components,  cylinders  are  also  structural  elements  and  consequently 
need  to  be  appraised  for  structural  integrity.  Structural  failure  of  hydraulic  cylinders 
may  occur  due  to  static  overloading  (usually  buckling  failure)  or  due  to  cyclic  loading 
(fatigue).  The  joint  flexure  theory,  developed  by  the  Fluid  Power  Research  Center, 
for  static  analysis  of  hydraulic  cylinders  has  been  computerized  and  used  for  analyzing 
a variety  of  cylinders,  not  only  at  the  Fluid  Power  Research  Center  but  also  industrial 
companies  sponsoring  the  Technology  Development  Program  on  Cylinder  Integrity 
Assessment.  Dynamic  analysis  of  locked-rod  (impulse)  testing  has  been  presented  and 
conputerized  for  digital  simulation.  Test  results  confirm  that  it  is  sufficient  to  measure 
pressure  and  deduce  the  cylinder  and  rod  stresses  by  applying  standard  formulas, 
rather  than  to  strain  gage  the  cylinders. 

Digital  simulation  of  the  test  setup  for  locked-rod  cylinder  tests  indicate  that  within 
certain  limits  the  pin  clearance  does  not  significantly  affect  the  pressure  rise  rate. 
However,  the  pressure  setting  of  the  system  relief  valve  significantly  influences  the 
pressure  rise  rate  and  consequently  caution  should  be  exercised  in  comparing  test 
results  on  cylinders  of  different  pressure  ratings. 
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A survey  of  cylinder  test  reports  listed  in  Ref.  [2],  Section  I , Appendix  B,  indi- 


cate that  pin  and  rod-eye  failures  are  the  more  common  modes  of  failure  in  impulse 
testing.  Since  the  effect  of  pin-eye  clearance  is  to  cause  cyclic  impact  loading  , 
these  elements  must  be  designed  for  dynamic  loading  rather  than  just  static  loads. 
Even  design  for  fatigue  stregth  using  classical  S/N  diagrams  may  be  inadequate  since 
such  curves  are  usually  generated  from  sinusoidal  loading  rather  than  impact  loading. 
Fatigue  tests  conducted  to  failure  on  rod-eye  arrangements  subjected  to  impact  loads 
are  recommended  as  a method  of  acquiring  a reliable  data  base.  Investigation  of 
stress  concentration  effects  on  both  static  as  well  as  dynamic  loading  is  also  recom- 
mended for  future  investigation. 
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APPENDIX  A 


DRAFT  TEST  PROCEDURE  FOR  IMPULSE  AND  STROKING  TESTS 
ON  HYDRAULIC  CYLINDERS 


IMPULSE  (LOCKED-ROD)  TEST 

1 . Scope:  This  test  procedure  is  intended  to  be  applicable  to  all  hydraulic  cylin- 

ders which  are  not  integral  parts  of  other  components. 

2.  Purpose:  The  objective  of  the  test  described  herein  is  to  ascertain  the  capability 

of  a hydraulic  cylinder  to  repeated  pressure  cycles. 

3.  Materials  & Test  Apparatus 

3.1  Test  fluid  shall  be  that  designated  for  use  in  the  hydraulic  system 
in  which  the  cylinder  is  used. 

3.2  The  test  frame  in  which  the  cylinder  is  installed  for  this  test 
should  be  sufficiently  stiff  so  that  the  relative  motion  between  the 
piston  and  the  cylinder  does  not  exceed  0.1%  of  the  stroke  of  the 
test  cyl  inder . 

3.3  The  test  cyl  inder  sha  1 1 be  mounted  usi  ng  the  mounts  (i . e . , pins  , 
trunnions,  etc.)  used  in  installing  the  cylinder  in  the  intended 
application.  Pin  material,  clearances  and  lubrications  arrange- 
ments shall  conform  to  that  used  in  the  intended  application. 

3.4  The  test  frame  shall  have  provisions  to  adjust  and  measure 
eccentricity  of  the  test  cylinder  to  the  center  line  of  loading. 

3.5  The  test  frame  shall  not  impose  side  loads  or  torque  on  the 
test  cyl  inder . 

3.6  The  hydraulic  circuit  shall  be  as  shown  in  Fig.  A-1  or  equivalent. 
The  system  should  bo  capable  of  maintaining  the  prescribed  pres- 
sure waveform  for  the  entire  test  duration. 

4.  Measurement  Accuracy 

All  physical  quantities  shall  be  measured  to  the  accuracies  nivon  below: 

4.1  Pressure:  + 2%.  The  pressure  transducer  shall  have  been 

a flat  frequency  response  from  0-6  kHz. 


PH hf. tail  ;;s  FAGS  HOT  FUJGD 
BLANK 


Fig.  A-l  . Hydraulic  Circuit  Schematic  for  Blocked-Rod  Test 

4.2  Flow:  ± 2%. 

4.3  Temperature:  ± 5°  F ( ±2.8°C). 

4.4  Eccentricity:  ± 0.05%  of  stroke  of  test  cylinder. 

4.5  Time:  ±0.05  seconds 
5.  Test  Conditions  & Definitions 

5.1  Test  system  fluid  shall  be  at  a temperature  of  1 80°F  ± 5°F 
(82°C  ± 2 . 8°C) . 

5.2  Test  fluid  shall  be  maintained  clean  and  shall  have  a particulate 
matter  content  of  not  more  than  1500  particles  per  milliliter  of  size 
greater  than  10  micrometers. 

5.3  Test  waveform  shall  conform  to  particulars  given  in  Fig.  A-2. 

5.4  Packing  drag  is  defined  as  the  pressure  needed  to  move  the  piston 
with  all  ex'ernal  forces  (mechanical  or  hydraulic)  removed  or  in 
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AVERAGE  PRESSURE  RISE  RATE 


Fig.  A-2.  Pressure  Waveform 
for  Cylinder  Impulse  Test 


T Time  Period  tot  One  Cycle  60  NMinc  i 
N Cycling  Ffwjuency  (cycles  per  minute! 

T Time  On  Ptetsu>«  isec  I 
P,  Mammum  feit  Pressure  psi  MPa 
P„  Minimum  Test  Pressure  ips.  MPai 
AP,  Pressure  Overshoot  (pti  MPa 
^ Average  P'essuie  Rise  Rate  between  15^ 
and  90%  at  P,  .psi  Sec  or  MPa  sec > 

N Shall  Be  Between  29  aryl  31  Cycles  Per 

0 046 T T 0 55 T 20  000  P 80  000  Ipsi  seel 
AP,  0 1P, 

P,  0 IP, 


balance . 

5.5  Piston  dirft  is  the  measure  of  internal  leakage  of  the  cylinder 
at  the  specified  pressure. 

5.6  Failure  is  any  condition  which  precludes  the  use  of  the  cylinder 
for  its  intended  application.  It  is  evidenced  by  excessive  leak- 
age— internal  or  external,  structural  deformation  or  breakage 
of  any  part  of  the  cylinder  and  its  mountings. 

Test  Procedure 

6.1  Install  the  test  cylinder  in  the  test  frame. 

6.2  Adjust  pump  flow  and  relief  valve  setting  to  prescribed  values. 

6.3  Adjust  the  timer  for  the  solenoid  directional  control  valve  so  as 
to  obtain  the  prescribed  waveform. 
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STROKING  TEST 


1 . Scope:  This  test  procedure  is  intended  to  be  applicable  to  all  hydraulic 

cylinders  which  are  not  integral  parts  of  other  components. 

2.  Purpose:  The  objective  of  the  test  described  herein  is  to  ascertain  the  capability 

of  a hydraulic  cylinder  to  piston  rod  stroking  under  realistic  loading 
conditions . 

3.  Materials  & Test  Apparatus 

3.1  Test  fluid  shall  be  that  designated  for  use  in  the  hydraulic 
system  in  which  the  cylinder  is  used. 

3.2  The  test  cylinder  shall  be  installed  in  a test  frame  and  test  circuit 
as  shown  in  Fig.  A-3.  The  slave  cylinder  may  be  identical  to 
the  test  cylinder.  If  it  is  different  it  should  have  a stroke  and 
piston  rod  length,  at  least  equal  to  the  test  cylinder's  and  shall 
not  exceed  110%  of  the  above  mentioned  values. 

3.3  The  connecting  piece  between  the  piston  rods  for  the  test  and 
slave  cylinders  shall  not  exceed  10%  of  the  test  cylinder  stroke. 

3.4  The  test  frame  shall  have  privisions  to  adjust  and  measure  eccen- 
tricity of  the  test  cylinder  to  the  center  line  of  loading. 

3.5  The  test  and  slave  cylinders  shall  be  mounted  so  as  not  to  impose 
any  side  loads  on  the  test  or  slave  cylinders. 

4.  Measurement  Accuracy 

All  physical  quantities  shall  be  measured  to  the  accuracies  given  below. 

4.1  Pressure:  ± 2%. 

The  pressure  transducer  shall  have  a flat  frequency  response 
from  0-5  kHz . 

4.2  Flow:  ± 2%. 

4.3  Temperature:  ± 5°F  ( ' ± 2.8°C) 

4.4  Eccentricity:  ± 0.05%  of  stroke  of  test  cylinder. 

4.5  Time:  ± 0.05  seconds. 

5.  Test  Conditions  & Definitions 

5.1  Test  system  fluid  shall  be  at  a temperature  of  1 80°F  ± 5°F 

(82°C  ± 2.8°C)  during  the  first  10%  of  total  cycles  and  1 50°F 
5°F  (63°±  2.8°C)  during  the  balance  of  the  test  run. 


5.2  Test  system  fluid  shall  be  maintained  clean  and  shall  have  a 
particulate  matter  content  of  not  more  than  1500  particles  per 
milliliter  of  size  greater  than  10  micrometers. 

5.3  The  average  pressure  rise  rate,  as  measured  between  15%  and 
90%  of  the  maximum  operating  pressure  of  the  cycle,  shall  not 
be  less  than  20,000  psi/sec. 

5.4  Packing  drag  is  defined  as  the  pressure  needed  to  move  the  piston 
with  all  external  forces  (mechanical  or  hydraulic)  removed  or  in 
balance. 

5.5  Piston  drift  is  the  rteasure  of  internal  leakage  of  the  cylinder  at 
the  specified  pressure. 

5.6  Failure  is  any  condition  which  precludes  the  use  of  the  cylinder 
for  its  intended  application.  It  is  evidenced  by  excessive  leak- 
age— internal  or  external,  structural  deformation  or  breakage  of 
any  part  of  the  cylinder  and  its  mountings. 

Test  Procedure 

6.1  Install  the  test  cylinder  in  the  test  frame.  Connect  the  test 
cylinder  to  the  slave  cylinder  and  check  and  adjust  alignment. 

6.2  Adjust  pump  flows  and  relief  valve  settings  so  as  to  attain  the 
prescribed  maximum  operating  pressure  and  the  pressure  rise 
rate.  The  latter  can  generally  be  adjusted  by  changing  the  pump 
flow  rate,  the  length  and  volume  of  connecting  lines  and  the  speed 
of  the  directional  control  valve. 

6.3  Perform  piston  drift  and  packing  drag  tests. 

6.4  Stroke  the  cylinder,  in  accordance  with  the  prescribed  schedule, 
by  switching  the  directional  control  valve.  The  pressure  rise 
rate  and  the  maximum  operating  pressure  shall  be  recorded  at 
least  twice  during  the  test;  once  at  t foe  beginning  of  the  test  and 
once  at  the  end.  Additional  records  of  pressure  waveform  should 
be  taken  when  any  irregular  operation  or  malfunction  is  observed. 
Stroking  shall  be  done  for  the  prescribed  number  of  cycles  or  until 
failure.  No  repairs  or  modifications  shall  be  carried  out  on  the 
test  cylinder  at  any  time  in  the  course  of  the  cycling. 

6.5  Ferform  piston  drift  and  packing  drug  tests  at  the  termination  of 
the  eye  I inf). 
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TEST  FRAME 


CYLINDER  UNDER  TEST 


CYLINDER  - SLAVE  & OR  TEST 


Fig.  A-3.  Test  Frame  and  Test  Circuit  for  Stroking  Test 
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APPENDIX  B 


LISTING  FOR  LOKROD 


RELEASE  2.0 


DATE  = 7701* 


15/09/38 


LOCKED  R00  CYLINDER  TEST  SIMULATION  770113 


IMPLICIT  REAL*8( A-H, Q-Z) 

C OMMON/ INIT/YYINT 110) 

C0MM0N/3R0ER /N 

COMMON/ BLK1/  OP , ST  IF F , DR  AG , AA , CL . V A , BE  TA , X I , CLRNCE , RE LFCR . XM , 
ITPRESS 
C 

DIMENSION  VI 8, 10) 

DIMENSION  PTITLE118) 

DIMENSION  TR AJ ( 10, 1200) 

DIMENSION  SAYEI10.10),  YMAXI10I.  ERRORUO),  PW<10,10) 

OIMENSION  DERI VIIO) 

OAT  A YMAX/ 10*1 .000/ . ERROR/ 1 0* 1.  00- 1/ . J /- 1 / 

RE  AD  1 5 • l ) PT  ITLE 

1 FORMA  Til 8AAI 
WRITE(6,21  PT1TLE 

2 FORHATIlHl  ,/////, 5X  , 1 8 Aa  , // // ) 

READ! 5. 10)  OP. ST  IFF, DRAG  . A A ,C L , V A , BE T A , X I , CL RNCE , RE LF CR,  X M 
WRITE  16, 3)  OP, ST  IFF, DR AG, A A, CL, V A 
WRITE  16, A)  BETA, XI  ,CL RNC E , K E LFC R . X M 

3 FORMAT  I 1H  , T 5, 1 7H S YS TE M P AR AME TE R S , / / 1 H , 

1 T10.6HINFL0W.TA0, 1PE13.A, 1 lH  C J . I N . / SEC , / 1H  , 

2 T 1 0, 1 5HFR AME  ST  I F F NE S S , T AO  , 1 PE  1 3. A , 1 0 H LBS.F/lN./lH  , 

3 T 10,  16HDR  AG  COE  FF  IC  I E N T , T A 0,  1PE  l 3 . A , 1 3H  L 8 S.  F * S EC / I N , / 1 H , 

A T 1 0 , 1 3HCY LINDER  ARE  A, T AO , 1 PE l 3. A, 8H  SO.  IN../1H  , 

5 T10, 2 OH LEAKAGE  COEF F I C I E NT , TAO, I PE 1 3 . A , 1 2 H I N ./ S EC *PS I , / 1 H , 

6 T10, 15HUPSTREAM  VOL  UME  , T AO  , IPE  1 3.  A , 8H  CU.  IN.) 

A FORMATUH  , T 10,1  8HFLU  I 0 BULK  MODUL  US  , T AO  , IP  E 1 3 . A , AH  PS1./1H  , 

1 T 10,  18HPI  STON  ROD  I NE  R T I A , TA  0 , 1 PE  l 3 . A , 1 8 H 1 BS -F  *S  E C* * 2 / I N,  / IH  , 

2 T10.20HT0TAL  SLACK  IN  F R AME , T AO , 1 PE  1 3 . A , AH  IN../IH  , 

3 T10, 30HRELI EF  VALVE  CRACKING  PRES S U R E . T AO , 1 P E 1 3 . A, AH  PSI./1H  , 

A T 10,  30HREL  IEF  VALVE  CRACMNG  GR  AO  I E NT  , T A 0 , 1 PE  1 3 . A , 

5 15H  CU.  IN./SEC*PS I,//// ) 

INTEGRATION  PARAMETERS: 

N ORDER  OF  SYSTEM 

MF  METHOO  NUMBER 

MAXDER  MAXIMUM  ORDER  OF  INTEGRATION 

MULTSP  PRINT-OUT  INTERVALS 

H INITIAL  STEP  SUE 

HMIN  MINIMUM  STEP  SIZE 

HMAX  MAXIMUM  STEP  SIZE 

EPS  ALLOWABLE  ERROR 

NST EPS  MAXIMUM  NUMBER  OF  STEPS 

TINTL  INITIAL  TIME 

TPRESS  TIME  AT  PRESSURE 

TFINAL  FINAL  TIME 

READ! 5,9)  N , MF , MA XDE R , MUL T S P 

9 F ORM  AT (815) 

REAO (5,10)  H , HMI N, HMAX, EPS, TINTL , T PR ES S , T F INAL 
READ! 5,  10)  (YYINT(J)  ,J=1  , N) 

10  FORMAT (7F10.3) 

TRAJ ll,l)>  TINTL 
NP  1=  N* 1 


00001000 

ooooloio 

00001020 

**00001030 

00  00 1 0 AO 
00001050 
00001060 
00001070 
00001080 
00001090 
00001 100 
00001 1 10 
00001120 
0000 11 30 

00001  1 AO 
000011 50 
00001 1 60 
000011 70 
00001180 
00001190 
00001200 
000012 10 
00001220 
00001230 
00  00 1 2 AO 
00001250 
00001260 
000012  70 
00001280 
00001290 
00001300 
00001310 
00001 320 
00001330 
000013A0 
00001350 
00001360 
00001370 
00001 3bO 
00001390 

00  00 1 A 00 
000C1A10 

OOOOlA/O 

00001  A 30 
00001 AAO 

00  00 1 A 50 
00001A60 

00001 A70 
0000 1 A 80 

00001  A 90 
00061 500 
00001510 
00001520 
000C1 530 
0000 1 5 AO 
00001550 
0000  l 5 C>G 
000015  70 
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RELEASE  2.0 


15/09/38 


INITIALIZE  TRAJECTORY  VECTORS 
00  50  JJ=2.  NP1 
TRAJIJJ,1I=  YYINTCJJ-LI 
T=TI NTL 

IFIHMAX.GT.TF INAL » HMA  X=  TF I NAL 

NSTEPS=(TFINAL*HMAX)/H 

J STAR  T«0 

00  30  JJ=1,  N 

Y ( 1 i JJJ  = YY INT ( J J» 

IF(HF.EQ.O)  WRITEI6.15) 

FORM  AT ( 1 H ,T5,27HINTEG RATION  METHOD:  ADAMS, ////I 

IFIMF.EQ.l)  WRITE(6«16) 

FORMAT ( 1H  ,T5,46HI NTEGRATION  METHOD:  GEAR'S  (WITH  ANALYTICAL  , 

112HDER1VAT IVES 1,/lH  , T 2 7 , 24HST l F F SYSTEM  INTEGRATION,////! 
IFIMF.EQ.2)  WRITE(6«17 I 

FORMAT ( 1H  .T5.45HINTEGRAT10N  METHOO:  GEAR'S  (WITH  NUMERICAL  , 

116HDIFFERENT I AT  ION ) . / 1 H ,T  2 7, 24HST IFF  SYSTEM  INTEGRATION,////) 

WRI TE ( 6,20)  N.MAXDER.MULTSP 

FORMAT  1 1H  .T5.22HINTEGRATION  PARAME TER S, // 1H  , 

1 T 1 0 , 15H0RDE  R OF  S YST EM, T43 , 15, / 1H  , 

2 T10.28HMAX! MUM  ORDER  OF  I N TEGR A T I ON , T 43 , 1 5 , / 1 H , 

3 T 10, 19HPRIN T-OUT  I N T ER V AL S , T43 , 1 5 I 
WRI TEI6,21)  H , HM IN.HMAX, EPS 

F0RMATI1H  ,T10,17H INITIAL  STEP  S I Z E ,T3 0 , i PEI 8 . 4 , / 1H  , 

1 T 10 , 17 HM I N I MUM  STEP  S I Z E , T 30, 1PE 1 8. 4 , / 1H  , 

2 T10,l 7HMAXI MUM  STEP  S I Z E ,T 30 , 1 PEI 8 .4, / 1 H , 

3 T10, 15HALL0WABLE  ERROR , T3 0 . IPE 18 . 4) 

WRIT  E (6,22  ) NSTEPS.T I NTL , T P RES S , T F l NAL 

FORMA T ( 1 H , T 10  ,2  OHMA  X . NUMBER  OF  S T E PS , T 3 3 , 1 1 5 , / 1 H , 

1 T10,  1 2H  INITIAL  TI ME , T 30 , IPE 18.  4 ,/ 1H  , 

2 T10.16HTIME  AT  PR ES SU RE , T 3 0, 1 P E 18 .4 , / 1H  , 

2 T10.10HFINAL  TI ME  ,T 30 , 1 PE l 8. 4 , / // / ) 

WRITEI6.23)  (YVINTIJJI,  JJ»1,  N) 

FORMATUH  ,T5,20HINIT  IAL  STATE  VECTOR,//, 1H  , 

1 T 10, 8HPRE SSURE , T2 7 , 1 PE 1 2. 4 ,1 X , 3 HP S I ,/lH  , 

2 T 10 , 16HR00  01 SPLACEMENT ,T2 7, IPE 12.4, IX, 3HIN. ,/lH  , 

3 T 1 0 , 1 2HR00  VELOCITY, T27.1PE12. 4, IX, 6H IN/S  EC, /// / ) 

MULT S P = NUMBER  OF  STEPS  BETWEEN  PRINTOUT 
KOUNT  = MULT  SP 
00  100  NNN= 1 ,N  STEP  S 
NNNP 1 = NNN*l 

IF(T.GT.TFINAL)  GO  TO  9120 
IFINNN.GT.l)  J ST  AR  T= 1 
CONT I NUE 

CALL  01FSUBI N, T, Y , SAVE, H,HM IN, HMAX.EPS, MF , YMAX.E  RROP , KF LAG , J ST  ART  , 
1 MAX  OE  R , PW  J 

IFIKFL AG.NE  . 1 ) GO  TO  9120 

IFIKOUNT.LT. MULTSP)  GO  TO  109 

ICOL  = 1+  NNNP 1 / MULT  S P 

IF ( IC0L.GE.1200)  GO  TO  9120 

TRAJI 1 , ICOL) =T 

00  111  JJ=2,  NP1 

TRAJI JJ,  ICOL )=  Y I l.JJ-l) 

CONTI NUE 
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0064 

00 65 

KOUNT  = KOUNT  -l 

IFUOUNT  .60.0)  KOUNT=  MULTSP 

0D0U216J 
000021 70 

0066 

100 

CONT I NUE 

00002180 

C 

00002 190 

C 

00002200 

c 

00002210 

0067 

9120 

CONTINUE 

00002220 

0066 

WRITE16.9130)  ICOL 

000022  30 

0069 

9130 

FORMAT!///, 1H  ,T5, 12HTRAJECT0RI ES, 5X, 1 7HNUMBER  OF  P 01 N T S: , I 8 • / / / / . 

i 00002240 

11H  ,T9,4HTIME,T25,8HPRESSURE,T41 .16HRO0  D 1 S PL  ACE M EN T , T 62 , 

00002250 

212HR00  VELOCITY, /IH  , T9 , 5H 1 SEC » , T 2 6 , 5H ( P SI)  , T46 , 5H ( I N S ! . 

00002260 

3T63.9H! INS/SECI ,/  ) 

00002270 

0070 

00  9210  JJ=1 , ICOL 

00002280 

0071 

MR  I T E ( 6, 9220  ) ( TR A J( NR , 0 J I , NR-l,  NP  11 

00002290 

0072 

9210 

CONTINUE 

00002300 

0073 

9220 

FORMAT! T5, IP El  3. 4, T21 , 1 PE  13. 4, T42, L PEI  3.4 ,T59,1PE13.4) 

00002310 

0074 

IFIKFLAG.EQ.  1)  GO  TO  9999 

00002320 

0075 

WRITE<6,9140)  T 

00002330 

00  76 

9140 

FORMAT (////  1H  ,T5,44H****  SIMULATION  STOPPED  ABNORMALLY  AT  TIME  =, 

00002340 

1 2X  , 1 PD1 5» 4 ) 

00002350 

0077 

IFIKFLAG.EQ.-l  ) M R I T E ( 6 , 9 1 4 1 J 

00002360 

0078 

9141 

FORMAT!  1H  ,T10,46HSTEP  MAS  TAKEN  MlTH  H = H Ml N , BUT  REQUESTED  ERROR, 

000023  70 

118H  WAS  NOT  ACHIEVED.) 

00002380 

0079 

IF! KFLAG.EQ.-2)  MR  I T E ( fa ,9 1 4 2 ) 

00002390 

0080 

9142 

FORMAT ( IH  , T 10 , 44HTHE  MAXIMUM  ORDER  SPECIFIED  WAS  FOUND  TO  BE  , 

00002400 

110HT00  LARGE.) 

00002410 

0081 

IF(KFLAG.EQ.-3)  WRI TEI 0,9143) 

00002420 

0082 

9143 

FORMAT ( IH  ,T  10.44HCORREC TOR  CONVERGENCE  COULD  NOT  BE  ACHIEVED  , 

00002430 

113HP0R  H.GT.HMIN) 

00002440 

0083 

■ 

IF! KFLAG.E0.-4)  WR I T E ( fa , 9 l 44) 

UJ002450 

0084 

9144 

F0RMAT11H  .T10.43HTHE  REQUESTED  ERROR  IS  SMALLER  THAN  CAN  BE  , 

00002460 

1 2 5HHANDL  ED  FOR  THIS  PROBLEM.) 

000024  70 

0085 

9999 

CONTINUE 

00002480 

0086 

WRITE(6,9980) 

00002490 

0087 

9980 

FORMA  T ! / // / / ) 

00002500 

0088 

STOP 

00002510 

0089 

END 

00002520 

1 
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10 

11 


c 

c 

20 


SUBROUTINE  D IFFUNI T , V ,0Y ) 

IMPLICIT  REAL*8lA-H,Q-Z) 

THIS  SUBROUTINE  EVALUATES  THE  DERIVATIVES  OF  THE  DEPENDENT  VAR  I AB 

LES  STORED  IN  YU, I)  FOR  1=1  TO  N,  AND  STORES  THE  DERIVATIVES  IN 

ARRAY  DY.  N=  NUMBER  OF  FIRST  OROER  DIFFERENTIAL  EQUATIONS 

T = 1 NDEPE NDENT  VARIABLE 

COMMON/INI T/YYINTI 101 

COMMON/ORDER/N 

COMMON/B  LK1 / QP , ST  I FF , DR  AG , AA, CL , V A, BET  A , X l , CLRNC 6,  RELFCR, XM, 
1TPRESS 

DIMENSION  Y ( 8, 3) , DY t 101 

EQUATIONS  IN  ' D I F F UN 1 FOR  SIMULATING  LOCKED  ROD  FATIGUE  TEST 
MOOELEO  AS  A THIRD-ORDER  SYSTEM 

STATE  VARIABLES: 

Y(l,l)  PRESSURE 

Y( 1,21  ROO  DISPLACEMENT 

Y ( l , 3 1 ROD  VELOCITY 

INPUTS: 

QP 

PARAMETERS 
STIFF 
XKF 
DRAG 
AA 
CL 
VA 

BETA 
XI 

CLRNCE 
RELFCR 
XM 


INFLOW 


ULTIMATE  FRAME  STIFFNESS 

FRAME  STIFFNESS 

DRAG  COEFFICIENT 

CYL  INOER  AREA 

LEAKAGE  COEFFICIENT 

UPSTREAM  VOLUME 

FLUID  BULK  MODULUS 

PISTON  ROD  INERTIA 

TOTAL  SLACK  IN  FRAME 

RELIEF  VALVE  CRACKING  PRESSURE 

RELIEF  VALVE  GRADIENT 


GO  JO 
YI  1 ,1)  .LT. RELFCR) 


20 

GO  TO 
RELFCR) 


10 


IFIT.GT.TPRESSI 
IF  ( ' ' . 

RVFLO=  XM*IY( 1, I I - 
QA=  QP  - RVFLO 
GO  TO  11 
QA=  QP 
CONTI NUE 

D YI 1)=  I QA— AA* I YI 1,3) 

DY  12 )=  Y 1 1,  3) 

I F ( YU  ,21-GE. CLRNCE) 

XKF  = (YI 1,2) /CLRNCE)**4  * STIFF 
DY  1 3 ) = I Y(1,1)*AA  —DR AG*Y (1,3) 


RETURN 
CONTINUE 
QA=  0.00*0 
RETURN 
ENO 


♦ CL*Y( 1 ,1)1)  * BETA/VA 


XKF=ST IFF 

* 


- XKF  *Y I 1,2)  I /XI 


00002530 
00002540 
-00002550 
00002560 
00002570 
00002560 
00002590 
00002600 
00002610 
00002620 
00002630 
00002640 
00002650 
00002660 
000026  70 
00002680 
00002690 
00002700 
00002710 
00002720 
00002730 
00002740 
00002750 
00002760 
00002770 
00002780 
00002790 
00002800 
00002810 
30002820 
00002830 
00002840 
00002850 
00002660 
00002870 
3 0 002d  80 
00002890 
00002900 
00002910 
00002920 
00002930 
00002940 
00002950 
00002960 
000029  70 
00002980 
00002990 
00003000 
00003010 
00003020 
000030  30 
00003043 
00003050 
00003060 
000030  70 
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ElEASE  2.0 


?Hl>S0SUBPROGRAMJyNVEftTiMiNJN  BY  N MATRIX  STORED  IN  ARRAY  Pw 
THE  FOLLOWING  SSP  SUBROUTINES  ARE  CALLED  BY  THIS  SUBROUTINE: 


ARRAY  INTERCONVERSION  OF  SINGLE  AND  DOUBLE  DIMENSION  ARRAYS 
M I NV  MATRIX  INVERSION 


THEY  MAY  BE  REPLACED  BY  APPROPRIATE  SUBROUTINES 


DIMENS  ION  PR (M.M) 

DIMENSION  DUMMY  1 10,1 01 ,S11001,LI23 ),MM(20! 
DATA  OUMMY/IOOYO. 0/, S/100*0.0/ 

DO  10  NROW  = 1 , N 
DO  10  NC OL  = 1 , N 

OUMMYlNROw.NCOL 1=  P M ( NRO M , NCOL ) 

10  CONTINUE 
MODE  = 2 

CALL  ARRAY (MODE ,N, N,  10 , 10, S .DUMMY) 

CALL  Ml NVl S . N.DET , L, MM) 

IF ( DE  T.EO.O. 0)  GO  TO  BO 
J=  1 

MODE  =1 

CAL  l ARRAY ( MODE , N , N , 1 U , 1 0, S .OUMMY) 

DO  20  NRnw= i , N 
DO  2 0 NC  OL  =1  , N 

Pw< NROw.NCOL )=  DUMMY! NRUW.NCOL ) 

20  CONTINUE 

RE  TURN 

50  J=-l 

RETURN 

END 


00003080 

00003040 

J0J031O3 

000031 10 

00003120 

00003130 

00003140 

000031 50 

00003160 

000031 70 

00003180 

00003140 

00003200 

000032  10 

0000  3220 

00003230 

00003240 

00003250 

000032 

00003270 

00033283 

000032  40 

00003300 

00003310 

00003320 

00003333 

00003340 

0 0003350 

0000  3360 

00003373 

00003380 
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IMPLICIT  REAL*8< a-h.o-z J 
DIMENSION  PRIM, Ml 
RE  TURN 
END 
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000034 13 
03003423 
00003430 
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00003453 
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SUBROUTINE  0 I FSU  B I N,  T . Y , S AV  E.  H,  HM  I N,  HM  AX  , EPS  , MF,  YMA X , E RROR . KF L AS , 00003460 
IJSTART.MAXDER.PW)  00003470 

C********* ******** ****** ************** ********************************** 000034 80 
C THIS  IS  THE  IMPLICIT  INTEGRATION  PROGRAM  AS  GIVEN  BY  C .W. GEAR  IN  00003490 

C "NUMERICAL  INITIAL  VALUE  PROBLEMS  IN  ORDINARY  DIFFERENTIAL  00003500 

C EQUATIONS".  PUB.  PR ENT  IC E-HALL . I NC . 1971.  00003510 

C •*•**•***• ******* *************** ***************** ****************  ******00003520 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS  SUBROUTINE  INTEGRATES  A SET  OF  N ORDINARY  DIFFERENTIAL  FIRST  00003530 
ORDER  EQUATIONS  OVER  ONE  STEP  OF  LENGTH  H AT  EACH  CALL.  H CAN  BE  00003540 
SPECIFIED  BY  THE  USER  FOR  EACH  STEP,  BUT  IT  MAY  BE  INCREASED  OR  00003550 
DECREASED  BY  OIFSUB  WITHIN  THE  RANGE  HMI N TO  HMA X IN  ORDER  TO  00003560 
ACHIEVE  AS  LARGE  A STEP  AS  POSSIBLE  WHILE  NOT  COMMITTING  A SINGLE  00003570 


STEP  ERROR  WHICH 
COMPONENT  OF  THE 


IS  LARGER  THAN  EPS  IN  THE  L-2  NORM,  WHERE  EACH 
ERROR  IS  DIVIDED  BY  THE  COMPONENTS  OF  YMAX. 


THIS  PROGRAM  REQUIRES  THREE  SUBROUTINES  NAMED 
DIFFUNIT, Y,DY) 

MATINVIPW ,N,M, J) 

PEDERVIT.Y.PW.M) 

THE  FIRST,  DIFFUN,  EVALUATES  THE  DERIVATIVES  OF  THE  DEPENDENT 
VARIABLES  STORED  IN  Yll, II  FOR  1 = 1 TO  N,  AND  STORES  THE 
DERIVATIVES  IN  THE  ARRAY  DY.  THE  SECOND  IS  CALLED  ONLY  IF  THE 
METHOD  FLAG  MF  IS  SET  TO  1 OR  2 FOR  STIFF  METHODS 
THE  N BY  N MATRIX  STORED  IN  THE  ARRAY  PwIM.M).  IF 
SUCCESFUL,  J SHOULD  BE  SET  TO  1 ~ 

PEDERV  IS  USED  ONLY  IF  MF  IS  1, 

DERIVATIVES  OF  THE  DIFFERENTIAL 
MF  PARAMETER. 


00003580 
00003590 
00003600 
00003610 
00003620 
00003630 
00003640 
00003650 
00003660 
00003670 
IT  MUST  INVERT00003680 
THE  INVERSION  100003690 
OTHERWISE  IT  SHOULD  BE  SET  TO  -1.00003700 


AND  COMPUTES  THE  PARTIAL 
EQUATIONS  AS  DESCRIBED  UNDER  THE 


THE  PROGRAM  USES  DOUBLE  PRECISION  ARITHMETIC 
POINT  VARIABLES  EXCEPT  THOSE  STARTING  WITH  P 
SINGLE  PRECISION  TO  SAVE  TIME  AND  SPACE. 


FOR  ALL  FLOAT ING 
THE  FORMER  ARE 


THE  TEMPORARY  STORAGE  SPACE  IS  PROVIDED  BY  THE  CALLER  IN  THE 
SINGLE  PRECISION  ARRAY  PW  AND  THE  COJBLE  PRECISION  ARRAY  SAVE. 

THE  ARRAY  PW  IS  USED  ONLY  TO  HOLD  THE  MATRIX  OF  THE  SAME  NAME, 

TO  HOLD  SEVERAL  ARRAYS.  THE  REGIONS  USED  ARE 

1.LE.J.LE.8  AND  l.LE.I.LE.N  IS  USED  TO  SAVE  THE 
VALUES  OF  Y IN  CASE  A STEP  HAS  TO  BE  REPEATED. 

IS  USED  MAINLY  TO  HOLD  THE  CORRECTION  TERMS  IN  THE 
CORRECTOR  LOGP. 

IS  USED  TO  SAVE  THE  VALUES  OF  THE  SUMS  OF  ALL  OF 


SAVE  IS  USED 
SAVE!  J . I > 

SAVE  1 9 , I ) 

SAVE!  10, I I 


00003710 
00003720 
00003730 
00003740 
00003  750 
00003760 
00003770 
00003780 
00003790 
00003800 
BUT00003810 
00003820 
00003830 
00003840 
00003850 
00003860 
T 0000  3 8 70 


THE 

THE 


CORRECTION  TERMS  IN  THE  PREVIOUS  STEP  AFTER  THEY  00003880 
HAVE  BEEN  ACCUMULATED  IN  THE  ARRAY  ERROR  IN  THE  00003890 
CURRENT  STEP.  THIS  ENABLES  THE  BACKWARDS  D I F FE RE NC 00003900 
OF  ERROR  TO  BE  FORMED.  IT  IS  USED  TO  ESTIMATE  THE  00003910 
STEP  SUE  FOR  ONE  ORDER  HIGHER  THAN  CURRENT.  00003920 

SAVEINl+1,1)  IS  USED  TO  STORE  THE  DERIVATIVES  WHEN  THEY  ARE  00003930 
COMPUTED  BY  DIFFJN.  IT  IS  ALSO  ACCESSED  AS  00003940 

SAVEIN2,!)  AS  A COMPLETE  ARRAY.  00003950 

SAVEIN5+I.II  HOLDS  THE  DERIVATIVES  DURING  JACOBIAN  E VAL UAT I QNS00003960 
IT  IS  REFERENCED  AS  SAVE!N6,l)  AS  A COMPLETE  ARR A Y . 0000 3 9 70 

00003980 

PARAMETERS  TO  THE  SU8R0UT INE  OIFSUB  HAVE 
FOLLOWING  MEANINGS.. 


N 


THE 

MAY 


NUMBER  OF  FIRST  ORDER  DIFFERENTIAL  EQUATIONS. 
3F  DECREASED  ON  LATER  CALLS  IF  THE  NUMBER  OF 


N 


00003990 

00004000 

00004010 

00004020 

00004030 
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T 

Y 


SAVE 

H 


HMI  N 

HR  AX 
EPS 

MF 


VMA  X 

ERROR 

KFLAG 


1TI0NS  REDUCES.  BUJ 

INCREASED  WITHOUT  CALLING 


ACT  IVE  EQUAl 


...  IT  MUST  NOT  BE 
W ITH  J ST  ART  = 0. 

THE  INDEPENDENT  VARIABLE. 

AN  8 8Y  N ARRAY  CONTAINING  THE  DEPENDENT  VARIABLES  AND 

THEIR  SCALED  DERIVATIVES.  Y ( J*l , I ) CONT A INS 

THE  J-TH  DERIVATIVE  OF  YIN  SCALED  BY 

H**J/FACTORIAUJ)  WHERE  H IS  THE  CURRENT 

STEP  SIZE.  ONLY  Tll.ll  NEED  BE  PROVIOED  BY 

THE  CALLING  PROGRAM  ON  THE  FIRST  ENTRY. 

IF  IT  IS  DESIRED  TO  INTERPOLATE  TO  NON  MESH  POINTS 
THESE  VALUES  CAN  BE  USED.  IF  THE  CURRENT  STEP  SIZE 
IS  H AND  THE  VALUE  AT  T t E IS  NEEDED,  FORM 
S = E/H,  AND  THEN  COMPUTE 
NO 

Y II  I I T ♦ E I = SUM  Yl J*l. I 1*S**J 

J=0 


A BLOCK  UF  AT  LEAST  12*N  FLOATING  POINT  LOCATIONS 
USED  BY  THE  SUBROUTINES. 

THE  STEP  SIZE  TO  BE  ATTEMPTED  ON  THE  NEXT  STEP. 

H MAY  BE  ADJUSTEO  UP  OR  DOUN  BY  THE  PROGRAM 
IN  ORDER  TO  ACHE  I V E AN  ECONOMICAL  INTEGRATION. 

HOWEVER,  IF  THE  H PROVIDED  BY  THE  USER  DOES 
NOT  CAUSE  A LARGER  ERROR  THAN  REQUESTED,  IT 
WILL  BE  USED.  TO  SAVE  COMPUTER  TIME,  THE  USER  IS 
ASVISED  TO  USE  A FAIRLY  SMALL  STEP  FOR  THE  FIRST 
CALL.  IT  WILL  BE  AUTOMATICALLY  INCREASED  LATER. 

THE  MINIMUM  STEP  SIZE  THAT  WILL  BE  USED  FOR  THE 
INTEGRATION.  NOTE  THAT  ON  STARTING  THIS  MUST 
BE  MUCH  SMALLER  THAN  THE  AVERAGE  H EXPECTED  SINCE 
A FIRST  ORDER  METHOD  IS  JSED  INITIALLY. 

THE  MAXIMUM  SIZE  TO  WHICH  THE  STEP  WILL  BE  INCREASED 
THE  ERROR  TEST  CONSTANT.  SINGLE  STEP  ERROR  ESTIMATES 
DIVIDED  BY  YMAXll)  MUST  BE  LESS  THAN  THIS 
IN  THE  EUCLIDEAN  NORM.  THE  STEP  AND/OR  ORDER  IS 
ADJUSTED  TO  ACHE  I V E THIS. 

THE  METHOD  INDICATOR.  THE  FOLLOWING  ARE  ALLOwEO. 

0 AN  ADAMS  PREDICTOR  CORRECTOR  IS  USED. 

1 A MULTI-STEP  METHOD  SUITABLE  FOR  STIFF  SYSTEMS  IS  USED. 
IT  WILL  ALSO  WORK  FOR  NON-STIFF  SYSTEMS.  HOWEVER 

THE  USER  MUST  PROVIDE  A SUBROUTINE  PEDERV  WHICH 
EVALUATES  THE  PARTIAL  DERIVATIVES  OF  THE  DIFFERENTIAL 
EQUATIONS  WITH  RESPECT  TO  THE  Y ' S.  THIS  IS  DONE  BY 
CALL  PEDERV <T,Y,PW,M) . PW  IS  AN  N BY  N ARRAY  WHICH 
MUST  BE  SET  TO  THE  PARTIAL  OF  THE  l-TH  EQUATION  WITH 
RESPECT  TO  THE  J DEPENDENT  VARIABLE  IN  PWU.J).  PW  IS 
ACTUALLY  STORED  IN  AN  M BY  M ARRAY  WHERE  M IS  THE 
VALUE  OF  N USED  ON  THE  FIRST  CALL  TO  THE  PROGRAM. 

2 THE  SAME  CASE  AS  1,  EXCEPT  THAT  THIS  SUBROUTINE 
COMPUTES  THE  PARTIAL  DERIVATIVES  BY  NUMERICAL 
DIFFERENCING  OF  THE  DERIVATIVES.  HENCE  PEDERV  IS  NOT 
CALLED. 

AN  ARRAY  OF  N LOCATIONS  WHICH  CONTAINS  THE  MAXIMUM 
OF  EACH  Y SEFN  SO  FAR.  IT  SHOULD  NORMALLY  BE  SET  TO 
1 IN  EACH  COMPONENT  BEFORE  THE  FIRST  ENTRY.  (SEE  THE 
DESCRIPTION  OF  EPS. I 

AN  ARRAY  OF  N ELEMENTS  WHICH  CONTAINS  THE  ESTIMATED 
ONE  STEP  ERROR  IN  EACH  COMPONENT. 

A COMPLETION  CODE  WITH  THE  FOLLOWING  MEANINGS.. 


00004040 
00004D5D 
00004060 
00004070 
00004080 
00004090 
00004100 
00004110 
00004120 
000041 30 
00004140 
00004150 
00004160 
00004170 
00004180 
00004190 
00004200 
00004210 
00004220 
00004230 
00004240 
00004250 
00004260 
000042  70 
00004280 
00004290 
00004300 
00004310 
00004  320 
00004330 
00004340 
00004350 
00004360 
00004370 
00004380 
00004390 
00004400 
000044 10 
00004420 
00004430 
00004440 
00004450 
00004460 
000044  70 
00004480 
00004490 
00004500 
000045 10 
00004520 
000045  30 
00004540 
00004550 
00004560 
OC  004  5 70 
00004580 
00004590 
00004600 
00004610 
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c 

c 

c 
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c 
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♦ 1 
-1 


_ _ SUCCESSFUL. 

THE  STfP  WAS  TAKEN  WITH  H=HMIN,  BUT 


THE 


TO 


THE  STFP  WAS 

THE  STFP  WAS  I ak  cn  min  n=r,mn<,  ot 
REQUESTED  ERROR  WAS  NOT  ACHIEVED. 

THE  MAXIMUM  ORDER  SPECIFIED  WAS  FOUNO 
BE  TOO  LARGE. 

CORRECTOR  CONVERGENCE  COULD  NOT  BE 
ACHIEVED  FOR  H.GT.HMIN. 

THE  REQUESTED  ERROR  IS  SMALLER  THAN  CAN 
BE  HANDLED  FOR  THIS  PROBLEM. 

JSTART  AN  INPUT  INDICATOR  WITH  THE  FOLLOWING  MEANINGS.. 

-1  REPEAT  THE  LAST  STEP  WITH  A NEW  H 

PERFORM  THE  FIRST  STEP.  THE  FIRST  STEP 
MUST  BE  DONE  WITH  THIS  VALUE  OF  JSTART 
SO  THAT  THE  SUBROUTINE  CAN  INITIALIZE 
ITSELF. 

TAKE  A NEW  STEP  CONTINUING  FROM  THE  LAST. 

JSTART  IS  SET  TO  NQ,  THE  CURRENT  ORDER  OF  THE  METHOD 
AT  EXIT.  NQ  IS  ALSO  THE  ORDER  OF  THE  MAXIMUM 
DERIVATIVE  AVAILABLE. 

THE  MAXIMUM  DERIVATIVE  THAT  SHOULD  BE  USED  IN  THE 
METHOD.  SINCE  THE  ORDER  IS  EQUAL  TO  THE  HIGHEST 
DERIVATIVE  USED,  ThIS  RESTRICTS  THE  ORDER.  IT  MUST 
8 E LESS  THAN  8 FOR  ADAMS  AND  7 FOR  STIFF  METHODS. 

A BLOCK  OF  AT  LEAST  N**2  FLOATING  POINT  LOCATIONS. 


-3 

-4 


0 


♦ I 


MAXDER 


*****************  ***********  * **,*  *******  ****************** 


IMPLICIT  REAL*d( A-H, Q-Z) 
RE  AL  *4  AM  AX  1 


DIMENSION  Y<  8,  1 ) .YMAXl  1 ) , 5 A VE I 10 , 1 ) .ERROR!  ll.PWIll, 

I AI81 .PERTST (7 ,2 ,31 

THE  COEFFICIENTS  IN  PEP.TST  ARE  USED  IN  SELECTING  THE  STEP  AND 
ORDER,  THEREFORE  ONLY  ABOUT  ONE  PERCENT  ACCURACY  IS  NEEDED. 

OAT A PERTST/2. 0,4. 5, 7. 333,10. 42, l 3. 7,17. 15,1.0, 

1 2. 0,1 2. 0,24. 0,37. 89, 53. 33, 70. 08, 8 7. 9 7, 

2 3. 0,6.  0,9.  167, 12. 5, 15. 98, 1. 0,1.0, 

3 12. 0,24.0, 37. 89, 53. 33,  70.  08, 87. 97,1.  0, 

4 l.  ,1  . , 0.5,0. IBB  7, 0.04133. 0.00826 7,  1.0, 

5 1.0,  1.0, 2. 0,1. 0,-3157, .074 07  , .0139/ 

DATA  A I 2 I / — 1 .00  + 0/ 

I RE  T = l 
KFL AG=  1 

IF  I JSTART .L E. 01  GO  TO  140 

BEGIN  BY  SAVING  INFORMATION  FOR  POSSIBLE  RESTARTS  AND  CLANGING 
H BY  THE  FACTOR  R IF  THE  CALLER  HAS  CHANGED  H.  ALL  VARIABLES 
DEPENDENT  ON  H MUST  ALSO  BE  CHANGEO. 

E IS  A COMPARISON  FOR  ERRORS  OF  THE  CURRENT  ORDER  NQ.  EUP  'S 
TO  TEST  FOR  INCREASING  THE  ORDER,  EOWN  FOR  DECREASING  THE  ORDER. 
HNEW  IS  THE  STEP  SIZE  THAT  mAS  USED  ON  THE  LAST  CALL. 

100  00  1 10  I =1  , N 

DO  110  J*l,  K 
110  SAVE (J, I )=  Y (J,  I I 


00  004o  20 
00004630 
00004640 
00004650 
00004660 
00004670 
000046  80 
00004690 
30004700 
00004710 
00004720 
00004730 
00004740 
00004750 
00004760 
00004770 
00304780 
00004790 
00004800 
000048  13 
0C004820 
00004830 
00004840 
00004850 
**00004863 
00004870 
00004880 
00004890 
00004900 
00004910 
00004920 
00004930 
00004940 
00004950 
00004960 
000049  70 
00004980 
00004990 
00005000 
000050 10 
00005020 
00005030 
00005040 
00005050 
00005060 
00005070 
00005080 
00005090 
00005100 
000051  10 
00005120 
00005130 
00005140 
00005150 
00005160 
00005170 
00005180 
00005190 


! 


1-48 


RELEASE  2.0 


D1FSUB 


DATE  = 77014 


15/09/38 


120 

130 


140 

C 

C 

C 

c 


150 


C 

c 

c 

160 


C 

c 

c 

c 

c 

c 

c 

c 

170 


180 

190 

C 

C 

c 

c 

c 

c 

c 


H0L0  = HNEW 

I F ( H. E Q. HOLD ) GO  TO  130 

RACUM=  H/HOLD 

I RET  1 = 1 

GO  TO  750 

NUOLD=NQ 

TOLD  = T 

RACUM=  1 

IF  ( JSTART.GT.OI  GO  TO  250 
GO  TO  170 

IF<  JSTART.EQ.-l)  GO  TO  160 

ON  THE  FIRST  CALL.  THE  OROER  IS  SET  TO  1 AND  THE  INITIAL 
DERIVATIVES  ARE  CALCULATED. 

N0=  1 
N3=  N 
Nl=  N* 10 
N2=  N1+  1 
N4=  N**2 
N5=  N 1 + N 
N6=  N5  ♦ 1 

CALL  0IFFON1  T . Y , SA VE ( N2 , l ) ) 

DO  150  1=1.  N 

Yl 2 . I ) = S A VE ( N 1 ♦ 1,1)  *H 

HN  EW  = H 

K=2 

GO  TO  100 

REPEAT  LAST  STEP  BY  RESTORING  SAVED  INFORMATION. 

IMNO.EQ.NQOLD)  J S TAR  T=  1 
T = TOLD 
NU=  NOOLD 
K=  NO  ♦ 1 
GO  TO  120 

SET  THE  COEFFICIENTS  THAT  OETERMINF  THE  ORDER  AND  THE  METHOD 
TYPE.  CHECK  FOR  EXCESSIVE  ORDER.  THE  LAST  TOO  STATEMENTS  OF 
THIS  SECTION  SET  I WE  V A L.  G T . 0 IF  PW  IS  TO  BE  RE-EVALUATED 
BECAUSE  OF  THE  ORDER  CHANGE,  AND  THE  REPEAT  THE  INTEGRATION 
STEP  IF  IT  HAS  NOT  YET  BEEN  DONE  ( IRET  = 1)  OR  SKIP  TO  A FINAL 
SCALING  BEFORE  EXIT  IF  IT  HAS  BEEN  COMPLETED  (IRET  = 2). 

IF(MF.EQ.O)  GO  TO  180 
IFINQ.GT .6)  GO  TO  190 
GOTO  (221,222,223,224, 225, 226) , NO 
IFINQ.GT. 7)  GO  TO  190 

GO  TO  (211.212,213,214,215,216,217) , NO 

KFL  AG=-2 

RETURN 

THE  FOLLOWING  COEFFICIENTS  SHOULD  BE  DEFINED  TO  THE  MAXIMUM 
ACCURACY  PERMITTED  BY  THE  MACHINE.  THEY  ARE,  IN  THE  ORDER  USED. 
-1, -1/2, -1/2, -5/1 2, -3/4, -1/6, -3/8, -11/12, -1/3,  -1/24 
-251/ 720,-25/24,-35/  72, -5/ 48, -1/120, -95/ 288,-1  37/120, -5/8 
-17/96.-1/40 ,- 1/ 72 0.-190B 7/ 60480 ,-49/40,-203/270 ,-49/ 192 
-7/144,-7/1440,- 1/5040,- 1,-2/ 3,-1 /3.-6/1 1 ,-6/11,-1/11,-12/25 


00005200 
00005210 
00005220 
00005230 
00005240 
00005250 
00005260 
00005270 
000  052  oO 
00005290 
00005300 
00005310 
00005320 
00005330 
00005340 
00005350 
00005360 
00005370 
00005380 
00005390 
00005400 
00005410 
00005420 
00005430 
00005440 
00005450 
00005460 
00005470 
00005480 
00005490 
00005500 
00005510 
00005520 
000055  30 
00005540 
00005550 
00005560 
00005570 
00005580 
00005590 
00005600 
00005610 
00 005620 
000056  30 
00005640 
30005650 
00005660 
00005c*  70 
00005680 
00005690 
00005700 
00005  7 10 
00005720 
00005730 
00005740 
00005750 
00005  760 
00005770 


! I -49 


d 


OoOcg 


RELEASE  2.0 


DIFSUB 


OATE  = 77015 


15/09/38 


-7/10.-1/5.-1/50 .-120/279,-225/275.-85/279,-1 
-180/551.-58/63.-15/36,-25/252,-3/252,-1/1765 


15/2  75,-1/275 


11 

212 

213 

215 

215 

216 

217 

221 

222 

223 

225 

225 


A ( 1 » = -1  .0  0*0 

GO  TO  230 

All)—  - . 5 000  00  000  <-0 
A 1 3 1 = -.5000000000*0 
GO  TO  230 

A ( 1 I = -.91666666666666670*0 
A ( 3 ) = -.75000000000*0 
A ( 5 1 = -.16666066666666670*0 
GO  TO  230 

A 1 1 I = — 0.375  000000D*0 
A(3)=  -0.91666666666666670*0 
A 1 5 I = -0  .33333333333333330*0 
A ( 5 ) = -0.051666666666666670+0 
GO  TO  230 

A ( 1 1 * -0.358611111111111  10  + 0 
A 1 3 ) = -1.05166666666666670  + 0 
A l 5 1 = -0.58611111111111110*0 
A15)  = —0.10516666666666670+0 
At  6) =-0. 0083  33  33  33  33  3 33  3 330+0 
GO  TO  230 

A ( l ) = -0.32986111111111110*0 
A ( 3) = -1.15166666666666670+0 
A 1 5 I = -0.6250000000+0 
A < 5 1 = -0.17708333333333330*0 
A ( 61 = -0.0250000000*0 
A < 7 I = -0.001  3888888888888890*0 
GO  TO  230 

A<11  = -0.31559193121693120*0 
A13>  = -l  .235000000D+0 
A( 5) =-0. 75185185185185190*0 
A ( 5 1 = -0.25520833333333330*0 
A(61 =-0.0586 11111111111110*0 
A ( 7) =-0.0058611111111111110*0 
A ( 8 )=  -0.00019851269851269850*0 
GO  TO  230 

At  11=  -1.0000000000+0 
GO  TO  230 

A(l)=  -0.66666666666666670*0 
A1 3) = -0.33333333333333330+0 
GO  TO  230 

A( II =-0. 55  55  5555  55555555  0*0 
A ( 3 ) = A l 11 

A ( 5 1 = — 0 . 0909  0909  CT909  0909  10*0 
GO  TO  230 

A ( 1 1=  -0.580000000000*0 
A 13  I =-0. 7000000000000*0 
A( 51 =-0. 2000000000000*0 
A<  5 I = -0.0200000000000*0 


GO  TO  230 

A 1 1I--0.  5379562053795620*0 
A(3I= -0.82116 78832116 78 80+0 
A 151  =-0.31021897810218980*0 
A(  51 =-0.  0557  5552  555  755  52  60*  0 
A<6 I =-0.0036 4963 50 365963 5050*0 
GO  TO  230 


00005780 
00005790 
00005800 
000058 10 
00005820 
00005830 
00005850 
00005650 
00005860 
00005670 
00005880 
00005890 
00005900 
00005910 
00005923 
00005930 
30005950 
00005950 
00005960 
00005970 
000059  80 
00  005  99 J 
00006000 
00006010 
00006020 
00006030 
00006050 
00006050 
00006060 
00006070 
00006080 
00006090 
00006100 
000061 13 
00006120 
000061 30 
00006 153 
0C0061 50 
00006160 
000061 70 
000061 80 
00006190 
00006200 
000062  10 
00006223 
00006230 
00006250 
00006250 
000062  60 
00006270 
00006280 
00006290 
00006300 
00006313 
00006320 
00006330 
00006350 
00006350 


I I -SO 


RELEASE  2-0 


01 FSUB 


DATE 


7 70  14 


15/09 /3b 


226 


230 


240 

C 

c 

c 

c 

c 

250 


260 

C 

C 

c 

c 

c 

c 

c 

c 

c 

270 


C 

C 

C 

c 

c 

c 


280 

C 

c 


All) 
A ( 3 
A ( 4 ) 
A(  5) 
A (6  ) 
A ( 7) 


0.  4081  63265306122  50*- 0 
-0 .92063492063492060+0 
0. 41666666666666670+0 
0. 09920634920634920+0 
0.01190476190476190+0 
0. 0005  66  89  3 4240  3 62  d2  0 + 0 


K=NQ+1 

I00U8=K 

MTYP-- ( 4- ME  ) /2 

ENU2  = -5/DFL0AT (N0+ 1) 

ENU3=.5/DFLOAT (NQ+2) 

ENQ1=0. 5/OFLQA i i NO) 

PE  PS  H=  EP  S 

EUP=  < PERTSTl NO, MTYP.2 ) *PfcPS  F)**2 
E=(PERTST(NG,MTYP.l)*PEPSH)**2 
EOWN= (PERT  ST  4N0, MTYP.3 )*PEPSH)**2 
IFIEDWN.EQ.OIGO  TO  780 
8ND=EPS*EN03/DFL0A TIN) 

IWEV AL  = MF 

GO  TO  (250,680)  , I RE  T 


THIS  SECTION  COMPUTES  THE  PREDICTED  VALUES  BY  EFFECTIVELY 
MULTIPLYING  THE  SAVEO  INFORMATION  BY  THE  PASCAL  TRIANGLE 

MATRIX. 


T=  T + H 

on  260  J = 2,K 
00  260  J l = J , K 
J2  = K- J 1+ J- 1 
00  260  I = 1 , N 

Y 4 J2 ,l)=Y( J2 ,1 )+Y( J2+1  , I ) 

UP  TO  3 CORRECTOR  ITERATIONS  ARE  TAKEN.  CONVERGENCE  IS  TESTED 
BY  REOUIRING  CHANGES  TO  BE  LESS  THAN  BN 0 WHICH  IS  DEPENDENT  ON 
THE  ERROR  TEST  CONSTANT. 

THE  SUM  OF  THE  CORRECTIONS  IS  ACCUMULATED  IN  THE  ARRAY 
ERRORIII.  IT  IS  EUUAL  TO  THE  K-TH  DERIVATIVE  OF  Y MULTIPLIED 
BY  H**K /< F AC  TOR  I AL I K-l) »AI K ) ) , AND  IS  THEREFORE  PROPORTIONAL 
TO  THE  ACTUAL  ERRORS  TO  THE  LOWEST  POWER  OF  H PRESENT..  (H**M 


00  270  1 = 1, N 
ERROR l I ) =0  .0 
DO  430  L=1  ,3 

CALL  OIFFUNI  T, Y,  SAVEIN2,  1) ) 

IF  THERE  HAS  BEEN  A CHANGE  OF  OROER  OR  THERE  HAS  BEEN  TROUBLE 
WITH  CONVERGENCE,  Pw  IS  RE-E  VALUA  TfcO  PRIOR  TO  STARTING  THE 
CORRECTOR  ITERATION  IN  THE  CASE  OF  STIFF  METHODS.  IWEVAL  IS 
THEN  SET  TO  -1  AS  AN  INDICATOR  THAT  IT  HAS  BEEN  DONE. 


IFIIWEVAL.LT.il  GO  TO  350 
IFIMF.E0.2)  GO  TO  310 
CALL  PEDERVI T,Y,Pw,N3) 

R = A I 1 )*H 


00  280  1=1  , N4 

PW(  I )=PW<  I ) *R 


A DO  THE  IDENTITY  MATRIX  TO  THE  JACOBIAN  AND  INVERT  TO  GET  PW . 


03006363 
000063  70 
00006380 
00006390 
00C06400 
000064 10 
00006420 
00006430 
00006440 
00006450 
00006460 
00006470 
00006480 
00006490 
00006500 
00006510 
00006520 
00006530 
00006540 
0 000  6 5 50 
000J6560 
00006570 
00006580 
00006590 
00006600 
000066 10 
00006620 
00006630 
00006640 
00006650 
00006660 
000066  70 
00006660 
00006690 
00006700 
00006710 
00006720 
0000673 0 
00006  740 
00006750 
00006760 
000067  70 
00006780 
00006790 
00006800 
00006810 
33006820 
00006830 
00006840 
00006850 
00006860 
000068  70 
00006880 
00006890 
00006900 
00006910 
00006920 
00006930 


1-51 


RELEASE  2.0 


01 F SUB 


DATE 


7701'* 


15/09/38 


C 

290  00  300  I =1  ,N 

300  PWt  I *1  N3*l )-N3l  = l. O+PWI l *( N3*l I-N3  I 
I WEV  AL  = - 1 

CALL  MATINVIPW,N,N3, J1 I 
IF(Jl.GT.O)  GO  TO  350 
GO  TO  440 
C 

C EVALUATE  THE  JACOBIAN  INTO  PW  BY  NUMERICAL  0 1 FF E RENC I NG.  R IS  THE 
C CHANGE  MADE  TO  THE  ELEMENT  OF  Y.  IT  IS  EPS  RELATIVE  TO  Y WITH 
C A MINIMUM  OF  EPS**2. 

C 

310  00  320  1=1, N 

320  SAVE  I 9, 1 I = VI 1 , 1 1 
DO  3 AO  J=l,N 

R=EPS*OMAXll EPS, DABS! SAVE! 9.J) i I 
Y I 1 , J ) = Y 1 1 ,J)*R 
0=  A ( l ) *H/R 

CALL  0 IFFUNIT , Y .SAVE! N6, 1 II 
| 00  330  I =1  , N 

330  PWt I ♦ ! J-  1)*N3)=( SAVE! N5*l, ll-SAVE IN  1*1 ,1) 1*0 
340  Y(1 , JJ=SAVEI9, J) 

GO  TO  290 

350  IF  (MF.NE.OI  GO  TO  370 
00  360  I =1  ,N 

360  SAVE  I 9, I )=YI  2, I I -SAVE  I Nl  + l , l ) «H 
GO  TO  410 

370  DO  380  1=  1,N 

380  SAVEIN5H  , It  =Y12,  I I-SAVEINIH  ,1  I *H 
00  400  1=1, N 
0=0.0 

00  390  J=1,N 

390  0 = 0* PW I I ♦( J-l I *N3 ) *SAVEIN5*J, It 

400  SAVE (9, 1 1=0 
410  NT=N 
C 

C CORRECT  ANO  SEE  IF  ALL  CHANGES  ARE  LESS  THAN  BN E RELATIVE  TO  YMAX. 
C IF  SO,  THE  CORRECTOR  IS  SAID  TO  HAVE  CONVERGEO. 

C 

00  420  1=1  , N 

YI1,  I ) = Y ( 1 , I )*Al  1 1 *$ AVE ( 9, I I 
YI2 , I )=Y  12  ,1  I-SAVEI9,  l I 
ERROR  I I ) =ERROR! I 1 A SAVE (9, 1 I 

IF  I DABS  I SAVE  19, I ) ) . L E . ( 8N0* YM AX < I ) ) I N T =N  T-  1 
420  CONTINUE 

IF  INT.LE.OI  GO  TO  490 
430  CONTINUE 
C 

C THE  CORRECTOR  ITERATION  FAILED  TO  CONVERGE  IN  3 TRIES.  VARIOJS 
C POSSIBILITIES  ARE  CHECKEO  FOR.  IF  H IS  ALREADY  HMI N AND 
C THIS  IS  EITHER  ADAMS  METHOD  OR  THE  STIFF  METHOD  IN  WHICH  THE 
C MATRIX  PW  HAS  ALREAOY  BEEN  TE - E VAL U AT E 0 , A NO  CONVERGENCE  EXIT 
C IS  TAKEN.  OTHERWISE  THE  MATRIX  PW  IS  RE-EVALUATED  AND/OR  THE 
C STEP  IS  REDUCED  TO  TRY  AND  GET  CONVERGENCE. 

C 

440  T=TOLD 

IF! I H.LE.I HMIN*1 .00001 » I .AND. I I IWEV AL-MTYP).LT .- I I I GO  TO  450 
I F ( (MF.EO.O) .OR. I I WE VAL .NE . 01 ) R AC UM  = R AC UM* 0. 2500 


00005940 
00006950 
00006960 
00006970 
00006980 
00006990 
00007000 
00007010 
00007020 
00007030 
00007040 
00007050 
00007060 
00007070 
00007080 
00007090 
00007 100 
0000711 0 
00007120 
000071 30 
00007140 
00007150 
00007160 
000071 70 
000071 80 
00007190 
00007200 
00007210 
00007220 
00007230 
00007240 
00007250 
00007260 
00007270 
00007280 
00007290 
0030  7300 
00007310 
00007320 
00007330 
00007340 
00007350 
00007360 
00007370 
00007380 
00007390 
00007400 
00007410 
00007420 
00007430 
00007440 
00007450 
00007460 
00007470 
30007480 
00007490 
00007500 
00007510 


DIF  SUB 


DATE 


7 7014 


15/09/38 


RELEASE  Z.C 


imw 

GO  TO  750 
460  KFLAG=-3 
470  00  480  I=ltN 

DO  480  J=1.K 
480  Y 4 J , I 1 =S AVEl  J. I » 
H=H0LD 
NO=NOOLO 
JST  ART  = NQ 
RETURN 


C 

C 

c 

c 


THE  CORRECTOR  CONVERGED  ANO  CONTROL  IS  PASSED  TO  STATEMENT  520 
IF  THE  ERROR  TEST  IS  O.K.,  AND  TO  540  OTHERWISE. 

IF  THE  STEP  IS  0.  K.  IT  IS  ACCEPTED.  IF  I DOUB  HAS  BEEN  REOUCEO 

C TO  ONE,  A TEST  IS  MADE  TO  SEE  IF  THE  STEP  CAN  BE  IN  REASED 

C AT  THE  CURRENT  ORDER  OR  BY  GOING  TO  ONE  HIGHER  OR  ONE  LOWER. 

c such  a change  is  only  made  if  the  step  can  be  increased  by  at 

C LEAST  l.l.  IF  NO  CHANGE  IS  POSSIBLE  ID0UB  IS  SET  TO  10  TO 

C PREVENT  FUTHER  TESTING  FOR  10  STEPS. 

C IF  A CHANGE  IS  POSSIBLE,  IT  IS  MADE  ANO  IDOUb  IS  SET  TO 
C NQ  ♦ 1 TO  PREVENT  FURTHER  TESTING  FOR  THAT  NUMBER  OF  STEPS. 

C IF  THE  ERROR  WAS  TOO  LARGE,  THE  OPTIMUM  STEP  SIZE  FOR  THIS  OR 
C LOWER  OROER  IS  COMPUTED,  ANO  THE  STEP  RETRIED.  IF  IT  SHOULD 
C FAIL  TWICE  MORE  IT  IS  AN  INDICATION  THAT  THE  DERIVATIVES  THAT 
HAVE  ACCUMULATED  IN  THE  Y ARRAY  HAVE  ERRORS  OF  THE  WRONG  OROER 
SO  THE  FIRST  DERIVATIVES  ARE  RECOMPUTED  ANO  THE  ORDER  IS  SET 
TO  1. 


C 
C 
C 
C 

490 


500 


0=0.0 

DO  500  I =1 • N 

0=0* (ERROR (I l/YMAXI I >1**2 
I WE V AL  = 0 
IF  (D.GT.E)  GO  TO  540 

tFU.LT.3l  GO  TO  520  + 

COMPLETE  THE  CORRECTION  OF  THE  HIGHER  ORDER  DERIVATIVES  AFTER  A 
SUCCESFUL  STEP. 


DO  510  J=3,K 
DO  510  I = l.N 

510  Y(  J,  I ) = Y ( J,  I )«-A(  J ) *ERROR(  I ) 

520  KFLAG=*1 
HNEW=H 

IF ( IDOUB. LE. 1 ) GO  TO  550 
I DOUB= IDOUB  -1 
IF!  IOOUB .GT.  1 ) GO  TO  700 
DO  530  I =1  ,N 

530  SAVE!  10,  I ) =E  RR  OR  ( I ) 

GO  TO  700 

C REDUCE  THE  FAILURE  FLAG  COUNT  TO  CHECK  FOR  MULTIPLE  FAILURES. 

C RESTORE  T TO  ITS  ORIGINAL  VALUE  AND  TRY  AGAIN  UNLESS  THERE  HAVE 
C THREE  FAILURES.  IN  THAT  CASE  THE  DERIVATIVES  ARE  ASSUMED  TO  HAVE 
C ACCUMULATED  ERRORS  SO  A RESTART  FROM  THE  CURRENT  VALUES  OF  Y IS 

C TRIED. 

C 

540  KFLAG*KFLAG-2 


00007520 
00007530 
00007540 
00007550 
00007560 
00007570 
00007580 
00007590 
00007600 
00007610 
00007620 
00007630 
00007640 
00007650 
00007660 
00007670 
00007680 
00007690 
00007700 
000077 10 
00007720 
00007730 
00007  740 
00007750 
00007760 
00007770 
00007760 
00007790 
00007800 
00007610 
00007820 
00007830 
00007840 
00007850 
00007860 
00007870 
00007880 
00007890 
00007900 
00007910 
00007920 
00007930 
00007940 
00007950 
00007960 
00007970 
00007980 
0000  7990 
00008000 
00008010 
00008020 
00008030 
00008040 
00008050 
00008060 
00008070 
00008080 
00008090 


I I -S  3 


RELEASE  2.0  DIFSUB  DATE  = 77014  15/09/38 

IF  ( H.LE .( HMIN*l. 00001) I GO  TO  740 
T=TOCO 

IF  CKFLAG.LE.-5)  GO  TO  720 
C 

C PR1.PR2.AN0  PR3  WILL  CONTAIN  THE  AMOUNTS  BY  WHICH  I HE  STEP  SIZE 
C SHOULD  BE  DIVIDED  AT  ORDER  ONE  LOWER,  AT  THIS  ORDER,  AND  AT  ORDER 
C ONE  HIGHER  RESPECTIVELY. 

C 

550  PR2  =IO/E) **ENQ2*1.2 
PR3= 1 . E+20 

IF  I ( NQ.GE.MAXDERI .OR. (KFLAG.LE.-1 ))  GO  TO  570 
D=0.  0 

DO  560  1 = 1, N 

560  0=DM  (ERROR!  I ) -S AV El  1 0 , 111 / YM AX  ( I I > **2 

PR  3=  ( 0 / E UP ) **E  NO  3*  1 . 4 
570  PRl=l.E+20 

IF(NQ.LE.l)  GO  TO  590 

0=0.0 

DO  580  1=1, N 

580  D=DMY(K»I  1/YMAXl  I 11**2 
PR1= (D/EDwN)**ENQl*1.3 
590  CONTINUE 

IF  (PR2.LE.PR3)  GO  TO  650 
IF  (PR3.LT. PRl)  GO  TO  660 
600  R=1.0/AMAXl(PRl,l.E-4) 

NEW0=NQ-  1 
610  I DOU  B= 10 

IF( (KFLAG.EO.l 1 . AND. (R.LT.( 1.1)1  I GO  TO  700 
1F1NEW0.LE.N0)  GO  TO  630 
C 

C COMPUTE  ONE  ADDITIONAL  SCALED  DERIVATIVE  IF  ORDER  IS  INCREASED. 

C 

DO  620  1=1, N 

620  YI  NE  WO*  1 ,1  I =E  RROR  ( I ) *A(  M/DFLOAT  (K) 

630  K=NEW0+1 

IFIKFLAG.EO.l)  GO  TO  670 
R AC  UM  = R AC  UM*  R 
IRET 1=  3 
GO  TO  750 

640  IFINEWO.EO.NQ)  GO  TO  250 
N0=  NE  WO 
GO  TO  170 

650  IF  (PR2.GT.PR1)  GO  TO  600 
NEW  0 = NQ 

R=1.0/AMAX l ( PR2 , 1 . E-4) 

GO  TO  610 

660  R=  1 . 0/ AM  AX  1 1 PR  3,  l.E— 41 
NEWO=NO» 1 
GO  TO  610 
670  IRET  = 2 

R = DMI N1 I R.HMAX/OABSIH)  I 

H = H*R 

HNEW=H 

IF  (NO.Ea.NEWO)  GO  TO  680 
NO=NEWQ 
GO  TO  170 
680  R 1= 1 . 0 

00  690  J=  2 » K 


00QO8 l 00 
00008110 
00008120 
000081 30 
00008140 
00008150 
00008160 
000081 70 
00008180 
000081 90 
00008200 
00008213 
00008220 
00008230 
00008240 
00008250 
00008260 
00008273 
00008280 
00008290 
00008300 
00008310 
00008320 
30008330 
00008340 
00008  353 
00008360 
00008370 
00008380 
00008390 
00008400 
00008413 
00008420 
00008430 
00008443 
00008450 
00008460 
00008473 
00008480 
00008490 
03008500 
00008510 
00008520 
00  00(35  30 
00008540 
00008550 
00008560 
00008570 
00008580 
00008590 
00008600 
00008610 
00008620 
00008630 
00008640 
00008650 
00006660 
000086  70 


I I -54 


RELEASE  2.0 


R1=R1*R 
00  690  1 =1 ,N 
690  rij, i ) = r ( j , i i*ri 

IDOUB-K 

700  00  710  1 =1 (N 

710  YMAXI  1 ) = ON  AX  II  YMA  X l l 1 , 0 AB  S 1 Y(  1.  I ) 1 1 

jstart=no 

RE  TURN 

720  IFlNQ.EO.il  GU  TO  780 

CALL  D I F FU  NI  T .Y»SAVE<N2. 1)1 

R=H/HOLD 

DO  7 JO  l = 1 , N 

Y(1 , I 1 =S AVEl  1 , I) 

SAVE (2,1  )=HOLD*SAVE<  Nl  + I ,1) 

730  Y12.  I |=SAVE12,  1 1 *R 

N0  = 1 
KFL  AG= 1 
GO  TO  170 
790  KFLAG=-1 
hnEw=h 
JSTART=NO 
RETURN 

C THIS  SECTION  SCALES  ALL  VARIABLES  CONNECTED  WITH  H AND  RETURNS 
C TO  THE  ENTERING  SECTION. 

C 

750  RACUM  = 0MAX11 OABSIHMI N/HOLO)  ,RACU*1 

R ACUM=  ON  I N 1 ( RACUM,  DABS  I HMAX/HOL J 1 1 
Rl=l .0 

00  760  J = 2 « K 
R 1 = R 1*RACUM 
DO  760  I =1  ,N 

760  Y( J,  I 1 = S A VE I J, I I *R l 

h=hold*racum 
DO  770  1=1, N 
770  Y(  1,  I 1 = S A V E I l,  I 1 

I DOUB=K 

GO  TO  ( 130.250,6501 , IRET1 
780  KFL  AG= -5 

GO  TO  570 
END 


mmt 

00008700 
00008710 
00008720 
00008730 
00008750 
00008750 
00008760 
00 008770 
00008780 
00008790 
00008800 
00008810 
00008820 
00008830 
00008850 
00008850 
00008860 
00008870 
00008880 
00008890 
00008900 
00008910 
00008920 
000089  30 
00008950 
00008950 
00008960 
00008970 
00008980 
00008990 
00009000 
000090 10 
00009020 
000090  30 
00009050 
00009050 
0C009060 
00009070 


I l-SS 


5 I F I LI) 

Security  Classification 


DD  1473 


UNCLASS  If  II  D 

Security  Classification 


Operational  Severity 
System  Performance 
Rel  ief  Valves 
Contaminant  Degradation 
Dynamic  Behavior 
Data  Acquisition 
Data  Analysis 


INSTRUCTIONS 

1.  ORIGINATING  ACTIVITY:  Enter  the  name  and  address  imposei 

of  the  contractor,  subcontractor,  grantee,  Department  of  De-  such  as 

fense  activity  or  other  organization  (corporate  author)  Issuing  ^ ^ 

the  report. 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over-  (2) 

all  security  classification  of  the  report.  Indicate  whether 
“Restricted  Data”  is  included.  Marking  is  to  be  in  accord- 
ance with  appropriate  security  regulations.  (3) 

2b.  GROUP:  Automatic  downgrading  is  specified  in  DoD  Di- 

rective 5200.  10  aqd  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3 and  Group  4 as  author-  ^ 

ized. 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  “Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC." 

(2)  “Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  “U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


3.  REPORT  TITLE.  Enter  the  complete  report  title  in  ail 
capital  letters.  Titles  in  all  casea  should  be  unclassified. 

If  a meaningful  title  cannot  be  selected  without  classifies* 
tion,  show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 

report,  e.g. , interim,  progress,  summary,  annual,  or  firm]. 

Give  the  inclusive  dates  when  a specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  nume(s)  of  authors)  as  shown  on 
or  in  the  report.  Entei  last  name,  first  name,  middle  initial. 

If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  O;«thor  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report  as  day, 
month,  year;  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  l.  e.,  enter  the 
number  of  pages  containing  information. 

76.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 

references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 

the  applicable  number  of  ’he  contract  or  grant  under  which 
the  report  was  written. 

86,  8c,  & 8<7.  PROJECT  NUMBER.  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a  ORIGINATOR’S  REPORT  NUMBER(S):  Enter  the  offi- 

cial report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

96.  OTHER  REPORT  NUMBER! S):  If  the  report  has  been 
assigned  any  other  report  numbers  (either  h\  the  originator 
or  by  the  sponsor),  also  enter  this  number!  s). 

10.  AVAILABILITY  LIMITATION  NOTICES:  Enter  any  lim- 
itations on  further  dissemination  of  the  report,  other  than  those 


(4)  “U.  S.  military  agencies  may  obtain  copies  of  this 

report  directly  from  DDC  Other  qualified  users 
shall  request  through 


(5)  “All  distribution  of  this  report  is  controlled.  Qual- 
ified DDC  users  shall  request  through 


If  the  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi- 
cate this  fact  and  enter  the  price,  if  known. 

U,  SUPPLEMENTARY  NOTES:  Use  for  additional  explana- 

tory notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  ''pay- 
inf, } for)  the  research  and  development.  Include  address. 

15  ATTRACT  Enter  an  abstract  giving  a brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re- 
port If  additional  space  is  inquired,  a continuation  sheet  shajl 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in 
formation  in  the  paragraph,  represented  as  (TSj  (Si  (C).  or  fU)  I 

There  is  no  limitation  on  the  length  of  the  abstract  How-  i 
ever,  the  suggested  length  is  from  150  to  225  words 

14  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a -eport  and  may  be  used  as  ^ 
index  entries  for  cataloging  the  rcpoi  Key  words  must  be 
selected  so  that  no  security  classification  is  required  Identi- 
fiers. such  as  equipment  model  designation,  trade  name  military 
project  < 'dc  name,  geographic  location,  mav  be  used  as  key 
words  but  will  he  followed  by  an  indication  of  teenmeal  con- 
text. The  assignment  of  links,  rules  and  weights  is  optional 


UNCLASSIF  IED 


Security  Classification 

k. 

__ 

SECTION  I I I 


HYDRAULIC  OPERATIONAL  SEVERITY  ASSESSMENT 


PROJECT  STAFF 

S.  K.  R.  Iyengar,  Co-Project  Manager 
R.  L.  Decker,  Co-Project  Manager 
R.  F.  Sharp,  Project  Engineer 
R.  L.  Brown,  Project  Engineer 
M.  T.  Yokley,  Project  Engineer 


i 


FOREWORD 

This  section  presents  the  results  of  the  pro  ject  on  'w  , • , system  operational 
severity  assessment . The  principal  object iv es  of  this  year'  effort  were  t<  examine 
the  changes  in  component  performance  parameters  due  to  their  .)■  jradation  and  relate 
them  to  cumulative  data  of  the  type  acquired  by  the  Statistical  Analou  Monitor  (ST  AM). 
Effort  was  also  to  be  expended  on  collecting  ST  AM  data  on  selected  machines  to  exa- 
mine the  influence  of  operational  duty  cycles  on  STAM  profiles.  This  report  presents 
a summary  of  data  so  collected. 
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CHAPTER  I 


INTRODUCTION 


The  assessment  of  performance  degradation  is  of  vital  importance  to  both  equip- 
ment manufacturers  and  users  of  hydraidic  equipment,  since  the  drop  in  performance 
capabilities  of  a machine  can  almost  always  be  traced  to  degradation  of  one  or  more 
components.  Significant  improvement  in  overall  system  reliability  can  be  achieved 
by  monitoring  the  performance  of  individual  components  and  repairing  or  replacing 
them  when  needed.  Even  though  the  degradation  of  a component  can  usually  be  ascer- 
tained by  disassembly  and  inspection,  this  process  is  generally  time-consuming,  and 
often  impracticable  for  field-maintenance  units.  Consequently,  the  development  of 
non-intrusi ve  diagnostics,  wherein  the  state  of  a component  is  assessed  without  re- 
moving components  from  a system  and  if  possible,  even  disturbing  normal  operation  of 
the  equipment,  hold  much  promise  for  improving  overall  system  reliability. 

Non-intrusive  diagnostics  depends  on  the  measurement  of  physical  quantities  like 
pressure,  temperature,  flow  rate,  etc.,  by  permanently  installing  the  requisite  trans- 
ducers in  the  system.  Hitherto,  such  instrumentation  had  been  accomplished  only  on 
experimental  vehicles  and  needed  the  use  of  elaborate  recording  or  transmitting  equip- 
ment. The  development  of  the  Statistical  Analog  Monitor  (S  1 AM)  has  obviated  the 
need  for  much  real-time  data  acquisition.  Details  of  construction  and  operation  of 
ST  AM  have  been  presented  in  earlier  reports  [1  ,2,  f]  as  well  as  a number  of  papers 
[4,5,6]. 
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Even  though  the  collection  of  cumulative  data  is  a challenging  task,  such  data,  by 
themselves,  reveal  little  .oformation  about  the  state  of  a specific  component  on  a mach'ne 
The  development  of  data  interpretation  algorithms  which  can  furnish  qualitative  or 
quantitative  information  about  the  state  of  the  system  or  any  component  therein  needs, 
as  a prerequisite,  an  analysis  of  the  system.  Examples  of  such  analysis  have  been 
presented  in  earlier  reports  [7,8].  Results  of  some  STAM  data  interpretation  are 
given  in  Ref.  [2].  The  effort  reported  in  this  document  is  a continuation  of  the  earlier 
work  and  focuses  on  individual  component  degradation. 

Chapter  I I examines  the  behavior  of  a system  in  which  components  are  exchanged 
for  those  of  larger  or  smaller  size.  It  is  shown  that  both  power  output  and  overall 
efficiency  may  be  significantly  altered  by  such  component  changes.  Chapter  I I I 
reports  on  component  performance  changes  due  to  degradation.  A relief  valve  was 
chosen  for  experimentation,  since  such  valves  are  critical  components  on  many  mobile 
hydraulic  systeims. 
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CHAPTER  I I 


EFFECT  OF  COMPONENT  CHANGES  ON  OPERATIONAL  SEVERITY 
INTRODUCTION 

Fluid  power  equipment,  like  any  other  piece  of  machinery,  is  designed  to  perform 
in  a prescribed  manner  under  the  range  of  specified  operating  conditions.  Even  though 
the  end-item  company  is  generally  responsible  for  the  call-out  of  individual  components, 
the  equipment  buyer  should  exercise  some  control  over  component  selection  [9].  In 
order  to  exercise  such  control  effectively  and  efficiently,  and  also  to  be  able  to  under- 
take or  guide  "commercial  modification"  to  suit  specific  operational  requirements,  it  is 
desirable  that  the  effect  of  changing  circuit  components  be  evaluated.  In  the  case  of 
simple  hydraulic  circuits,  the  effect  of  sizing  changes  can  be  estimated  by  quick  hand 
calculations.  Thus,  in  a circuit  comprised  of  a fixed  displacement  pump,  a cylinder, 
and  a four-way,  two-position  directional  control  valve,  enlarging  the  pump  displace- 
ment or  reducing  the  cylinder  will  increase  the  piston  velocity  proportionately.  How- 
ever, if  the  circuit  becomes  complex  due  to  aggregation  of  a number  of  actuators  and 
use  of  flow  modulation  through  the  use  of  various  types  of  flow  and  pressure  controls, 
hand  calculations  become  cumbersome  and  time-consuming,  if  not  altogether  impossible; 
under  such  circumstances,  it  is  advantageous  to  develop  appropriate  mathematical 
models  of  individual  components  and  assemble  them  to  obtain  a mathematical  model  for 
the  entire  system  under  investigation.  The  mathematical  model  for  the  system  can, 
then,  be  used  to  study  the  effect  of  changing  components  on  the  operational  severity 


of  the  system  under  prescribed  operating  conditions.  Performance  appraisal  using 
mathematical  models  of  components  can  be  very  efficiently  and  economically  done  with 
the  help  of  digital  computers.  An  example  of  such  analysis  is  presented  below, 
primarily  to  indicate  the  scope  of  effort  involved  and  the  type  of  appraisal  information 
generated  by  such  analysis.  Numerical  values  for  various  design  parameters  (e.g.  , 
pump  and  cylinder  sizes,  duty  cycles,  etc.)  were  chosen  so  as  to  illustrate  the  method- 
ology and  do  not  reflect  any  actual  operating  conditions. 


EXAMPLE  SYSTEM  DESCRIPTION 

The  open  center  system  consisting  of  a fixed  displacement  pump,  a directional 
control  valve,  and  two  hydraulic  cylinders  in  parallel  , as  shown  in  Fig.  2-1  , was 
chosen  as  the  test  circuit.  This  type  of  system  is  not  only  common  to  many  mobile 
hydraulic  systems  but  also  the  basic  subsystem  in  many  large  and  complex  hydraulic 
machinery.  It  also  illustrates  the  kind  of  trade-offs  which  system  designers  and 
equipment  users  need  to  consider  in  selecting  components. 


Fig.  2-1.  Circuit  Schematic  of  Open  Center  System 
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The  dotted  lines  in  F ig.  2-1  demarcate  subsystems;  each  of  which  in  turn  is  com- 
prised of  subsystems.  Details  of  modeling  the  subsystems  and  individual  components 
are  discussed  only  briefly,  since  they  have  been  presented  elsewhere  110,11]. 

The  positive  displacement  pump  is  assumed  to  operated  at  a constant  speed.  Its 
flow  rate  is  assumed  to  drop  linearly,  so  that  its  mathematical  model  consists  of  an 
equation  relating  the  flow  rate  to  the  delivery  pressure.  The  filter  and  heat  exchanger 
though  shown  in  Fig.  2-1  do  not  enter  into  the  mathematical  model  for  Subsystem  1 — 
an  assumption  which  can  be  lifted  if  the  f low-pressure  characteristics  of  the  component 
are  known . 

The  directional  control  valve  comprises  Subsystem  2 and  is  shown  to  manually 
operated  and  infinitely  positioning.  It  could  just  as  well  be  pilot  operated  or  electri- 
cally actuated.  Anti-cavitation  checks  and  cross-over  relief  valves,  which  are  often 
present  in  open  center  systems,  have  been  excluded  in  the  modeling.  The  valve  is 
described  mathematically  by  a set  of  graphs  which  show  how  the  various  metering 
orifices  change  with  spool  displacement,  Fig.  2-2.  Such  graphs  can  be  developed 
from  test  data  generated  in  metering  characteristics  tests  [12]. 

t 

The  actuators  comprise  Subsystem  3 and  are  mechanically  connected  so  that  they 
share  the  external  load  equally  and  move  in  unison.  It  is  also  assumed  that  the 
actuators  do  not  leak  and  do  not  cause  any  significant  drag.  Pressure  drops  in  the 
various  lines  are  also  considered  negligible. 
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In  order  to  keep  the  mathematical  models  as  simple  as  possible,  all  dynamic  con- 
siderations have  been  excluded.  Thus,  the  effects  of  fluid  compressibility  and  inertia 
of  moving  parts  in  the  cylinders  and  elsewhere  have  been  ignored.  The  end  result  is 
that  the  mathematical  models  of  all  components  are  in  the  form  of  algebraic  equations 
or  in  graphic  form,  in  the  case  of  the  directional  control  valve.  The  system  model 
obtained  by  assembling  the  subsystem  model  is  called  a static  model  to  distinguish  it 
from  one  which  portrays  dynamic  behavior.  It  should  be  noted  that,  even  though  it 
is  called  a "static"  model,  it  still  depicts  the  motion  of  the  actuators  and  the  changes 
in  the  various  flows  and  pressures  when  a time-history  of  inputs  is  given. 

The  mathematical  model  for  any  system  can  be  thought  of  as  a black  box  which 
accepts  certain  input  information  and  generates  the  desired  output  information.  In  the 
case  of  the  open  center  system  under  consideration,  the  following  quantities  are  consid- 
ered as  inputs:  pump  theoretical  flow,  valve  spool  position,  external  load,  and  tank 

pressure.  Of  those,  pump  flow  and  tank  pressure  are  treated  as  constants,  while  tie 
other  two  vary  with  time  in  accordance  with  the  prescribed  operational  cycle. 

Using  the  nomenclature  of  Table  2-1  , the  mathematical  models  for  the  subsystems 
can  be  shown  to  be  as  follows: 


Pump  Subsystem 


Q 


th 


- k P 
P s 


(2-1) 
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TABLL  2-1.  Nomenclature  for  Open  Center  byt>tem  Analysis 


SYMBOL 

COMPUTER 

NAME 

DESCRIPTION 

VALUE 

U S.  UNITS 

S.  1.  UNITS 

Ahc 

CYLH 

Cylinder  Effective  Area,  Head  End 

@ sq.  ins. 

mm^ 

Arc 

CYLR 

Cylinder  Effective  Area,  Rod  End 

@ sq  ins. 

mm^ 

a1 

A1 

Open-Center  Valve  Metering  Orifice 
Area 

sq.  in. 

" 

a2 

A2 

Open  Center  Valve  Metering  Orific*. 
Area 

sq.  in. 

■■ 

a3 

A3 

Open-Center  Valve  Metering  Orifice 
Area 

sq.  in. 

■■ 

a4 

A4 

Open-Center  Valve  Metering  Orifice 
Area 

sq.  in. 

" 

a5 

A5 

Open-Center  Valve  Metering  Orifice 
Area 

sq.  in. 

- 

HP.n 

HPIN 

Input  Horsepower 

HPout 

HPOUT 

Output  Horsepower 

^in 

Average  Input  Horsepower 

^out 

Average  Output  Horsepower 

k 

_ , 

XK1 

Pump  Slip  Coefficient 

Pa 

PA 

Port  A Pressure 

psi 

MPa 

Pb 

PB 

Port  B Pressure 

psi 

MPa 

Ps 

Supply  Pressure 

psi 

MPa 

Pt 

PTANK 

Tank  Port  Pressure 

50  psi 

0.345  Pa 

Qhc 

QHC 

Flow  to  Head  End  of  Cylinders 
(inflow  positive) 

gpm 

litres/minute 

Ore 

QRC 

Flow  from  Rod  End  of  Cylinders 
(outflow  positive! 

gpm 

litres/minute 

Qs 

FLOWIN 

Pump  Flow 

gpm 

litres/minute 

°t 

QTA  NK 

Flow  Through  a^ 

gpm 

litres/minute 

°th 

- TFLOW 

Theoretical  Pump  Flow 

gpm 

litres/minute 

V 

VELOC 

Actuator  Velocity  Extension  Positive 

in/sec 

mm/s 

w* 

ELOAD 

Load  on  Cylinders  (positive  if  resisting 
extension) 

Ibsf 

N 

x * 

SPOOL 

Spool  Displacement 

ins. 

mm 

P 

RHO 

Fluid  Density 

Ibsf  sec^/in* 

kg/t 

CYLDIA 

Cylinder  Diameter 

@ ins. 

mm 

ROD2 

Rod  Diameter 

2.5  ins. 

63.5  mm 

0 Parameter  (constant  for  a pwan  operational  cycle) 

* Venable 
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This  equation  indicates  that  at  any  time  the  actual  outflow  from  the  pump  is  equal 
to  the  theoretical  flow  less  the  slip  flow  which  is  proportional  to  the  supply  pressure. 

Valve  Subsystem 

Figure  2-2  presented  the  valve  metering  characteristics  graphical ly . Since  it 
is  rarely  practicable  to  describe  such  characteristics  in  the  form  of  equations,  a 
numerical  representation  is  used  in  the  analysis.  Such  representation  takes  the  form 
of  arrays  of  numbers  from  which  the  metering  areas  corresponding  tc  a given  spool 
displacement  can  be  read  off  either  directly  or  by  a process  of  interpolation.  The 
numerical  form  of  representation  is  particularly  well-adapted  to  computerization  [7,dj. 

Actuator  Subsystem 


Three  equations  are  needed  to  describe  the  actuator  operation.  The  first  is  a 
force  balance  equation  and  relates  the  external  load  to  the  pressures  on  either  side 
of  the  actuator;  the  other  two  equations  relate  the  actuator  velocity  to  the  inflow  and 
outflow  of  the  cylinder. 


PA,  - P,  A = W/2  (2-2) 

a he  b rc 


v 


Q /( 2 A 
he  he 


Q /( 2 A 
rc 


rc 


) 


) 


(2-3) 


(2-4) 
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It  may  be  noted  that  these  equations  are  algebraic  and  imply  that  the  actuator 
subsystem  has  been  described  by  a static  model. 


TOTAL  SYSTEM  REPRESENTATION  AND  SIMULATION 

Equations  (2-2)  through  (2-4),  in  conjunction  with  the  numerical  model  for  the 
directional  control  valve,  comprise  the  mathematical  model  for  the  entire  system.  It 
should  be  noted  that  the  equations  are  "coupled";  i .e.  , some  quantities  appear  in  more 
than  one  equation.  Consequently,  the  entire  set  of  equations  has  to  be  solved  simul- 
taneously for  each  set  of  inputs.  Since  a typical  operational  duty  cycle  will  involve 
continuously  changiny  spool  displacement  and  loads,  it  is  necessary  to  solve  the  set 
of  equations  representing  the  total  system  for  a number  of  sets  of  inputs.  By  doing 
so,  it  will  be  p"ssible  to  establish  the  output  quantities — namely,  actuator  velocity, 
cylinder  port  pressures,  flows,  and  supply  pressure  corresponding  to  each  set  of 
inputs.  Once  the  values  of  the  above-mentioned  outputs  are  established,  it  is  relatively 
easy  to  calculate  supplementary  quantities  which  may  be  of  interest.  Thus: 

Hydraulic  Power  Input  (HP.  ) = P Q HP 

in  s s 

6599 

Output  Power  (HP  .)  - VV . v HP 
out  6600 

Hydraulic  System  Efficiency  - LIP  ./HP. 

7 7 out  in 
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It  may  be  noted  that  all  of  the  above  quantities  change  from  instant  to  instant 
during  an  operational  duty  cycle,  and  average  values  need  to  be  calculated  by  account- 
ing for  the  time  period  at  each  level . 

Apart  from  calculating  instantaneous  values  of  pressures,  flows,  velocities, 
power,  etc.  , the  digital  computer  can  be  used  to  generate  cumulative  data,  such  as 
that  acquired  by  the  Statistical  Analog  Monitor  (STAM)  [10].  Simulation  of  STAM 
is  especially  useful,  since  it  can  be  used  for  the  rational  selection  of  transducer 
range  and  location,  gate  levels,  and  scale  factors  for  counters. 

PARAMETER  SENSITIVITY 

Once  a component  has  been  selected,  the  parameters  associated  with  it  are  fixed, 
no  matter  what  the  operational  duty  cycle.  Thus,  selection  of  a pump  fixes  the  theore- 
tical flow  rate,  while  selection  of  a cylinder  determines,  once  for  all,  the  ratio  of  head 
end  to  rod  end  flow  and  the  rod  velocity  for  a given  flow.  If  the  system  designer  or 
the  equipment  buyer  wishes  to  ascertain  how  a change  in  component  specifications  will 
affect  the  overall  performance,  one  method  is  to  actually  carry  out  the  change  on  the 
hardware  and  measure  output  quantities  under  both  circumstances.  The  expense  and 
time  involved  in  such  experimentation  is  not  inconsiderable.  On  the  other  hand, 
substantially  the  same  information  can  be  obtained  by  performing  what  is  known  as 
parameter  sensitivity  analysis.  One  method  of  doing  such  analysis,  and  often  the 
only  one  when  the  system  is  complex,  consists  of  modifying  the  mathematical  models 
for  the  components  to  reflect  the  parameter  changes  and  compare  the  results  of  simu- 
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tion.  Due  to  the  general  availability  of  digital  computers  and  the  downward  trend  in 
computational  costs,  the  total  cost  of  such  computer-aided  analysis  is  generally  orders 
of  magnitude  less  than  equivalent  experimental  work.  The  remainder  of  this  section 
will  present  the  results  of  varying  cylinder  bore  and  the  effect  of  incorporating  a relief 
valve  in  the  open  center  system. 

In  order  to  examine  the  effect  of  changing  a parameter,  it  is  necessary  to  keep 
all  other  parameters,  as  well  as  the  operational  duty  cycle,  invariant.  It  is  import- 
ant to  emphasize  that  any  trends  observed  by  changing  a parameter  in  the  above  manner 
are,  strictly  speaking,  valid  only  for  the  duty  cycle  used  in  the  simulation.  For  the 
other  duty  cycles,  the  entire  analysis  has  to  be  repeated.  Such  iterative  analysis  is, 
however,  easily  performed  on  the  digital  computer. 

Figure  2-3  presents  one  of  the  input  time-histories  used  in  the  simulation.  Only 
two  input  quantities  are  depicted,  since  the  other  two — pump  theoretical  flow  and  tank 
pressure — are  invariant.  Figures  2-4a,  b,  and  c depict  the  instantaneous  values  of 
the  input  end  output  power  as  well  as  their  average  values  for  three  cylinder  bore 
diameters.  The  actual  average  input  power  is  shown  in  addition  to  that  indicated  by 
STAM.  It  is  seen  that  the  STAM  estimates  are  close  to  the  actual  average. 

The  sensitivity  of  any  specified  output  quantity  to  the  parameter  under  study  is 
best  presented  graphically,  plotting  the  parameter  value  on  the  x-axis  and  the  output 
on  the  y-axis.  Figure  2-5  is  such  a graph  and  shows  the  variation  of  input  and  output 
power  and  efficiency  with  cylinder  bore  for  a specified  operational  duty  cycle. 
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Fig.  2-3.  Typical  Input  Time-Histories  used  for  Simulation 
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Fig.  2-4a.  Time-Histories  of  Input  and  Output  Power  for  a System 
with  4"  (102mm)  Cylinders 
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Fig.  2-4b.  Time-Histories  of  Input  and  Output  Power  for  a System 
with  5"  (127mm)  Cylinders 
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-4c.  Time-Histories  of  Input  and  Output  Power  for  a System 
with  5j"  (140mm)  Cylinders 
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Figure  2-6  shows  another  operational  cycle — one  in  which  the  valve  spool  is 
maintained  at  one  or  the  other  extreme  position — with  the  result  that  all  pump  flow 
is  directed  to  either  the  head  end  or  rod  end  of  the  cylinder.  Three  different  pump 
flow  rates,  3.89,  7.79,  and  77.9  gpm  (14.7,  29.4,  and  294  liters  per  minute,  respec- 
tively), and  three  cylinder  bores,  4",  5",  and  5?"  (102  mm,  127  mm,  and  140  mm, 
respectively),  were  used.  Figures  2—7  and  2—8  depict  the  output  power  and  effi- 
ciency time-histories  for  a pump  flow  rate  of  7.7  gpm  (29.4  Ipm);  similar  trajectories 
were  obtained  for  the  other  flow  rates. 
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Fig.  2-6.  Operational  Cycle  with  Spool  in  Extreme  Position 
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Fiq.  2-8.  Efficiency  Trajectories  for  Open  Center  System  with 
Three  Different  Cylinder  Sizes 
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For  this  particular  operational  duty  cycle,  the  power  and  efficiency  trajectories 


have  the  same  general  shape  as  the  load  trajectory.  This  is  not  surprising  in  view 
of  the  fact  that  the  valve  spool  was  moved  to  one  extremity  or  the  other;  and,  conse- 
quently, the  cylinder  velocity  was  dependent  solely  on  the  pump  flow  rate  and  the 
cylinder  bore.  Consequently,  for  a given  flow  rate  and  cylinder  size,  the  output 
power  would  increase  with  load.  Since  the  power  loss  in  the  valve  is  constant,  the 
efficiency  correlates  with  output  power.  It  is  also  of  interest  to  note  that  the  smaller 
cylinder  size  results  in  higher  output  power  and  efficiency.  This  can  be  ascribed  to 
higher  load  velocities  and,  therefore,  higher  output  power  for  the  same  losses. 
Figures  2-9  and  2-10  depict  the  sensitivity  of  cycle  efficiency  and  output  power  to 
both  cylinder  bore  and  pump  flow  rate.  Such  graphs  are  useful  in  appraising  the 


F'  i.  2-9.  Variation  of  Open  Center  System  Cycle  t fficiency  due 
to  Changes  in  Cylinder  Size  and  Pump  h low  Rate 
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Fig.  2-10.  Variation  of  Open  Center  System  Power  Output  due  to 
Changes  in  Cylinder  Size  and  Pump  Flow  Rate 
(See  Fig.  2-9  for  symbol  explanation.) 

effects  of  simultaneous  changes  in  two  variables.  It  should  be  emphasized  that  such 
graphs  are  valid  only  for  the  operational  duty  cycle  used  in  the  analysis;  and,  if  a 
machine  is  subjected  to  a number  of  distinct  duty  cycles,  such  an  analysis  may  be 
repeated  for  every  cycle. 


The  above  discussion  illustrates  the  methodology  for  systematically  analyzing 
the  effect  of  component  changes  on  operational  severity.  In  view  of  the  number  of 
variables  and  parameters  involved  in  the  mathematical  model  , the  results  expose 
only  some  of  the  tr  ade-offs  in  component  selection  and  system  operations.  Thus, 
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SPOOL  DISPLACEMENT 


for  example,  if  the  circuit  shown  in  Fig.  2-1  was  augmented  with  an  adjustable  circuit 
relief  valve,  the  mathematical  model  for  the  relief  valve  is  all  that  is  needed  to  appraise 


its  impact  on  the  system  behavior. 

Figure  2-11  shows  an  operational  duty  cycle  in  which  the  spool  displacement  has 
a large  enough  absolute  value,  such  that  all  the  pump  flow  is  directed  to  the  cylinder. 


Fig.  2-1  1 . Operational  Duty  Cycle  for  an  Open  Center  System 
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Fig.  2-13.  Time-Histories  for  a System  with  a Relief  Valve 

(7.38  Imp)  at  substantial  ly  the  cracking  pressure.  Consequently,  when  the  relief 
valve  opens,  no  flow  is  available  to  the  actuator  and  the  output  work  is  zero.  The 
net  result  on  the  operational  cycle  is  a drastic  reduction  in  outpower  as  well  as 
efficiency,  as  summarized  below: 
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TABLE  2-2. 


Comparison 

of  Average  1 nput 

and  Output  Power 

and  Efficiency 

System 

Average 
1 nput  Power 

Average 
Output  Power 

Cycle  Efficiency 

Without  relief 
valve 

0.54  HP 

0.48  HP 

89% 

With  rel  ief  val  ve 
set  at  1100  psi 
(7.59  MPa) 

0.65  HP 

0.24  HP 

38% 

A method  has  been  developed  to  examine  the  effect  of  component  changes  on  opera- 
tional severity  of  hydraulic  systems  used  on  mobile  hydraulic  equipment.  It  does  not 
require  extensive  testing  of  actual  hardware,  since  it  uses  mathematical  models  which 
can  be  developed  from  standard  performance  evaluation  tests.  By  consolidating  the 
mathematical  models  for  the  components  or  subsystems  comprising  a system,  it  is 
possible  to  simulate  its  behavior  under  any  specified  operating  conditions  and  for  any 
specified  operational  duty  cycle.  Since  it  is  generally  far  easier  to  change  or  insert 
a component  mathematical  model  than  to  implement  the  change  on  the  actual  hardware, 
it  is  possible  to  ascertain  the  impact  of  changing  component  parameters  even  on  complex 
systems  efficiently  and  quickly  by  simulation.  An  open  center  hydraulic  system  has 
been  analyzed  for  sensitivity  of  power  input  and  output  and  efficiency  to  changes  in 
actuator  size,  as  well  as  introduction  of  a relief  valve  in  the  pump  delivery  line. 

Similar  sensitivity  analysis  can  be  performed  on  any  system  for  which  component 
mathematical  models  are  available.  Such  models  should  include  in  them  the  specific 
parameter  (e.g.,  cylinder  bore,  relief  valve,  cracking  pressure,  etc.)  relative  to  which 
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sensitivity  analysis  is  to  be  performed.  Since  the  results  of  such  analysis  are, 
strictly  speaking,  valid  only  for  the  specific  operational  duty  cycle  used  in  the  simula- 
tion, it  is  advisable  in  the  case  of  multi-task  machines  to  iterate  such  analysis  for  the 
different  operational  duty  cycles  the  machine  is  likely  to  encounter . 


CHAPTER  I I I 


EXPERIMENTAL  ASSESSMENT  OF  OPERATIONAL  SEVERITY 
INTRODUCTION 

Since  relief  valves  are  among  the  most  critical  components  in  hydraulic  systems  and 
their  degradation  could  lead  to  poor  overall  performance  or  even  malfunctioning  of  the 
equipment,  it  was  decided  that  experimental  effort  would  be  focussed  on  them.  The  ob- 
jective was  to  develop  a methodology  for  correlating  performance  degradation  under  con- 
trolled operating  conditions  with  real  time  and  STAM  histograms. 

OPERATIONAL  SEVERITY  EFFECTS  ON  CYCLIC  TEST  RESULTS 

Figure  3-1  presents  the  circuit  schematic  for  the  test  system.  The  fixed  displace- 
ment pump  is  operated  at  a constant  speed  of  1750  RPM.  The  cycling  valve  is  a relief 
valve  set  at  a high  pressure  and  vented  through  a solenoid  directional  control  valve. 

The  solenoid  is  energized  by  a cycle  timer  so  as  to  control  the  "on"  time  and  "off"  time 
portions  of  the  cycle.  The  test  valve  has  a needle  bypass  which  is  used  to  control  the 
flow  going  through  the  test  valve.  When  the  bypass  is  completely  closed,  all  the  pump  tlow 
is  directed  through  the  test  valve.  A flowmeter  and  fluid  conditioning  equipment  (filters 
and  heat  exchangers)  make  up  the  remainder  of  the  test  circuit. 

A two-stage  relief  valve  with  adjustable  pressure  setting  was  used  as  the  test  valve. 
Cycling  was  performed  with  a "new"  valve  and  an  "old"  valve.  The  new  and  old  valves 
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Fig.  3-1.  Hydraulic  Circuit 
Schematic  for  Valve  Performance 
Assessment 


CYCLING  VALVE 


were  obtained  by  assembling,  respectively,  a good  and  a degraded  main  stage  in  the 


stage  was  performed  without  disturbing  either  the  pilot  stage  or  the  pressure  setting 


Consequently,  any  change  in  the  performance  can  be  attributed  to  the  main  stage 


degradation.  Visual  inspection  of  the  degraded  main  stage  revealed  significant 


wear  of  the  seating  cone.  Other  surfaces  were  macroscopically  unaffected.  Figure 


3-2  shows  the  seating  cone  and  the  region  of  wear 


Wear  Region 
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Figures  3-3  through  3-8  depict  the  pressure  wave  forms  obtained  for  the  new 


and  old  valves  at  two  different  flow  rates.  The  lower  flow  rate  was  obtained  by 


adjusting  the  bypass  valve  so  that  only  the  requisite  flow  passes  through  the  test  valve 


for  the  "on"  time  of  the  cycle.  Table  3-1  summarizes  the  salient  features  of  the  pres 


sure  wave  forms  for  the  old  and  new  valves 
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Pressure  Waveform  for  Old  Relief  Valve  with  Cycling 
at  Rated  Flow 
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F ig.  3-8.  Pressure  Waveform  for  Old  Relief  Valve  with  Cycling 
at  Low  Flow 


TABLE  3-1  . Summary  of  Cyclic  Performance  of  Relief  Valves 
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A scrutiny  of  the  pressure  wave  forms  and  the  table  leads  to  the  following 
conclusions: 

(1)  The  system  pressure  at  the  high  flow  rate  is  slightly  reduced  due  to  wear  of 
the  main  stage.  There  is  no  reduction  for  the  low  flow  rate.  This  is  not 
surprising  considering  that,  at  the  low  flow  rate,  the  main  stage  does  not 
open — i.e.,  all  the  flow  goes  through  the  pilot  stage  only.  Since  the  pilot 
stage  is  exactly  the  same  for  both  new  and  old  valves,  the  invariance  of  the 
pressure  at  low  flow  is  understandable. 

(2)  The  oressure  rise  rate  remains  substantially  the  same,  whether  the  valve  is 
new  or  old.  The  difference  between  the  rise  rates  for  low  and  high  flows 
can  be  attributed  in  part  to  the  ,^ad  valve  being  partly  open  for  the  low  flow 
cycle.  Since  the  load  valve  functions  as  a nonlinear  orifice,  the  flow  through 
it  bears  a complex  relationship  to  the  pressure;  and,  consequently,  the  test 
valve  is  not  subjected  to  a step  change  in  flow  like  that  encountered  at  high 
flow  cycl ing. 

(3)  The  wave  forms  for  low  flow  cycling  are  substantially  the  same  for'  both  the 
old  and  new  valves;  whereas,  for  high  flow  rates,  they  are  significantly 
different.  For  example,  the  new  valve  does  not  exhibit  the  pressure  over- 
shoot shown  by  the  old  valve. 

(4)  The  pressure  ripple  as  measured  peak  to  peaK  at  the  test  valve  is  significantly 
less  for  the  new  valve  than  the  old  valve.  Inspection  of  Figs.  3-3  and  3-6 
reveals  that  the  fundamental  and  first  harmonic  of  the  pressure  ripple  are 
unchanged  in  frequency.  The  fundamental  frequency  of  29.17  Hz  corresponds 
to  the  shaft  speed  (1750  rpm)  while  the  first  harmonic  (233  Hz)  corresponds 

to  the  number  of  pumping  elements  (gear  teeth,  in  this  case). 
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OPERATIONAL  SEVERITY  EFFECTS  ON  STATIC  AND  DYNAMIC  BEHAVIOR 


The  "old"  valve  used  in  the  above-mentioned  tests  was  artificial  in  that  the  test 
valve  contained  a degraded  main  stage  operated  in  conjunction  with  a "new"  pilot 
stage,  and  this  combination  is  not  likely  to  be  met  under  normal  field  operating  condi- 
tions. Consequently,  additional  tests  were  conducted  using  the  degraded  main  stage 
with  the  degraded  pilot  with  which  it  was  originally  equipped.  Since  the  pressure 
adjustment  mechanism  had  to  be  dismantled  to  reach  the  pilot,  it  was  impractical  to 
ensure  that  the  pilot  spring  compression  would  be  maintained  the  same,  with  the  "old" 
and  "new"  pilots.  Hence,  it  was  decided  to  conduct  tests  with  the  pressure  setting 
screw  adjusted  so  that  the  valve  maintained  the  same  pressure  at  the  rated  flow  of 
27  gpm  (102  liters/m),  with  both  the  old  and  new  pilots. 

Table  3-2  and  Fig.  3-9  present  the  static  characteristics  of  valves  equipped 
with  new  and  old  pilot  stages  used  in  conjunction  with  a degraded  main  stage,  as  also 
that  of  a new  valve,  to  serve  as  reference.  It  is  seen  that  a degraded  main  stage 
does  not  manifest  itself  unless  the  flow  rate  is  fairly  low — i.e.  , 12%  of  the  rated  flow 
or  less.  A degraded  valve  (i.e.  , one  with  old  pilot  and  main  stage)  is  seen  to  lose 
its  regulation  characteristics  over  most  of  the  flow  range  of  the  valve.  This  valve 
also  exhibited  erratic  behavior  in  that  repeated  pressure  measurements  at  the  same 
flow  resulted,  in  one  case,  in  pressures  as  far  apart  as  100  psi  (690  kPa).  A more 
serious  failing  of  the  valve  with  degraded  main  and  pilot  stages  was  its  poor  dynamic 
characteristics . The  response  time  for  step  changes  in  flow  was  so  large  that  the 
one-second  on/one-second  off  cycle  which  was  used  for  earlier  tests  (See  Figs.  3-3 
through  3-8.)  did  not  allow  the  pressure  to  stabilize  in  the  "on"  part  of  the  cycle. 
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TAQLt  3-2.  Flow-Pressure  Characteristics 
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hig.  3-9.  Static  Characteristics  of  "Old"  and  "New"  Relief  Valves 


Consequently,  the  "on"  time  had  to  be  lengthened  to  1 . 5 seconds  to  permit  steady  state 
operation  for  one-third  of  the  "on"  time.  Figures  3-10  through  3-11  present  the  wave 
forms  obtained  for  cyclic  operation  with  both  old  and  new  pilot  stages  in  conjunction 
with  the  old  main  stage.  It  can  be  seen  that  the  wave  form  for  low  flow  operation  is 
not  significantly  different  from  that  obtained  earlier  for  the  new  valve,  though  it  should 
be  noted  that  the  maximum  pressure  is  much  less.  The  wave  form  for  high  flow,  how- 
ever, exhibits  overshoots  of  the  order  of  100%,  and  steady  state  operation  is  not 
reached  for  at  least  0.5  seconds. 
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It  is  instructive  to  examine  how  cumulative  data,  such  as  that  acquired  by  the 
Statistical  Analog  Monitor  or  equivalent  devices,  portray  the  change  in  performance 
due  to  degraded  pilot  and  main  stages.  Tables  3-3  through  3-5  present  the  time  above 
pressure  and  number  of  crossings  for  eight  pressure  levels,  just  as  they  would  be 
obtained  from  the  Statistical  Analog  Monitor.  It  is  seen  that  the  cumulative  data  do 
not  differ  appreciably  when  the  main  stage  only  is  degraded  but  deviate  significantly 
when  the  pikot  stage  is  degraded.  Hence,  it  is  concluded  that,  for  the  design  of 
valve  tested,  degradation  of  the  pilot  stage  is  more  readily  evidenced  in  cumulative 
data  than  that  of  the  main  stage.  Since  in  actual  practice  degradation  of  a valve 
in  the  field  entails  changes  in  both  the  pilot  and  main  stages,  cumulative  data  can  be 
used  to  discriminate  between  new  and  degraded  valves,  provided  bench-mark  profiles 
have  been  established. 
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CHAPTER  IV 


FIELD  DATA  INTERPRETATIONS 


INTRODUCTION 

The  assessment  of  performance  degradation  of  hydraulic  system  components  is 
of  vital  importance  to  equipment  users.  Examples  of  such  degradation  are: 

(i)  a pump  losing  its  output  flow  for  a given  speed  and  pressure,  due  to 
increase  of  slip  flow, 

(ii)  a relief  valve  cracking  at  a lower  pressure  than  the  original  setting,  and 

(iii)  a hydraulic  cylinder  moving  slower  than  designed  for,  due  to  internal 
leakage . 


A major  cause  of  component  degradation  is  contaminant  wear.  It  is  usually 
possible  to  assess  the  extent  of  such  wear  by  disassembling  the  component.  How- 
ever, such  a procedure  is  not  always  practicable  or  desirable,  and  hence,  the  need 
for  non-intrusi ve  diagnostics  (i.e.,  methods  of  assessing  component  degradation 
without  disturbing  the  system).  In  principle  this  can  be  done  by  installing  trans- 
ducers at  selected  locations  in  the  hydraulic  system  and  monitoring  pertinent  vari- 
ables. The  task  of  selecting  the  quantities  to  be  monitored,  locating  the  sites  for 
transducer  installation  and  interpreting  the  data  requires  analysis  of  the  specific 
system  and  cannot  be  discussed  in  general  terms.  However,  it  is  possible  to 
deduce  conclusions  as  to  the  operational  severity  of  a specific  task  or  mix  of  tasks 
from  such  data.  Chapter  I I 1 presented  results  on  operational  severity  assess- 
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merit  of  a type  of  component  (i.e., 


relief  valves)  under  laboratory  conditions.  The 


remainder  of  this  chapter  will  present  interpretation  of  test  data  acquired  on  actual 
machines  doing  different  tasks  or  mix  of  tasks.  The  data  was  acquired  bv  Statistical 
Analog  Monitors  (STAM)  which  were  developed  by  the  Fluid  Power  Research  Center, 
Oklahoma  State  University,  for  storing  cumulative  data  from  electrical  analog  signals 
of  a variety  of  transducers  [1  , 2,  7].  Details  of  data  collection  can  be  found  in 
Section  IV  of  this  report . The  first  set  of  data  pertains  to  a backhoe  loader  system 
while  the  second  pertains  to  a wheeled  tractor  assigned  for  various  tasks. 


BACKHOE  LOADER  SYSTEM 

Pump  delivery  pressure,  swash  plate  angle  and  outlet  temperature  were  the 
parameters  measured  on  a closed  center  backhoe  loader  hydraulic  system.  STAM 
units  for  each  parameter  were  capable  of  measuring  time  above  level  for  eight  levels. 
Duration  of  data  acquisition  was  90  hours  for  the  first  set,  which  included  74  hours 
of  operation  as  a loader  and  16  hours  as  a backhoe.  For  the  other  two  data  sets, 
test  duration  was  100  hours,  with  50  hours  each  of  backhoe  and  loader  operation 
for  the  second,  and  40  hours  loader  and  60  hours  backhoe  operation  for  the  third. 
Tables  4-1a,  4-1b  and  4-1c  include  data  as  read  off  the  STAM  units  installed  on  the 
system.  For  purposes  of  comparison,  data  has  been  normalized  and  presented  in 
Table  4-2  and  Figs.  4-1  through  4-3.  Normalization  is  performed  by  dividing  the 
time  above  a level  by  the  total  test  duration.  Histograms  convert  the  cumulative 
data  into  range  data  and  depict  the  level  of  the  parameter  at  eight  or  less  discrete 
steps,  which  may  be  unequal.  It  should  be  noted  that  if  the  lowest  level  does  not 
have  a 100%  reading,  the  shortfall  has  to  be  assigned  io  a transducer  output  of  zero. 
Such  is  the  case  in  the  pressure  data  on  all  three  operational  cycles,  and  hence, 
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TABLE  4-1a.  Data  Sheet  - Operational  Cycle  # 1 


Cycle  Description: 

Parameter: 
Transducer  Location: 

Loader  Operation  74  hrs 
Backhoe  Operation  16  hrs 

Pressure 

Pump  Outlet 

Channel 

Level 

T i me 

(psi) 

(hrs) 

1 

158 

80.5 

2 

560 

36.9 

3 

952 

34.6 

4 

1 ,352 

28.9 

5 

1 ,752 

28.2 

6 

2,168 

23.8 

7 

2,552 

21  .4 

8 

2,944 

21  .4 

Parameter: 

Angle 

Transducer  Location: 

Pump  Swash  Plate 

Channel 

Level 

Time 

(degrees) 

(hrs) 

1 

1.13 

132.6 

2 

3.09 

132.3 

3 

5.0 

60.5 

4 

6.91 

57.2 

5 

8.87 

57. 1 

6 

10.78 

54.5 

7 

12.7 

45.5 

8 

14.65 

1 .5 
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TABLE  4-1  b. 

Data  Sheet  - Operational  Cycle  ft 2 

Cycle  Description:  Loader  Operation  50  hrs 

Backhoe  Operation  50  hrs 

Parameter:  Pressure 

Transducer  Location:  Pump  Outlet 

Channel 

Level 

Time 

(psi) 

(hrs) 

1 

117.7 

88.6 

2 

510.1 

30.4 

3 

894.6 

22.9 

4 

1 ,279.0 

18.6 

5 

1,663.0 

13.7 

6 

2,048.0 

10.2 

7 

2,432.0 

7.2 

8 

2,817.0 

3.0 

Parameter: 

Angle 

Transducer  Location: 

Pump  Swash  Plate 

Channel 

Level 

Time 

(psi) 

(hrs) 

1 

1.5 

100 

2 

3.3 

100 

3 

5.2 

100 

4 

7.1 

100 

5 

9.0 

100 

6 

10.9 

98.1 

7 

12.8 

75.1 

8 

14.7 

52.1 
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TABLE  4-1c.  Datasheet 


Operational  Cycle  #3 


Cycle  Description: 

Parameter: 
Transducer  Location: 

Loader  Operation  40  hrs 
Backhoe  Operation  60  hrs 

Pressure 

Pump  Outlet 

Channel 

Level 

T ime 

(psi) 

(hrs) 

1 

158.0 

61.8 

2 

560.0 

42.7 

3 

952.0 

40.6 

4 

1 ,352.0 

36.2 

5 

1 ,752.0 

32.8 

6 

2,168.0 

25.8 

7 

2,552.0 

18.6 

8 

2,944.0 

17.3 

Parameter: 

Angle 

Transducer  Location: 

Pump  Swash  Plate 

Channel 

Level 

Time 

(denrees) 

(hrs) 

1 

1.1 

100 

2 

3.1 

100 

3 

5.0 

100 

4 

6.9 

100 

5 

8.0 

100 

6 

10.8 

100 

7 

12.7 

96.2 

8 

14.7 

72.9 
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TABLE  4-2.  Cumulative  Data  for  100  Hour  Operation  of  a Backhoe  Loader 


Pressure 

Percent  of  Total  Time  Spent  at  Value  for  Operational  Cycle 

(bars) 

1 

2 

3 

0 

19.5 

1 1.4 

38.2 

1 1 .0 

43.6 

58.2 

19.1 

38.6 

2.3 

7.5 

2.3 

65.7 

5.7 

4.5 

6.2 

93.3 

0.7 

4.7 

1 .4 

121  .0 

4.4 

3.5 

7.0 

150.0 

2.4 

3.0 

7.4 

176 

0 

4.2 

1 .1 

203 

21.4 

3.0 

17.3 

Swash  Plate 

Percent  of  T otal 

Time  Spent  at  Value  for  Operational  Cycle 

Angle 

(degrees) 

1 

2 

3 

3 

55 

0 

0 

5 

2 

0 

0 

9 

2 

2 

0 

11 

7 

23 

4 

13 

33 

23 

23 

5 

1 

52 

73 
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Histograms  for  Pump  Outlet  Pressure  Fig.  4-2.  Histograms  for  Swash  Plate  Angle 

for  Three  Operational  Cycles  for  Three  Operational  Cycles 


w 


Fig.  4-3.  Average  Pressures  and 
Swash  Plate  Angle  for  Different 
Mixes  of  T asks 


t~ » — — — l 

100  50  0 

PERCENT  TIME  BACKHOE 

the  inclusion  of  zero  as  a level  in  the  pressure  data  in  Table  4-2  and  Fig.  4-1  . 

From  the  above  mentioned  table  and  figure,  it  is  seen  that  Operational  Cycle  #2 
shows  a substantially  lower  average  pressure  than  Cycles  #1  and  #2.  It  is  also 
seen  that  in  the  latter  two  cycles,  a substantial  fraction  of  time  is  spent  at  2940  psi 
(203  bars).  From  Fig.  4-2  it  is  seen  that  as  the  fraction  of  time  the  machine  per- 
forms as  a backhoe  increases  the  average  flow  demand  increases,  as  evidenced  by 
the  trend  of  the  swash  plate  angle. 

From  the  data  on  pump  outlet  pressure  and  swash  plate  angle,  the  change  in 
mean  values  of  these  parameters,  as  the  mix  of  tasks  changes,  can  be  plotted  as 
shown  in  Fig.  4-3.  If  the  pump  is  assumed  to  rotated  at  constant  speed,  the  pump 
delivery  can  be  considered  proportional  to  the  swash  plate  angle,  and  consequently, 
the  product  of  the  average  pressure  and  the  average  swash  plate  angle  is  a good 
approximation  to  the  average  output  hydraulic  power.  F igure  4-4  is  a plot  of  this 
power , normal  ized  for  a 50:50  task  mi  x (i . e . , the  power  is  taken  to  be  unity  when 
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the  machine  works  as  a loader  50%  of  the  time  and  a backhoe  for  the  remaining  50%). 

It  is  seen  that  operation  of  the  machine  as  a backhoe  is  more  severe  from  an  energy 
standpoint  than  usage  as  a loader. 

It  needs  to  be  emphasized  that  all  data  were  collected  during  controlled  tests 
and  the  conclusions  drawn  above  are  valid  for  operating  conditions  which  resemble 
the  test  conditions.  Deviation  from  the  latter  could  result  in  significantly  different 
test  data. 

TRACTOR  SYSTEM 

Tables  4-3  and  4-4  present  data  collected  by  STAM  on  a tractor  system  engaged 
in  two  different  tasks.  Pressures  were  measured  at  the  pump  outlet  as  well  as  an 
actuator.  The  data  illustrates  the  need  for  selecting  the  location  of  transducers 
carefully.  Thus,  in  both  cycles  the  fractions  of  total  time  spent  by  pump  pressure 
above  29.42  bars  is  significant,  whereas  for  cylinder  pressure  the  valve  is  much 
less.  Consequently,  pump  delivery  pressure  measurement  would  not  give  a true 
picture  of  the  severity  to  which  the  cylinder  is  being  subjected.  In  the  case  of 
Cycle  #2,  the  counts  for  the  pump  pressure  are  significantly  higher  than  for  the 
cylinder  pressure.  Pump  ripple  would  be  part  of  the  explanation.  It  should  be 
noted  that  the  STAM  has  a round-off  error  which  results  in  a very  small  number  being 
read  as  zero.  This  is  the  explanation  for  the  zero  counts  for  cylinder  pressure  in 
Cycle  #2  even  though  the  machine  has  spent  a finite  time  above  each  pressure  level. 
Another  inference  from  Table  4-4  is  that  the  time  above  a pressure  of  62.35  bars  is 
probably  due  to  a large  transient  which  exceeds  258  bars  (3740  psi).  Such  periodic 
transients  would  result  in  cumulative  time  at  all  intermediate  levels  but  would  show 
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TABLE  4-4.  Cylinder  Pressure  Cumulative  Data  for  Two  Different  Operational 
Duty  Cycles 
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Cycle  #2  Duration  - 65.7  hrs 


an  identical  count  for  the  number  of  crossings. 


Figures  4-5  and  4-6  present  histograms  of  system  temperature  for  four  different 
operating  conditions.  Since  the  tractor  was  performing  two  types  of  jobs  for  each  data 
set,  it  is  not  possible  to  use  this  data  to  develop  typical  temperature  profiles  for  each 
kind  of  task.  It  is,  however,  seen  that  Operational  Cycle  #2  not  only  results  in  the 
highest  temperature  attained  but  also  the  highest  mean  temperature.  Additional  data 
on  ambient  conditions  and  power  input  and  output  are  needed  before  an  assessment  of 
thermal  behavior  of  the  system  can  be  done. 


Fig.  4-4.  Normalized  Average  Power  for  Different 
Mixes  of  T asks 
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Fig.  4-5.  Histograms  for  System  Temperature 
for  Operational  Cycles  #1  and  #2 


Fig.  4-6.  Histograms  for  System  Temperature 
for  Operational  Cycles  #3  and  #4 
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CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

Hydraulic  system  components  are  subject  to  degradation,  primarily  as  a result 
of  contaminant  wear.  It  is  often  impractical  or  undesirable  to  dismantle  systems 
in  the  field  to  examine  the  state  of  system  components  and  take  remedial  action. 

Since  component  degradation  is  manifested  as  a change  in  operational  performance 
parameters,  a monitoring  of  the  latter,  in  principle,  permits  one  to  assess  the  state 
of  any  component.  The  continuous  monitoring  and  recording  of  system  parameters 
as  a machine  operates  in  the  field  is  rarely  feasible  due  to  technical  or  economic 
reasons.  To  alleviate  this  situation  the  Fluid  Power  Research  Center,  Oklahoma 
State  University,  developed  a Statistical  Analog  Monitor  (STAM)  to  record  cumulative 
data.  Since  such  data  cannot,  by  themselves,  be  used  to  generate  time-histories  of 
the  various  parameters  over  the  period  of  data  acquisition,  alternative  schemes  of 
data  acquisition  and  alternative  schemes  of  interpretation  had  to  be  developed. 

These  techniques  depend  on  mathematical  models  in  simulating  operational  severity 
has  been  explained  in  Chapter  I I . The  chapter  also  indicates  how  modeling  can 
be  used  to  predict  the  effect  of  changing  components  on  a given  system.  The  sub- 
sequent chapter  has  presented  experimental  work  which  reveals  the  effect  of  con- 
taminant wear  on  operational  performance  of  a component.  Typical  field  data 
collected  on  machines  operating  under  controlled  conditions  has  been  presented  in 
Chapter  IV. 

This  phase  of  the  reliability  project  has  demonstrated  the  importance  of 
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mathematical  analysis  of  a system  in  order  to  evolve  a rational  non-intrusive  diagnos- 
tic method  based  on  field  data.  The  interpretation  of  such  data  depends  crucially 
on  a thorough  analysis  of  the  system.  Even  the  installation  of  sensors  for  measuring 
physical  variables  and  STAM  threshold  settings  should  be  based  on  such  analysis. 


111-60 


REFERENCES 


1 . "On-Board  Hydraulic  System  Monitor  Study,"  prepared  for  U.S.  Army  MERDC 
by  the  Fluid  Power  Reseaneh  Center,  Oklahoma  State  University,  Stillwater, 
Oklahoma,  AD  779  468,  1973. 

2.  "On-Board  Hydraulic  System  Monitor  Study,"  prepared  for  U.S.  Army  MERDC 
by  the  Fluid  Power  Research  Center,  Oklahoma  State  University,  Stillwater, 
Oklahoma,  FPRC  - 4M5,  1974. 

3.  "Hydraulic  System  Reliability  Program,"  prepared  for  U.S.  Army  MERADCOM 
by  the  Fluid  Power  Research  Center,  Oklahoma  State  University,  Stillwater, 
Oklahoma,  AD  A032642,  1976. 

4.  Branstetter,  R.  L.,  "The  Answer  to  Low-Cost  Data  Acquisition,"  Paper  P73-RQ-2, 
Basic  Fluid  Power  Research  Proqram,  Fluid  Power  Research  Center,  Oklahoma 
State  University , Stillwater,  Oklahoma,  1973. 

5.  Branstetter,  R.  L.  and  G.  A.  Roberts,  "The  Statistical  Analog  Monitor — An 
Update,"  Paper  P74-14,  Basic  Fluid  Power  Research  Program,  Fluid  Power 
Research  Center , Oklahoma  State  University,  Stillwater,  Oklahoma,  1974. 

6.  Roberts,  G.  A.  and  G.  E.  Maroney,  "A  Statistical  Monitoring  Concept  for 
Accurate  Prediction  of  Component  Reliability,"  Proc.  4th  International  Fluid 
Power  Symposium,  sponsored  by  British  Hydromechanics  Research  Associa- 
tion, Sheffield,  U.K.,  1975. 

7.  "Hydraulic  System  Controls  Study  (V),"  Annual  Report,  prepared  for  U.S. 

Army  MERADCOM  by  the  Fluid  Power  Research  Center,  Oklahoma  State 
University,  Stillwater,  Oklahoma,  AD  779  628,  1973. 

8.  "Hydraulic  System  Controls  Study,"  prepared  for  U.S.  Army  MERADCOM  by 
the  Fluid  Power  Research  Center,  Oklahoma  State  University,  Stillwater, 
Oklahoma,  FPRC  - 4M1  , 1974. 

9.  Hopler,  P.  D.,  "Hydraulic  Systems  and  Components  for  U.S.  Army  Mobile 
Equipment,"  29th  National  Conference  on  Fluid  Power,  Cleveland,  Ohio,  1973. 

10.  Iyengar,  S.  K.  R.  , "A  Prognostic  Simulation  Technique  for  Fluid  Power 
Systems,"  Trans.  SAE  Paper  750821,  1975. 

11.  I yengar  , S . K . R . and  D.  G.  Miller,  "Static  Analysi s of  Di rectional  Control 
Valves  Using  Analytical  and  Test  Data,"  Paper  P73-SP3,  Annual  Report  #7, 

Basic  Fluid  Power  Research  Program,  Oklahoma  State  University,  Stillwater, 
Oklahoma,  1973. 

12.  "Hydraulic  System  Controls  Study  (V),  Section  I I ,"  Annual  Report,  prepared 
for  U.S.  Army  MERADCOM  by  the  Fluid  Power  Research  Center,  Oklahoma 

F ate  University,  Stillwater,  Oklahoma,  AD  758  876. 


1-61 


V 


h 


UNCLASSIFIED 

Security  Classification 

DOCUMENT  CONTROL  DATA  ■ R&D  — — — — — 

fSecufKy  classllication  ol  title  body  ol  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  ia  classified) 

1 ORIGINATING  ACTIVITY  (Corporate  author)  |2a  REPORT  SECURITY  CLASSIFICATION 

Fluid  Power  Research  Center  Unclassified 

Division  of  Engineering,  Technology  & Architecture  **>  group 

Oklahoma  State  University.  Stillwater.  OK  74074 j 

3 REPORT  TITLE 

MERADCOM/OSU  HYDRAULIC  SYSTEM  RELIABILITY  PROGRAM 
On-Board  Monitor  Study 

4 DESCRIPTIVE  NOTES  (Type  ol  report  and  inclusive  dates) 

Annual  Report 1 February  1976  - 31  January  1977 

5 AUTHOR(S)  (Laat  name.  (Irat  name.  Initial) 

Staff  of  the  Fluid  Power  Research  Center 


6 REPO  RT  DATE 

February  1977 

7 a TOTAL  NO  OF  PAGES  | 7 b NO  OF  REFS 

72  I 0 

8 a CONTRACT  OR  GRANT  NO 

DAAK02-75-C-01  37 

b PROJEC  T NO 

c 

d 

9a  ORIGINATOR'S  REPORT  NUMBER^) 

E PRC-7M3 

96  OTHER  REPORT  NOfS.)  (A  ny  other  numbers  that  may  be  assigned 
this  report) 



10  * VA  IL  ABILITY/LIMITATION  NOTICES 

Approved  for  Public  Release:  Distribution  Unlimited 


11  SUPPLEMENTARY  NOTES 

12  SPONSORING  MILITARY  ACTIVITY 

Mobility  Equipment  Research  & 

Development  Command 

13  ABSTRACT 

This  report  presents  the  data  obtained  from  the  On-Board  Monitor  activities. 

These  data  are  a result  of  endeavors  on  two  projects — On-Board  Monitor  and 
Operational  Severity.  STAM  units  were  installed  on  a variety  of  vehicles. 

As  many  as  five  parameters  per  vehicle  were  monitored  with  several  tests  repeated 
with  multiple  runs. 


DD  1473 

L 


UNCLASS  If  ILD  

Security  Classification 

A 


UNCLASSIFIED 


Security  Classification 


Statistical 
Monitor 
Hydraul  ics 
Duty  Cycle 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  name  and  address 
of  the  contractor,  subcontractor,  grantee,  Department  of  De- 
fense activity  or  other  organization  ( corporate  author)  issuing 
the  report. 

2 a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over- 
all  security  classification  of  the  report#  Indicate  whether 
"Restricted  Data"  is  included.  Marking  is  to  be  in  accord- 
ance with  appropriate  security  regulations. 

26.  GROUP:  Automatic  downgrading  is  specified  in  DoD  Di- 
rective 5200.10  and  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3 and  Group  4 as  author- 
ized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  in  all  rases  should  be  unclassified. 

If  a meaningful  title  cannot  be  selected  without  classifies 
tion,  show  title  classificat;on  in  all  capitals  In  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  .oter  the  type  of 
report,  e.g. , interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  name(s)  of  authorfs)  as  shown  on 

or  in  the  report.  Entei  last  name,  first  name,  middle  initial. 

If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  «c.ithor  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report  as  day, 

month,  year;  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication.  « 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  i. e. , enter  the 
number  of  pages  containing  information. 

76.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 

references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

86,  8c,  & 8 d.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enler  the  offi- 

cial  report  number  by  v.’hicn  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

96.  OTHER  REPORT  NUMBFR(S):  If  the  report  has  been 
assigned  any  other  repon  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  number(s). 

10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  any  lim- 
itations on  further  dissemination  of  the  report,  other  than  those 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC.  " 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 


(4)  "U.  S.  military  agencies  may  obtain  copies  of  this 

report  directly  from  DDC  Other  qualified  users 
shall  request  through 


(5)  "All  distribution  of  this  report  is  controlled.  QuU- 
lfied  DDC  users  shall  request  through 


If  the  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi- 
cate this  fact  and  enter  the  price,  if  known. 

11.  SUPPLEMENTARY  NOTES:  Use  for  additional  explana- 
tory notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay" 
mg  for)  the  research  and  development.  Include  address. 

13.  ABSTRACT  Enter  an  abstract  giving  a brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re- 
port. If  additional  space  is  required,  a continuation  sheet  shall 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in- 
formation in  the  paragraph,  represented  as  rrSj.  (SJ  fC).  or  (U)  | 

There  is  no  limitation  on  the  length  of  the  abstract  How-  i 
ever,  the  suggested  length  is  from  150  to  225  words. 

14  KEY  WORDS:  Key  words  are  technically  meaningful  terms  I 
or  short  phrases  that  characterize  a eport  and  may  be  used  as  | 
index  entries  for  cataloging  the  repoi  Key  words  must  be 
selected  so  that  no  security  classification  is  required  Identi- 
fiers, such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con- 
text. The  assignment  of  links,  rules,  and  weights  is  optional. 


UNCLASSIFIED 


Security  Classification 


I 

1 


SECTION  IV 


ON-BOARD  MONITOR  STUDY 


i 


Project  Staff 

R.  L.  Decker,  Program  Manager 
M.  T.  Yokley,  Project  Engineer 
T.  Herron,  Project  Assistant 


FOREWORD 

This  section  presents  the  data  obtained  from  the  On-Board  Monitor  activities. 
These  data  are  a result  of  endeavors  on  two  projects — On-Board  Monitor  and  Opera- 
tional Severity.  As  many  as  five  Statistical  Analog  Monitor  (STAM)  units  were  in- 
stalled on  a single  vehicle  measuring  a variety  of  parameters.  Some  of  the  tests 
were  repeated  on  the  same  vehicle. 
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CHAPTER  I 


INTRODUCTION 

The  Statistical  Analog  Monitor  is  a device  developed  to  obtain  long-term  informa- 
tion regarding  the  operation  of  hydraulic  systems  in  actual  field  environments.  It  is 
a small,  low-cost  device  ideally  suited  for  installation  on  mobile  equipment.  The  de- 
vice operates  unattended  while  gathering  information,  typically  in  100-hour  tests.  The 
theory  of  operation,  design,  and  data  reduction  techniques  are  discussed  in  previous 
reports . 

A cooperative  program  was  initiated  in  1974  involving  15  major  companies  which 
are  producers  of  mobile  equipment  similar  to  the  type  utilized  by  the  U.S.  Army.  This 
was  directed  primarily  toward  obtaining  temperature  and  pressure  characteristics  of 
hydraulic  systems  on  construction  equipment.  A report  of  the  results  of  this  program 
through  January,  1975,  is  contained  in  Ref.  [1].  At  that  time  several  units  were  still 
in  the  field.  Rather  than  recall  these  units  and  possibly  invalidate  tests,  the  STAM 
monitors  were  allowed  to  remain  in  the  field. 

This  report  documents  the  results  obtained  from  activities  utilizing  the  STAM 
during  this  year.  Two  projects  utilized  the  STAM  this  year.  The  On-Board  Monitor 
Project  is  a continuation  of  the  cooperative  program  previously  initiated.  The  Opera- 
tional Severity  Project  utilized  tfie  STAM  to  collect  data  to  be  used  for  assessing  the 
operational  severity.  The  interpretation  of  these  results  is  contained  in  Section  I I 
of  this  report . 
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CHAPTER  I I 


DATA  ACQUISITION 

The  scope  of  the  programs  for  this  year  required  that  two  activities  be  simul- 
taneously pursued.  The  STAM  units  in  the  field  at  the  time  of  the  completion  of  the 
cooperative  program  were  to  be  allowed  to  complete  the  tests.  These  units  were, 
then,  to  be  interrogated  and  reconditioned.  In  addition,  sponsorship  for  a program 
involving  a concentrated  acquisition  effort  required  for  operational  severity  assess- 
ment was  to  be  obtained. 

Units  which  were  returned  from  companies  participating  in  the  cooperative 
program  were  interrogated  and  recalibrated.  Many  of  these  companies  had  been 
active  through  the  entire  program.  In  order  to  broaden  the  data  base,  units  were 
returned  to  selected  companies.  These  companies  had  demonstrated  an  aggressive 
role  in  the  utilization  of  the  STAM  monitor  for  the  acquisition  of  information.  During 
the  year,  several  of  these  companies  incorporated  the  STAM  into  an  integral  part  of 
their  own  programs. 

Sponsorship  for  acquisition  efforts  was  solicited  and  obtained.  An  extensive 
program  to  obtain  a larqe  amount  of  information  about  specific  machines  was  initiated. 
As  many  as  five  STAM  units  were  installed  on  a single  vehicle.  Tests  were  run  with 
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of  similar  model.  Because  of  the  commonalities  existing  between  the  hydraulic  systems 
both  industrial  equipment  and  agricultural  tractors  were  allowed. 


During  the  year,  more  than  45  units  were  shipped.  Companies  participating  included 
the  fol  lowing: 

Allis-Chalmers 

Massey-Ferguson 

Bucyrus-Erie 

International  Harvester 

J . I . Case 

Eaton 

In  addition,  units  were  supplied  to  various  areas  within  MERADCOM.  The  reliability 
of  the  STAM  was  excellent  during  this  program  with  only  one  electronic  failure  occur- 
ring. Previous  reports  documented  logistical  and  technical  problems  associated  with 
the  STAM  system.  The  two  important  ones  are  damaged  cables  and  tranducer  as  a re- 
sult of  abuse  and  electromagnetic  interference  (EMI)  susceptibility.  During  this  year, 
the  first  of  these  has  not  been  prevalent  while  the  latter  was  only  a problem  in  isolated 
cases.  These  cases  have  been  readily  identified  and  only  exist  on  certain  vehicles. 
Experience  has  shown  that  vehicles  utilizing  solenoids  are  a particular  problem. 

Appendix  A catalogs  selected  results  from  these  activities.  Some  of  the  data 
have  been  marked  with  a vehicle  numerical  designation.  In  these  cases  an  alphabetic 
suffix  has  been  added  to  distinguish  individual  tests.  Section  I I of  this  report  contains 
interpretations  of  these  results.  A few  additional  comments  are  in  order  regarding 
the  data  obtained  from  the  STAM.  The  following  paragraphs  discuss  several  of  these. 
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The  STAM  monitor  is  an  analog  process.  The  accuracy  of  the  STAM  is  a func- 
tion of  several  factors  including  transducer  accuracy,  internal  scaling,  changes  in 
internal  scaling  caused  by  time,  changes  in  internal  scaling  caused  by  temperature, 
and  read-out  accuracy.  An  overall  accuracy  of  one  percent  of  full  scale  for  the  entire 
system  (including  accuracy  of  all  intermediate  processes  as  a result  of  all  factors) 
would  be  considered  excellent  with  a five  percent  accuracy  considered  good. 

The  evaluation  of  the  data  must  include  considering  the  uncertainties  associated 
with  the  data.  The  standard  scaling  for  the  STAM  is  100  hours  and  75,000  counts 
full  scale.  Assuming  a one  percent  accuracy,  uncertainties  of  ± 1 hour  and  1 750 
counts  exist  in  the  data. 

An  often  overlooked  aspect  in  evaluating  the  data  is  considering  the  physical 
location  of  the  transducers  in  the  system.  A pressure  transducer  mounted  directly 
on  the  outlet  of  a pump  will  result  in  the  pump  ripple  affecting  the  data — both  time 
above  level  and  counts.  If  the  transducer  has  been  installed  at  the  end  of  the  line 
which  had  the  net  effect  of  attenuating  this  ripple,  the  data  could  possibly  be  signifi- 
cantly different.  The  "correct"  location  of  the  transducer  will  be  dependent  on  the 
needs  of  a particular  application.  Test  setups  which  count  pump  ripple  may  only 
provide  little  information  about  the  actual  work  of  the  pump. 

The  evaluation  of  the  data  presented  in  Appendix  A should  be  performed  with 
an  awareness  of  these  considerations . The  data  have  been  rounded  to  the  nearest 
hour  and  to  the  nearest  1 ,000  counts.  The  uncertainties  associated  with  the  data 
result  in  the  possibility  of  zero  time  with  a non-zero  number  of  counts  at  a given 
level . 
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CHAPTER  I I I 

SUMMARY  AND  CONCLUSIONS 

During  this  year  technical  and  logistical  support  was  provided  in  order  to  receive 
late  returns  from  the  cooperative  program.  Several  units  initially  committed  during 
the  program  were  returned  to  the  field  in  order  to  increase  data  bases.  Units  were 
supplied  to  satisfy  the  needs  of  several  areas  within  MERADCOM. 

During  this  year  support  for  the  program  was  solicited  and  obtained.  Extensive 
investigations  were  performed  on  several  vehicles.  The  results  of  this  were  excellent. 
Many  of  the  previous  logistical  problems  were  solved  as  a result  of  the  format  of  this 
program.  Participants  demonstrated  an  active  and  aggressive  role  which  resulted  in 
a mutually  beneficial  program. 

The  ST  AM  has,  in  the  past,  demonstrated  its  importance  in  gathering  long-term 
information  about  hydraulic  systems.  The  results  from  this  year  reinforce  these  con- 
clusions. The  STAM  has  been  utilized  reliably  to  obtain  information  needed  for  the 
evaluation  of  the  operation  of  vehicles  in  field  environments.  STAM  has,  indeed,  been 
accepted  by  industry  as  demonstrated  by  the  financial  support  it  has  received  this  year. 
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APPENDIX  A 


DATA  SHEET 


DATE:  May,  1976 


UNIT  NO 


COMPANY:  D 


UNIT  TYPE:  Pressure 


APPLICATION:  /Type  of  Vehicle  and  ideation  of  Sensors) 


Ag  Tractor  - Steering  and  Clutch  Pressure 

Open  Center  with  Low  Pressure  Regulator 


CHANNEL 


LEVEL 

(PSI) 


(HRS.) 


(Hrs.)  Total  Operation  Time 


REMARKS:  Vehicle  No.  1-A 


I V-1 1 


COUNTS 


>210,000 


>210,000 


>210,000 


76,000 


40,000 


16.000 


1 


FHkCTLlAi,a  PAGE  MOT  FILMED 

blame 


DATA  SHEET 


DATE  May,  1976 UNIT  NO.  1094 

COMPANY:  D 

UNIT  TYPE:  Pressure 

APPLICATION:  i Type  of  Vehicle  and  Locution  of  Sensors) 

Ag  Tractor  — Hitch  Pressure 

Hitch  Automatically  Adjusting  for  Constant  Draft  Load 


CHANNEL 

— 

LEVEL 

TIME 

(HRS.) 

COUNTS 

I 

400 

2 

144,000 

‘ 

710 

2 

91.000 

3 

1010 

i 

75.000 

4 

1320 

0.5 

39.000 

5 

1630 

0.5 

34.000 

6 

1940 

0.5 

27.000 

7 

2250 

0.5 

27.000 

8 

2550 

0 

0 

( Hrs.)  Total  Operation  Time 


REMARKS 


Vehicle  No.  1-A 


76-2 


DATA  SHEET 


DATE:  May,  1976  UNIT  NO  '080 

COMPANY:  D 

UNIT  TYPE:  Temperature 

APPLICATION:  / Type  of  Vehicle  and  location  of  Sensors) 

Ag.  Tractor  — Reservoir  Temperature 


\ 

CHANNEL  | LEVEL 

H 

TIME 

(HRS.) 

COUNTS 

1 81 

76 

41,000 

: ![  101 

69 

99,000 

3 121 

1 

64 

145,000 

4 i 141 

47 

173,000 

5 162 

II 

22 

183.000 

6 !j  182 

1 

3,000 

7 Ij  202 

1 

0 

0 

8 jj  222 

0 

0 

(Hrs.)  Total  Operation  Time 


REMARKS.  Vehicle  No.  1-A 


76-3 


DATA  SHEET 


DATE:  May.  1976 UNIT  NO  . : 10X0 

COMPANY:  [) 

UNIT  TYPE:  Pressure 

APPLICATION:  (Type  of  Vehicle  and  Location  of  Sensors) 

Ag  Tractor  - Auxiliary  Hydraulic  System 
Valve  Cutlet  Pressure 
Various  Operations 


1 

CHANNEL  LEVEL 

(PSD 

TIME 

(HRS.) 

COUNTS 

1 i 570 

23 

>210.000 

2 970 

21 

>210,000 

|| 

3 ;i  1380 

I 

14 

>210.000 

4 

1 

1780 

8 

>210.000 

5 

2180 

1 

145.000 

6 

2600 

0 

0 

7 

3000 

0 

0 

R 

3410 

0 

0 

90 (Hrs.)  Total  Operation  Time 


REMARKS:  Vehicle  No.  I -A  76-4 


I V — 1 7 


DATA  SHEET 


DATE  July.  1976  UNIT  NO  1080 

COMPANY  1) 

UNIT  TYPE:  Temperature 

APPLICATION:  /Type  of  Vehicle  a ml  Location  uf  Sensors) 


1 

CHANNEL  LEVEL 

TIME 
(HRS  ) 

COUNTS 

i 1 z 

99 

38.000 

2 101 

99 

93.000 

3 .21 

97 

1 1 9,000 

il 

4 HI 

77 

>210.000 

5 162 

34 

>210.000 

6 182 

1 

5 

86.000 

7 ij  202 

- I 

2 

10,000 

8 !!  222 

II 

2 

8.000 

100  (Hrs.)  Total  Operation  Time 


REMARKS  Vehicle  No.  1-B  76-5 


9 


DATA  SHEET 


DATE:  July,  1976 UNIT  NO  : 1080 

COMPANY:  D 

UNIT  TYPE.  Pressure 


APPLICATION:  (Type  of  Vehicle  and  Location  of  Sensors) 

Ag  Tractor  — Hitch  Cylinder  Pressure 


— 

CHANNEL 



— 

LEVEL 



TIME 

(HRS.) 

COUNTS 

1 470 

2 

*0 

2 |j  960 

1 

0 

| 

3 |j  1440 

1 

0 

4 

— 

1930 

1 

0 

s 

2400 

1 

0 

6 

2910 

1 

0 

7 

3400 

1 

0 

8 

3880 

1 

0 

100  (Hrs.)  Total  Operation  Time 


REMARKS:  Vehicle  No.  1-B  76-6 


IV-21 


COUNTS  (X10J) 


DATA  SHEET 


DATE:  July,  1976 UNIT  NO  : 1094 

COMPANY:  D 

UNIT  TYPE:  Pressure 

APPLICATION:  ( Type  of  Vehicle  and  Location  of  Sensors) 

Ag  Tractor  — Auxiliary  Pump  Pressure 


1 

CHANNEL  .j  LEVEL 

TIME 

(HRS.) 

COUNTS 

I 

50 

100 

>210,000 

__  2 1 

420 

21 

>210.000 

3 

790 

11 

208,000 

4 

1160 

7 

155.000 



S 1520 

6 

61.000 

6 1 880 

4 

74.000 

n 

/ 

2250 

3 

52.000 

8 

2620 

1 

5 1 .000 

100  (Hrs.)  Total  Operation  Time 


REMARKS  Vehicle  No.  1-B 


76-7 


DATA  SHEET 


DATE:  July,  1976 UNIT  NO. : _1094 

COMPANY:  D 

UNIT  TYPE.  Pressure 

APPLICATION:  /Type  of  Vehicle  and  Location  of  Sensors) 

Ag  Tractor  — Steering  and  Clutch  Pump 


100  (Hrs.)  Total  Operation  Time 


REMARKS:  Vehicle  No.  I-B  76-8 


IS 


DATA  SHEET 


DATE:  Sept..  1976 

COMPANY:  ^ 


UNIT  TYPE:  Temperature 


UNIT  NO 


1038 


APPLICATION:  (Type  of  Vehicle  and  Location  of  Sensors] 

Ag  Tractor 


[i — 

CHANNEL 

LEVEL 

TIME 

(HRS.) 

COUNTS 

1 

82 

66 

1.000 

104 

65 

22.000 

S 

125 

63 

78.000 

4 

146 

52 

91.000 

5 

167 

36 

196.000 

6 

189 

10 

7 

2 

1 2.000 

8 

11 

231 

i 

7.000 

66 (Hrv)  Total  Operation  Time 


REMARKS:  Vehicle  No.  1-C  76-9 


IV-I 


DATA  SHEET 


DATE:  Sept,  1976 UNIT  NO  1038 

COMPANY:  D 

UNIT  TYPE:  Pressure 

APPLICATION:  (Type  of  Vehicle  and  Location  of  Sensors 

Ag  Tractor  - Hitch  Cylinder  Pressure 


REMARKS  Vehicle  No.  1-C 


76-10 


DATA  SHEET 


DATE:  June,  1976 

COMPANY:  D 

UNIT  TYPE.  Pressure 

APPLICATION:  (Type  of  Vehicle  and  Location  of  Sensi'rs) 

Ag  Tractor  — Steering  and  Clutch  System  Pump  Pressure 

COUNTS 

0 

>210.000 
30.000 
3,000 
0 
0 
0 
0 


DATA  SHEET 


DATE:  June,  1976 UNIT  NO  1066 

COMPANY:  P 

UNIT  TYPE:  Pressure 

APPLICATION:  (Type  of  Vehicle  and  Location  of  Sensors) 

Ag  Tractor  - Draft  Control  Cylinder 


I 

CHANNEL  I LEVEL 

(PSD 

TIME 

(HRS.) 

COUNTS 

• U 

II  50 

100 

0 

2 i|  400 

65 

3 1 .000 

|| 

3 ||  760 

17 

5 1 ,000 

4 ||  1120 

2 

5,000 

5 ||  1470 

0 

0 

6 1830 

0 

0 

7 ' 2200 

- - - - II 

0 

0 

8 jl  2550 

0 

0 

I 00 (Hrs.)  Total  Operation  Time 


REMARKS:  Vehicle  No.  2-A 


76-13 


COUNTS  '(XIO  ) HOURS 


DATA  SHEET 


DATE:  June.  1976 UNIT  NO  1038 

COMPANY:  L) 

UNIT  TYPE:  Pressure  

APPLICATION.  (Type  of  Vehicle  and  Locution  of  Sensors) 

Ag  Tractor  - Hitch  Pressure 


CHANNEL 

LEVEL 

TIME 
(HRS. I 

COUNTS 

I 430 

3 

26.000 

1 

900 

2 

10.000 

3 

1380 

i 

4.000 

4 

1850 

0 

— 

0 

S 

2320 

0 

0 

6 

2790 

0 

0 

1 

3270 

0 

0 

8 

3740 

0 

0 

1 00 (Hrs.)  Total  Operation  Time 


REMARKS  Vehicle  No.  2-A 


76-15 


DATA  SHEET 


DATE:  Feb.,  1976  UNIT  NO.  1078 

COMPANY.  A 

UNIT  TYPE:  Temperature 

APPLICATION:  (Type  of  Vehicle  and  Location  of  Sensors) 

Ag  Tractor  — Auxiliary  liydraulic  System 
Load  Sensing  System 


CHANNEL 

LEVEL 

<°F) 

TIME 

(HRS.) 

COUNTS 

i 

78 

100 

6.000 

' 

99 

100 

13,000 

3 

120 

90 

45.000 

4 

140 

72 

209.000 

5 

160 

14 

183.000 

6 

180 

1 

2.000 

7 

202 

0 

0 

8 

223 

0 

0 

100 (Hrs.)  Total  Operation  Time 


RI  MARKS 


76- 1 6 


DATA  SHEET 


DATE  Aiu:.,  1976  UNIT  NO 

COMPANY:  C 

UNIT  TYPE  Pressure 

APPLICATION  1 rnc  of  Vehicle  jnJ  l.ocjlion  oj  Sen st  >n  / 

Loader  Backlme  Pump  Pressure 

Pressure  Compensated  System 


1 

CHANNEL  LEVEL 

(PSIl 

TIME 

(HRS.) 

COUNTS 

1 120 

89 

>412.000 

: 510 

30 

>412.000 

l1 

4 890 

23 

>412.000 

4 1 280 

18 

>412.000 

5 1660 

14 

>412,000 

<*  2050 

10 

>412.000 

2430 

7 

363.000 

s 2820 

3 

286.000 

1 00 (Hrv)  Total  Operation  Time 

REMARKS  Vehicle  No.  4-A  76-1 S 


I V -a5 


COUNTS  (X10J)  hours 


DATA  SHEET 


k 


DATE:  Aug.,  1976  UNIT  NO  1006 

COMPANY:  E 

UNIT  TYPE:  Rotational  Angle 


APPLICATION:  Type  of  Vehicle  ami  location  of  Sensors) 

Loader/Backhoe  Swash  Plate  Angle 

Pressure  Compensated  System 


REMARKS  Vehicle  No.  4-A 


76-20 


COUNTS  ( X 1 0 J ) hours 


DATA  SHEET 


DATE:  Sept.,  1976 UNIT  NO 

COMPANY:  

UNIT  TYPE:  Pressure 

APPLICATION:  tTvpe  of  Vehicle  and  Ixication  of  Sensors/ 

Industrial  Tractor 


CHANNEL 

— 

LEVEL 

(PSI) 

TIME 
(HRS  ) 

COUNTS 

i 

440 

1 

1 .000 

930 

0 

0 

3 

1420 

0 

0 

4 

1900 

0 

0 

5 

2390 

0 

0 

6 

2870 

0 

0 

7 

3360 

0 

0 

8 

3840 

0 

0 

(Hrs  ) Total  Operation  Time 


REMARKS  Vehicle  No.  3-A  76-23 


I V 


HOURS 


DATA  SHEET 

DATE:  Sept.,  1976 UNIT  NO  : '°36 

COMPANY:  E 

UNIT  TYPE:  Temperature 

APPLICATION:  t Type  of  Vehicle  and  Location  of  Sensors) 

Industrial  Tractor 


T 

CHANNEL 

LEVEL 

(°F) 

TIME 

(HRS.) 

COUNTS 

1 

81 

2 

8.000 

[I 

2 

102 

1 

1 0.000 

| 

3 1 

1 23 

i 

0.000 

4 

144 

0 

0 

5 

165 

0 

0 

6 

185 

0 

0 

*7 

206 

0 

0 

8 


227 


0 


0 


5 (Hrs.)  Total  Operation  Time 


REMARKS 


76-24 


t/-S7 


COUNTS  (X10  ) HOUR 


0 


1 


8 


DATA  SHEET 


DATE:  October.  1976 UNIT  NO  1024 

COMPANY:  F 

UNIT  TYPE:  Temperature 


APPLICATION:  (Type  of  Vehicle  and  location  of  Sensors) 

Lift  Truck 


CHANNEL 

— 

LEVEL 

(°F) 

TIME 

(HRS.) 

COUNTS 

. 

83 

50 

28.000 

2 

104 

4 

13,000 

3 j 

2 

4.000 

4 

144 

2 

6.000 

5 

164 

. 

2 

3.000 

6 

184 

2 

3.000 

7 

204 

2 

1.000 

8 

225 

2 

1 .000 

(Hrs.)  Total  Operation  Time 


REMARKS  76-25 


IV-59 


COUNTS  (X10J) 


0 1 2 3 4 5 6 

LEVEL 


0 1 2 3 4 5 6 7 8 

LEVEL 


\ -bZ 


DATA  SHEET 


DATE  November,  1976 unit  no 

COMPANY  D 

UNIT  TYPE:  Pressure 

APPLICATION  (Type  of  Vehicle  anil  Location  of  Sensors) 


CHANNEL 

— 

LEVEL 

TIME 

(HRS.) 

COUNTS 

I 

430 

26 

181.000 

1 

900 

2 

4,000 

3 

1380 

0 

0 

, 

4 

1850 

0 

01 

5 

2320 

0 

0 

6 

2790 

0 

0 

7 

3270 

0 

0 

8 

3740 

0 

0 

(Hrs.l  Total  Operation  Time 


REMARKS 


76-27 


DATA  SHEET 


DATE:  November.  1976  UNIT  NO  1066 

COMPANY:  [) 

UNIT  TYPE:  Pressure 

APPLICATION:  ( Type  of  Vehicle  and  Location  of  Sensors) 


CHANNEL  1 LEVEL 

(PSD 

TIME 

(HRS.) 

COUNTS 

1 Z5 

103 

269.000 

2 

440 

17 

>923.000 

3 

806 

3 

108,000 

4 

__ 

1160 

2 

6.000 

5 

1520 

1 

3,000 

6 

1890 

1 

- 

7 

2250 

■ 

0 

0 

8 

2620 

o 

0 

103  (Hrs.)  Total  Operation  Time 


REMARKS  76-28 


V-bS 


DATA  SHEET 


DATE:  November,  1976 UNIT  NO  1082 

COMPANY:  _E 

UNIT  TYPE:  Temperature 

APPLICATION:  (Type  of  Vehicle  and  location  of  Sensors) 

Loader/ Backhoe  - Pump  Inlet  Temperature 


— 

CHANNEL 

LEVEL 

TIME 

(HRS.) 

COUNTS 

1 

96 

52.000 

107 

96 

35.000 

3 

128 

92 

156.000 

4 

ISO 

76 

>356.000 

5 

172 

30 

>356.000 

6 

193 

4 

36.000 

7 

215 

3 

1 2.000 

8 

237 

0 

0 

102  (Hrs.)  Total  Operation  Time 

REMARKS  76-29 


V- 


COUNTS  (X10J)  HOURS 


AD-A0*»3  677  OKLAHOMA  STATE  UNlV  STILLWATER  FLUID  POWER  RESEARCH  —ETC  F/6  13/7 

MERADCOM/OSU  HYDRAULIC  SYSTEM  RELIABILITY  PROGRAM.  SECTION  II.  — ETC(U) 
FEB  77  DAAK02-75-C-0137 

JNCLASSIFIEO  0SU-FPRC-7m2  ML 


DATA  SHEET 


DATE:  December,  1976 UNIT  NO. : 1040 

COMPANY:  G 

UNIT  TYPE:  Temperature 

APPLICATION:  i Type  of  Vehicle  and  I /nation  of  Sensors) 

Construction  Machine  — Return  Line  Temperature 


TIME 

(HRS.) 

COUNTS 

98 

1,400 

98 

1,000 

83 

1,000 

31 

4,000 

1 

4,000 

1 

1,000 

0 

0 

0 

0 

100  (Hrs.)  Total  Operation  Time 


REMARKS: 


76-30 


V-69 


